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(material pair A) 

Figure 6-73: The influence of x 0 on the normalized stress intensity 221 
factor as a function of nondimensional time x for an 
edge crack of length /j/A^O.2 in Model II , A 2 //i 1 = 24.0 , 

R^L = 9.0 and x^/^2 =0 - 01108 » T 0 =t o D ll h \ > 

o 0 T =- a’j E 1 0 o /( 1 — Vj) (material pair A) 

Figure 6-74: The influence of x 0 on the normalized stress intensity 222 
factor k(bj) as a function of nondimensional time x for an 
edge crack of length l 1 /h 1 =0.9 in Model II , h 2 /h l = 24.0 , 

R i /L= 9.0 and x^/^2 =0 - 0n0S , x^^D^h] , 

Oq T =- cx’j £ 1 © 0 /(l — Vj) (material pair A) 

Figure 6-75: The influence of x 0 on the normalized stress intensity 223 
factor k{b^) as a function of nondimensional time x for a 
broken clad in Model II h 2 /h 1 = 3 . 0 , R i /L = 9.0 and 

X£/£ 2 =0.01108 > V* f 0^V*i » o o r =-a’ 1 £ 1 0 o /(l-v 1 ) 
(material pair A) 

Figure 6-76: The influence of x 0 on the normalized stress intensity 224 
factor k{b x ) as a function of nondimensional time x for a 
broken clad in Model II h 2 /h l = 9.0 , R i /L= 9.0 and 

X£/£ 2 = 0.01 10 8 , x 0 =r 0 D 1 /A 1 , Gq^=~ (Xj £^0q/(1— Vj) 

(material pair A) 

Figure 6-77: The influence of x 0 on the normalized stress intensity 225 
factor k(b\) as a function of nondimensional time x for a 
broken clad in Model II /i 2 //i 1 = 24.0 , RJL = 9.0 and 

X£/£ 2 =0.01 108 , x 0 =t 0 D x lh\ , o/s-aiFjOo/d-Vj) 
(material pair A) 

Figure 6-78: The influence of x 0 on the normalized stress intensity 226 
factor k(a 2 ) as a function of nondimensional time x for an 
under-clad crack of length / 2 / ^* i = 0-004 in Model II 

V A i=3. 0 , R i /L= 9.0 and x£/£ 2 =0.01108 , x 0 =t 0 D l /h \ , 

Oq T =- a’j £j 0 O /(1 -Vj) (material pair A) 

Figure 6-79: The influence of x 0 on the normalized stress intensity 227 
factor k(b 2 ) as a function of nondimensional time x for a 
fixed under-clad crack length l 2 /h^ = 0.004 in Model II 

V A i = 3- 0 , RJL= 9.0 and x£/£ 2 =0.01108 , x^^D^h] , 

Oq T =- (Xj E 1 0 o /(1 -Vj) (material pair A) 

Figure 6-80: The normalized stress intensity factor k(a 2 ) as a 228 
function of nondimensional time x for sin under-clad 
crack in Model II for x o =0.0 , h 2 /h 1 =9.0 , R i /L= 9.0 and 

X£/£ 2 = 0.01108 , x 0 =t Q D 1 /h 2 l , c 0 T =- a\ ^©o/d-Vj) 
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.(material pair A) 

Figure 6-81: The normalized stress intensity factor k(b 2 ) as a function 229 
of nondimensional time x for an under-clad crack in 
Model II for x o =0.0 , /i 2 //i 1 =9.0 , R i /L=9.0 and 

XL/£ 2 = 0.01108 , z 0 =t 0 D l /h] » °o r= “ a i E l 0 o/ C 1 _v i) 
.(material pair A) 

Figure 6-82: The normalized stress intensity factor k(a^) as a 230 
function of nondimensional time x for an under-clad 
crack in Model II for x 0 =10.0 , A 2 /Aj = 9.0 , /?;/£ = 9.0 and 

XL/E 2 =0.0im , z^D^h] , o o 7 ’=-ai£ 1 0 o /(l-v 1 ) 
.(material pair A) 

Figure 6-83: The normalized stress intensity factor k(b 2 ) as a 231 
function of nondimensional time z for an under-clad 
crack in Model II for x 0 =10.0 , h 2 /h l =9.0 , R i /L= 9.0 and 

XL/ E 2 =0.0110S , Tq=?qZ)^/Aj , Cq^=— ct'^ E^Qq/(1—\j) 

.(material pair A) 

Figure 6-84: The influence of x 0 on the normalized stress intensity 232 
factor k(a 2 ) as a function of nondimensional time z for an 
under-clad crack of length / 2 Mi = 1.0 in Model II 

/j 2 //i 1 =9.0 , R i /L= 9.0 and x£/£ 2 =0.01108 > ^o D i! h \ * 

o 0 T =- a’j E 1 0 O /(1 -V!) (material pair A) 

Figure 6-85: The influence of Zq on the normalized stress intensity 233 
factor k(b 2 ) as a function of nondimensional time z for an 
under-clad crack of length / 2 //ii = 1.0 in Model II 

h 2 /h l =9.0 , R^L -9.0 and x£/£ 2 = 0.01108 > t 0 =r 0 £ V /l l > 

Cq T =- a’j £i 0 O /(1 -V!) (material pair A) 

Figure 6-86: The influence of x 0 on the normalized stress intensity 234 
factor k(a 2 ) as a function of nondimensional time z for an 
under-clad crack of length l 2 /h ^3.0 in Model II 

h 2 /h 1 = 9.0 , R i /L = 9.0 and x£/£ 2 = 001108 > ^o=to D i/ h \ » 

g 0 t =- aj £j 0 O /(1 -V!) (material pair A) 

Figure 6-87: The influence of x 0 on the normalized stress intensity 235 
factor k(b 2 ) as a function of nondimensional time x for an 
under-clad crack of length 3.0 in Model II 

h 2 /h l =9.0 , £,•/£= 9.0 and x£/£ 2 =0.01108 , z^D^h] , 

o 0 T =- £ x © 0 / (1 - Vi) (material pair A) 

Figure 6-88: The normalized stress intensity factor k(a^) , k(b 2 ) as a 236 
function of nondimensional time x for under-clad crack 
in Model II for / 2 /A 1 = 9.0 , x o =0.0 , A 2 /ft 1 = 24.0 , RJL = 9.0 

and x£/£ 2 = 0.01 108 , z 0 =t 0 D l /h 2 l , o 0 r =- a\ ^©q/G-v^ 
.(material pair A) 
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Figure 6-89: The normalized stress intensity factor k(b 2 ) as a 237 
function of nondimensional time T for an under-clad 
crack in Model II for t Q =0.0 , h 2 /h 1 =24.0 , R i /L=9.0 and 

XL/E 2 = 0.01 108 , x 0 =t 0 D 1 /h\ , C 0 T =- a\ E x © 0 /(l - Vj) 
.(material pair A) 

Figure 6-90: The normalized stress intensity factor k(a 2 ) , k(b 2 ) as a 238 
function of nondimensional time x for an under-clad 
crack in Model II for l 2 lh^-9 . 0 , x o =20.0 , =24.0 , 

Ri/L- 9.0 and x^/^2 =0 - 01108 > *o='o D i/ A i » 

o 0 T =- a’j © 0 /(l -Vj) .(material pair A) 

Figure 6-91: The normalized stress intensity factor k(b 2 ) as a function 239 
of nondimensional time T for an under-clad crack in 
Model II for x o =20.0 , A 2 /A 1 =24.0 , £,/£= 9.0 and 

XL/£ 2 =0.01108 , x 0 =f 0 Dj/fcj , o 0 r =-ai ^/(l- v 2 ) 
.(material pair A) 

Figure 6-92: The influence of x 0 on the normalized stress intensity 240 
factor k(a£ as a function of nondimensional time x for an 
under-clad crack of length / 2 Mi= 0.5 in Model II, 

h 2 /h l = 24.0 , R i /L= 9.0 and x£/£ 2 =0.01108 » ' c 0 =f 0 Z) l/ /t l * 

c 0 T =- a’j E 1 0 O /(1 -Vj) (material pair A) 

Figure 6-93: The influence of x 0 on the normalized stress intensity 241 
factor k(b 2 ) as a function of nondimensional time x for an 
under-clad crack of length l 2 lh^= 0.5 in Model II, 

A 2 /A 1 = 24.0 , R i /L= 9.0 and %LIE 2 = 0.01108 , x 0 =f 0 ZV A i , 

g 0 t =- a’j £j © 0 /(l -Vj) (material pair A) 

Figure 6-94: The influence of x 0 on the normalized stress intensity 242 
factor kia^) as a function of nondimensional time x for an 
under-clad crack of length / 2 //i 1 =4.0 in Model II, 

h 2 /h l = 24.0 , RJL= 9.0 and xL/E 2 =0.01108 , x 0 =t 0 D 1 /h\ , 
C 0 r =- cx’^ £ 1 0 O /(1 -Vj) (material pair A) 

Figure 6-95: The influence of x 0 on the normalized stress intensity 243 
factor k(b 2 ) as a function of nondimensional time x for an 
under-clad crack of length / 2 //i 1 =4.0 in Model II, 

/^/A 1 =24.0 , Ri!L= 9.0 and x£/£ 2 = 0.01 108 > > 

Gq T =- © 0 /(l -Vj) (material pair A) 

Figure 6-96: The normalized stress intensity factor k(b 2 ) as a 244 
function of nondimensional time x for an edge crack 
crossing the interface in Model II for x Q =0.0 , h 2 /h l = 3.0 , 

RJL= 9.0 and x^/^2 =001108 » ' t o =f o I V /l i » 

g 0 t =- oc^ E x 0 O /(1 -Vj) (material pair A) 
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Figure 6-97: The normalized stress intensity factor k(b 2 ) as a 245 
function of nondimensional time x for an edge crack 
crossing the interface in Model II for x 0 =6.0 , h 2 /h l = 3. 0 , 

R i /L= 9.0 and x£AE2 = 0.01108 » x 0 =t 0 ^i/^i > 

<3q T =- a’j 0 O /(1 -Vj) (material pair A) 

Figure 6-98: The influence of x 0 on the normalized stress intensity 246 
factor k(b 2 ) as a function of nondimensional time x for an 
edge crack of length l/h l = l.2 crossing the interface in 
Model II , /i 2 //* 1 = 3.0 , R i /L = 9.0 and x^/F 2 =0.01108 , 

x 0 =t 0 D l /h l , c 0 r =- a\ Fj© q/CI-Vj) (material pair A) 

Figure 6-99: The influence of x 0 on the normalized stress intensity 247 
factor k(b 2 ) as a function of nondimensional time x for an 
edge crack of length l/h { = 2.0 crossing the interface in 
Model n , h 2 /h l =3. 0 , R i /L=9.0 and xh/£ 2 =0.01108 , 

x 0 = t Q D 1 /h 1 , o 0 T =-a’ 1 E 1 Q 0 /(l-v 1 ) (material pair A) 

Figure 6-100: The normalized stress intensity factor k(b 2 ) as a 248 
function of nondimensional time x for an edge crack 
crossing the interface in Model II for x Q =0.0 , h 2 /h 1 = 9.0 

, R i /L= 9.0 and yLIE 2 =0.0n0% > » 

o 0 r =- a’j E l 0 o /( 1 -Vj) (material pair A) 

Figure 6-101: The normalized stress intensity factor k(b 2 ) as a 249 
function of nondimensional time x for an edge crack 
crossing the interface in Model II for x 0 =10.0 , 

A 2 /A 1 = 9.0 , R i /L= 9.0 and x£/E 2 =0.01108 , , 

o 0 T =- a’j £ x 0 O /(1 -Vj) (material pair A) 

Figure 6-102: The influence of x 0 on the normalized stress intensity 250 
factor k(b 2 ) as a function of nondimensional time x for 
an edge crack of length l/h^- 1.5 crossing the interface 
in Model II , A 2 /A 1 = 9.0 , R i /L = 9.0 and x£/£ 2 =0.01108 , 
X 0 =t Q D l /h* , a 0 r =- a’j £j 0 O /(1 -Vj) (material pair A) 

Figure 6-103: The influence of x 0 on the normalized stress intensity 251 
factor k(jb 2 ) as a function of nondimensional time x for 
an edge crack of length l//t 1 =4.0 crossing the interface 
in Model II, /* 2 /A 1 =9. 0 , R i /L=9.0 and x£/£ 2 = 0.01 108 * 

x 0 =t 0 D i/ h2 i > ^o T ~~ a i E l -v i) (material pair A) 

Figure 6-104: The normalized stress intensity factor k(b 2 ) as a 252 
function of nondimensional time x for an edge crack 
crossing the interface in Model II for x o =0.0 , 

A 2 /A 1 =24.0 , £ ( /L= 9.0 and xL/E 2 =0.0U0% , x^^D^h] , 
o 0 T =- E l O 0 /(l -Vj) (material pair A) 
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Figure 6-105: The normalized stress intensity factor k(b 2 ) as a 253 
function of nondimensional time x for an edge crack 
crossing the interface in Model II for x Q =20.0 , 

A 2 Mi = 24.0 , R i /L=9.0 and xL/E 2 =0.0im , t 0 =r 0 D 1 /^ , 

o q t =- a’j E x 0 o /( 1 -Vj) (material pair A) 

Figure 6-106: The influence of x 0 on the normalized stress intensity 254 
factor k(b 2 ) as a function of nondimensional time x for 
an edge crack of length /// ^ = 1.5 crossing the interface 
in Model II , A 2 //i 1 =24.0 , R i /L= 9.0 and y,LIE 2 = 0.01108 , 

x 0 = t 0 D 1 /h l , g q t =- a’j £ : O 0 /(l -v^ (material pair A) 

Figure 6-107: The infl uence of x 0 on the normalized stress intensity 255 
factor k(b 2 ) as a function of nondimensional time x for 
an edge crack of length 0 crossing the interface 

in Model II , /i 2 /A 1 =24.0 , R i /L=9.0 and x£/F 2 =0.01108 , 

x Q =t 0 D l /h 2 1 , g 0 t =- aj ^©o/d-V!) (material pair A) 

Figure 6-108: The normalized transient temper ature distribution in 256 
Model II for x o =0.0 , h 2 /h l =9.0 , x 0 = f 0 Z? 2 , (material 
pair B) 

Figure 6-109: The normalized transient stress distribution o^/o J in 257 
Model II for x 0 =0.0 , h 2 /h 1 = 9.0 , x 0 =t 0 D 1 /h\ , 

c 0 T =- a\ E l 0 o /( 1 — Vj)* (material pair B) 

Figure 6-110: The normalized transient temperature distribution in 259 
Model II for x o =6.0 , h 2 /h l = 9.0 , x 0 =t Q D l /h 2 l , (material 
pair B) 

Figure 6-111: The normalized transient stress distribution c^/Cq 1 in 260 
Model II for x o =6.0 , h 2 /h 1 = 9.0 , i 0 =t 0 D 1 /h 2 1 , 

g 0 t =- a’j £j 0 O /(1 — Vj)- (material pair B) 

Figure 6-112: The normalized transient temperature distribution in 262 
Model II for x o =10.0 , /t 2 //i 1 = 9.0 , x 0 =t Q D l /h 2 1 , (material 
pair B) 

Figure 6-113: The normalized transient stress distribution QyyloJ in 263 
Model II for x 0 =10.0 , /i 2 //i 1 = 9.0 , x 0 =t 0 D l /h 2 , 
a\ E { 0 o /( 1 - Vj). (material pair B) 

Figure 6-114: The normalized transient temperature distribution in 265 
Model II for x o =20.0 , h 2 /h l = 9.0 , x 0 =t 0 D x /h 2 , (material 
pair B) 

Figure 6-115: The normalized transient stress distribution in 266 

Model II for x o =20.0 , h 2 /h l =9.0 , x 0 =t 0 D l /h 2 , 

Gq T =- oc’j £ ! © 0 /(l -Vj). (material pair B) 

Figure 6-116: The normahzed stress intensity factor k(bj) as a 268 
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function of nondimensional time t for an edge crack in 
Model II for x o =0.O , Aj/Aj- 9.0 , £,/£ = 9.0 and 

XL/£ 2 =0.01185 , t 0 =r 0 Z) 1 // i ; , o 0 r — aiEjGo/a-Vj) 
(material pair B) 

Figure 6-117: The normalized stress intensity factor k{b j) as a 269 
function of nondimensional time x for an edge crack in 
Model II for x o =10.0 , A 2 /^ 1 = 9.0 , R i /L = 9.0 and 

XL/£ 2 =0.01185 , x 0 =f 0 £> 1 M? , c 0 r =- a\ 

(material pair B) 

Figure 6*118: The influence of t 0 on the normalized stress intensity 270 
factor k(bj) as a fimction of nondimensional time x for 
an edge crack of length lJh^O.2 in Model II , 

V/t!=9. 0 , /?,/£= 9.0 and x^/E 2 =0.01H5 , x 0 =r 0 D 1 /Aj , 

Oq T =- cx’j E l 0q/(1 -Vj) (material pair B) 

Figure 6-119: The influence of x 0 on the normalized stress intensity 271 
factor k(b j) as a function of nondimensional time x for a 
broken clad in Model II for P L = 0.552538 , h 2 /h i =9.0 , 

R i fL=9.0 and x^/^2 =0 - 01185 > ^o=h D J h ] > 

c 0 T =- a’j £j 0 O /(1 -Vj) (material pair B) 

Figure 6-120: The normalized stress intensity factor k(a 2 ) as a 272 
function of nondimensional time x for an under-clad 
crack in Model II for 02=0.4512416 , x o =0.O , h 2 /h 1 = 9.Q , 

R i /L=9.0 and x^/^ 0 - 01185 » x o =f O £ V /l i . 

o 0 T =- ct\ E 1 Qq/(1-\ 1 ) .(material pair B) 

Figure 6-121: The normalized stress intensity factor k(b 2 ) as a 273 
function of nondimensional time x for an under-clad 
crack in Model II for 02=0.4512416 , x 0 =0.0 , /i 2 //z 1 = 9.0 , 

/?,/£= 9.0 and x^/^2 =0 - 01185 > x o =t o D \! h \ > 

a 0 r =- a’j E j 0 q/(1 -Vj) .(material pair B) 

Figure 6-122: The normalized stress intensity factor k(a 2 ) as a 274 
function of nondimensional time x for an under-clad 
crack in Model II for 02=0.4512416 , x 0 =10.0 , h 2 /h i =9.0 

, RJL = 9.0 and x^/^2 =0 - 01185 > J h l > 

o/=- o’j E l O 0 /(l — Vj) .(material pair B) 

Figure 6-123: The normalized stress intensity factor k(b 2 ) as a 275 
function of nondimensional time x for an under-clad 
crack in Model II for 02 = 0.4512416 , x 0 =10.0 , A 2 //i 1 = 9.0 

, R i /L= 9.0 and x^/^2= 001185 » T 0 =f 0 Z V A i » 

Gq T =- o’j £ 1 0q/(1 -v^ .(material pair B) 

Figure 6-124: The influence of x 0 on the normalized stress intensity 276 
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factor k(a 2 ) as a function of nondimensional time x for 
an under-clad crack of length / 2 //t 1 = 1.0 in Model II, 
02 = 0.4512416 /i 2 /* 1 =9.0 , R i /L= 9.0 and x£/£ 2 =0.01185 * 

x 0 =t Q D 1 /h 2 l , O 0 r =- aj E x 0 o /( 1 -VjXmaterial pair B) 
Figure 6-125: The influence of x 0 on the normalized stress intensity 
factor ifc(& 2 ) as a function of nondimensional time x for 
an under-clad crack of length / 2 /A 1 = 1.0 in Model II 
02=0.4512416 /i 2 //* 1 =9.0 , /?,/£= 9.0 and xL/£ 2 =0.01185 , 

x 0 =<o^i/^i > o 0 r =- aj £ x 0 O /(1 -v ^(material pair B) 
Figure 6-126: The normalized stress intensity factor Jfc x as a function 
of nondimensional time x for an edge crack crossing the 
interface in Model II for p 1 =a 2 =0.0187223 ,x o =0.0 , 

V*i=9- 0 , R t /L= 9.0 and x£/F 2 =0.01185 , x^^D^h] , 

O 0 r =- aj £j 0 q/(1 -Vj) . (material pair B) 

Figure 6-127: The normalized stress intensity factor as a function 
of nondimensional time x for an edge crack crossing the 
interface in Model II for Pi= 02=0.0187223 , x o =0.0 , 

A 2 /A 1 =9 .0 , RJL= 9.0 and x£/£ 2 =0.01185 > t 0 =f 0 £) l/ /l i > 

Q(f=- aj 0 O /(1 -Vj) . (material pair B) 

Figure 6-128: The normalized stress intensity factor k(b 2 ) as a 
function of nondimensional time x for an edge crack 
crossing the interface in Model II for P 1 = 0^=0.0187223 , 
x 0 =0.0 , h 2 /h l =9.0 , R i /L= 9.0 and x E / E 2 =0 - 0m5 > 

x 0 =t 0 D l /h 1 , g 0 t =- aj E x 0 O /(1 -Vj) . (material pair B) 
Figure 6-129: The normalized stress intensity factor k x as a function 
of nondimensional time x for an edge crack crossing the 
interface in Model II for Pi =0^=0.0187223 , x o =10.0 , 

h 1 /h l =9.0 , R^L =9.0 and x£/£ 2 =0.01185 , x^^Dy/h] , 

Gq^=— aj E j 0 O /(1 — Vj) . (material pair B) 

Figure 6-130: The normalized stress intensity factor k^ as a function 
of nondimensional time x for an edge crack crossing the 
interface in Model II for p i =ot 2 =0.0187223 , x o =10.0 , 

A 2 /A 1 =9.0 , R i /L= 9.0 and x£/£ 2 =0.01185 , x 0 =t 0 D 1 // , 

oJ=- aj E l © 0 /(l -Vj) . (material pair B) 

Figure 6-131: The normalized stress intensity factor k(b 2 ) as a 
function of nondimensional time x for an edge crack 
crossing the interface in Model II for p i =ot 2 =0.0187223 , 
x o =10.0 , h 2 /h { = 9.0 , R i /L= 9.0 and x E / E 2 =0 - 01185 , 

x 0 =t Q D l /h l , o 0 r =- aj Ey 0 o /(1 -Vj). (material pair B) 
Figure 6-132: The influence of x 0 on the normalized stress intensity 
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factor k x as a function of nondimensional time x for an 
edge crack of length //A 1 = 1.5 crossing the interface, in 
Model H for ^=02 = 0.01 872238 , ^^=9.0 , R t /L = 9.0 

and x£/£ 2 =0.01185 , T Q =t 0 D 1 /h 2 1 , o o r —aj£ 1 0 o /(l-v 1 ) 

. (material pair B) 

Figure 6-133: The influence of x 0 on the normalized stress intensity 285 
factor Jfc as a function of nondimensional time x for an 
edge crack of length llh x = 1.5 crossing the interface , in 
Model H for p 1= 02 = 0 . 01 872238 , A 2 /A 1 = 9.0 , /?,/L= 9.0 

and x l / £ 2 =0 - 01185 » x 0 =<o £) i/*i » 

Cq T =- a’j ^©q/CI -V j) . (material pair B) 

Figure 6-134: The influence of x 0 on the normalized stress intensity 286 
factor k(b 2 ) as a function of nondimensional time x for 
an edge crack of length l/h^ = 1.5 crossing the interface , 
in Model II for p 1= 0^=0.01872238 , h 2 /h 1 = 9.0 , R i /L = 9.0 

and xL/E 2 =0.011&5 , x Q =t Q D l /h 2 1 , 

a/=- a’j E l 0 o /( 1 -Vj) . (material pair B) 
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Table 6-14: The maximum of normalized stress intensity factors for 133 
an under-clad crack subjected to transient thermal 
stresses, h 2 /h l = 24.0, Rj/L=9.0 and xL/E 2 =0.01108 . 
(material pair A). 

Table 6-15: The maximum of normalized stress intensity factors for 134 
an edge crack subjected to transient thermal stresses, 
h^lhy^'i. 0 , R^L-9.0 and %h/£ 2 = 0.01 108 . (material pair 

A). 

Table 6-16: The maximum of normalized stress intensity factors for 135 
an edge crack subjected to transient thermal stresses, 
h^!h^= 9.0, R i /L= 9.0 and xL/£ 2 = 0.01108 . (material pair 
A). 

Table 6-17: The maximum of normalized stress intensity factors for 136 
an edge crack subjected to transient thermal stresses, 

/i 2 //i 1 = 24.0, R i /L=9.0 and %Z./£ 2 = 0.01 108 . (material pair 
A). 
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ABSTRACT 


In many applications, such as nuclear pressure vessels, micro-electronics, 
and the chemical industry, the study of transient thermal stresses in layered 
cylinders and plates in the presence of a crack is required. In this study first the 
cladded pressure vessel under thermal shock conditions which is simulated by 
using two simpler models. The first model (Model I) assumes that, if the crack 
size is very small compared to the vessel thickness, the problem can be treated 
as a semi-infinite elastic medium bonded to a very thin layer of different 
material. However, if the crack size is of the same order as the vessel thickness, 
the curvature effects may not be negligible. In this case it is assumed that the 
relatively thin walled hollow cylinder with cladding can be treated as a 
composite beam on an elastic foundation (Model II). In both models, the effect of 
surface cooling rate is studied by assuming the temperature boundary condition 
to be a ramp function. Among the crack geometries considered are: the edge 
crack in the clad, the broken clad, the edge crack going through the interface, 
the under-clad crack in the base material, and an internal crack crossing the 
interface. 

The calculated results include the transient temperature, thermal stresses 
in the uncracked medium and stress intensity factors which are presented as a 
function of time, and the duration of cooling ramp . The stress intensity factors 
are also presented as a function of the size and the location of the crack. The 
problem is solved for two bonded materials of different thermal and mechanical 
properties. The mathematical formulation results in two singular integral 
equations which are solved numerically. The results are given for two material 
pairs, namely an austenitic steel layer welded on a ferritic steel substrate, and a 
ceramic coating on ferritic steel. In the case of the yielded clad, the stress 
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intensity factors for a crack under the clad are determined by using a plastic 
strip model and are compared with elastic clad results. When the results 
obtained are compared with the corresponding elasticity solution for the thick- 
walled cylinder, the agreement was found to be quite satisfactory. 

In the composite plate with finite thickness the calculated results indicate 
that the peak values of the stress intensity factors decrease with decreasing 
cooling rate of the surface. On the other hand, in the case of semi-infinite 
medium (Model I), since the medium is fully constrained at infinity, the 
maximum values of the stress intensity factors are attained at steady state. 

The technique developed and the dimensionless results given in this study 
are , of course, applicable to any composite plate that is locally subjected to 
thermal shock. 
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Chapter 1 
INTRODUCTION 

It is known that the cladding, when bonded to the inner wall of pressure 
vessels or pipes, can be very effective in protecting the base metal from severe 
corrosion. In the nuclear pressure vessels the cladding serves the additional 
purpose of protecting the base metal from radiation damage. In some other 
applications, as in micro-electronics and the chemical industry where coating is 
used very widely, the study of residual stresses and transient thermal stresses 
in layered cylinders and plates is required. 

Our aim in this work is to study the structural integrity of cladded 
materials under severe thermal transient stresses in the presence of pre- 
existing flaws. These flaws may be due to imperfections such as voids, 
inclusions, or weld defects, which are generally treated as cracks. In case of 
pressure vessels with cladding, the cracks have been observed underneath the 
clad oriented in a plane perpendicular to the cylinder axis and terminating at 
the interface. 

In analyzing the subcritical crack growth in homogeneous materials, it is 
generally accepted that the stress intensity factor can be used quite effectively 
as a correlation parameter. The objective of this study is therefore to investigate 
the effect of the clad on the stress intensity factors in the case of a crack 
perpendicular to the interface under transient thermal stresses. Once we know 
the stress intensity factor, we could also determine whether catastrophic failure 
will occur in brittle materials due to unstable crack propagation. 

The thermal shock problem for a circumferentially cracked hollow cylinder 
with cladding was considered by Nied [1] for two special crack configurations. 
First, the circumferential crack embedded in the base metal, where the crack tip 
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is located just underneath the clad. Second, the edge crack which completely 

passes through the clad. His basic assumptions were 

1. The problem is axisymmetric for both crack geometry and 
temperature boundary conditions. 

2. The change in surface temperature is a unit step function. 

3. The transient thermal stress is quasi-static, i.e, the inertia effects 
are negligible. 

4. The materials are assumed to be isotropic, homogeneous, and 
linearly elastic such that the clad and the base material have the 
same elastic modulus and Poisson’s ratio, i.e, =£ 2 , v i =v 2 > 
non-homogeneous in thermal properties. 

5. If the cladding is sufficiently ductile, during a thermal transient, 
the cladding may yield. In this case the clad is examined by using a 
plastic strip model, which assumes that the clad is perfectly 
plastic. 

Now, by assuming an extremely severe temperature change at the inner 
wall of the cylinder, i.e., a step change in the inner wall temperature, he 
obtained a peak value of thermal stresses. These stresses should be lower if one 
assumes a more realistic temperature boundary condition. Erdogan [2] 
compared two sets of data to show the effect of the temperature boundary 
condition at the inner wall of a vessel without cladding, first set for a 
homogeneous vessel which contains a circumferential edge crack on the inner 
wall obtained from [3] where a step change in temperature is used as the 
boundary condition, and the second set for the same problem where a smoother 
temperature boundary condition obtained from [4] was used. He showed a 
si gnifi cant difference between the two sets. Hence, a reduction in the stress 
intensity factors in Nied’s results would occur if we repeat his work with more 
realistic temperature boundary condition. 

The main objective of this work is to solve the crack problem under 
transient thermal stresses with temperature boundary conditions at the inner 
wall of the cylinder that are somewhat more realistic than the step function 
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temperature used in previous studies as shown in Figure 1.1. c , and provided 
that the clad and the base materials have different thermal and mechanical 
properties. 

The actual problem is a very complicated three-dimensional problem and 
seems to be analytically intractable. So, we are going to simulate the problem by 
using two simpler models depending on the crack size. 

Model I 

Since the clad has a relatively small thickness, the base material has a 
large bulk, and the radius of the cylinder is large, the total strains in axial and 
circumferential directions would be nearly zero during transient heating or 
cooling. Also, since the surfaces of the cylinder are stress free, the radial 
stresses would also be negligible small. Thus, in case of a crack that is very 
small in size compared to the other dimensions, the problem can be treated as a 
semi-infinite elastic medium (base) bonded to a surface layer of different 
material (clad) of thickness h, shown in Figure 1.1. a (see for example [2]). In 
some other applications, this model can also be applied to certain micro- 
electronic devices and ceramic coated metal parts consisting of a relatively thin 
coating bonded to an elastic substrate and subjected to rapidly changing 
thermal environments. 

Model H 

In the case of a relatively large crack size, the more realistic model for 
thin walled hollow cylinders with cladding would be a composite beam on an 
elastic foundation in which the modulus of the foundation x is a function of the 
thickness of the cylinder wall, the radius of the cylinder, and the modulus of 
elasticity, see [5], The problem of interest is depicted in Figure l.l.b. This 
model can also be used in other cases such as, for example, a fully or partially 
constrained composite plate in which the effect of the constraints can be 
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Figure 1-1: a-The geometry for Model I, b-The geometry for Model II, c-The 
temperature boundary condition at the inner wall. 

represented by an elastic foundation. 

In recent years the results of various studies on the fracture mechanics of 
layered composites appeared in literature. These studies were concerned with 
interface cracks, cracks perpendicular to the interface, cracks terminating at the 
interface, and cracks crossing the interface. The singular behavior of stresses 
around a crack at the interface of two elastic half-planes of dissimilar material 
are examined by Williams [6]. Rice and Sih [7] considered a finite crack at an 
interface of two joined materials subjected to both symmetric and skew- 
symmetric in-plane loading. Also, the interface crack problem under "bending” 
was reported by Sih and Rice [8]. Stress distribution of two half-planes bonded 
to each other containing cracks along the bond was considered by England [9] 
and Erdogan [10]. The interface crack problem in multi-layered material has 
also been studied by Erdogan and Gupta [11]. Stress distribution in case of 
penny-shaped cracks lying between two bonded half-spaces are reported by 
Lowengrub and Sneddon [12], Lowengrub [13], and Willes [14]. A more 
generalized situation corresponds to the case of interface cracks lying on a 
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circular arc was considered by Perlman and Sih [15 ,16]. 

Khrapkov [17], and Cook and Erdogan [18] have considered the problem of 
a crack normal to a bimaterial interface. The case of a laminated composite 
containing a broken lamina was investigated by Ashbaugh [19] and Gupta [20]. 
Erdogan [21] has investigated the singular nature of the crack tip stress field in 
a non-homogeneous medium having a shear modulus with a discontinuous 
derivative for the antiplane shear loading of two bonded half spaces in which 
the crack is perpendicular to the interface. The fracture problem of a single 
layer of dissimilar material with a crack normal to the interface between two 
other layers of infinite height was considered by Hilton and Sih [22], Bogy [23], 
and Arin [24], while a crack parallel to the interface was considered by Hilton 
and Sih [25]. The anti-plane shear case was treated by Chen and Sih [26]. The 
same problem with penny-shaped crack under normal extension was solved by 
Arin and Erdogan [27] and torsion by Sih and Chen [28]. 

A crack going through the interface in two bonded half planes was 
investigated by Erdogan and Biricikoglu [29]. Erdogan and Cook [30] considered 
antiplane shear cracks terminating at and going through the interface. The 
fracture problem of a composite plate which consists of a bonded parallel load 
carrying laminates and buffer strips, and its limiting case of the collinear cracks 
joi ning and fo rming a stress-free end have been considered by Erdogan and 
Bakioglu [31 , 32]. Also, Delale and Erdogan [33] considered the same problem 
but for orthotropic materials. Goree and Venezia [34] have investigated the 
bonded elastic half-planes with an interface crack and a crack perpendicular to 
the interface. The problem of two bonded semi-infinite elastic media with a 
crack in one of the half-planes lying parallel to, and at an arbitrary distance 
from the interface has been investigated by Erdogan [34]. A summary and 
discussion of the various modes of cracks in composite material has been 
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reported by Erdogan [36]. 

The problem considered can be solved by the superposition technique as 
shown in Figure 1.2. Once the temperature distribution is known, the thermal 
stresses for uncracked problem can be obtained. The perturbation problem 
(mixed boundary value problem) can then be formulated by using the thermal 
stresses from uncracked problem, but with opposite sign, as the crack surface 
traction along the line of the crack. 


thermal stresses 


thermal stresses 


a T 


+ 

O 

r 


uncracked problem perturbation problem 


Figure 1-2: The superposition 

First, the problem of cracks fully embedded into the homogeneous strips, 
perpendicular to the interface is considered. A general formulation of the 
perturbation problem is given for plane strain and generalized plane stress 
cases by the use of Fourier Integral Transform. The singular behavior of the 
stresses for different crack geometries is studied in some detail by using the 
Muskhelishvili Technique [37]. In cases where the crack terminates at the 
interface, or goes through the interface, the stress singularity is a function of 
material constants as well as the intersection angle. 

The resulting system of singular integral equations is solved numerically 
by using the technique described in [38 , 39]. The stress intensity factors are 
calculated as functions of time and dimensions concerning the size and location 
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of the crack for various crack geometries and for two materials combinations 
described in Table 6-1. Also, in order to study the influence of the thermal 
boundary conditions on the stress intensity factor, the increase in time of the 
surface temperature is assumed to be an additional variable. 
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Chapter 2 

ANALYSIS OF THE PROBLEM 


2.1 Model I- Semi-infinite elastic medium (base metal) 
bonded to a surface layer of different material (clad) 

2.1.1 Temperature distribution 

2.1.1.1 Unit step function at the boundary 



Figure 2-1: a-Geometry of the problem (Model I) b-Temperature boundary 
condition (unit step function) 

Figure 2-l(a) represents the problem of interest, where x> is measured 
from the interface of the two materials, and Figure 2-l(b) represents the 
temperature change at the boundary. Let Dy , k'y and jD 2 > ^2 3X6 ^ ierma ^ 
diffusivity and thermal conductivity for the clad (1) and base (2), respectively. 


So, the differential equations for materials 1, 2 are : 


cP-Ty(x',i) ~ 1 dTy(x\t ) 
dx ’ 2 D \ dt 


-h<X’<0 ,t>0 , 


(a) 


3 2 r 2 (xM)_ 1 BT 2 (r,t) 
dx ’ 2 £*2 dt 


0<X’<oo ,t > 0 . 


( 2 . 1 ) 

(b) 
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where T 2 and T 2 are the temperatures in materials 1, 2 respectively. The initial 
and boundary conditions are 


o 

li 

O 

II 

8^ 

(a) 

T l (0,t) = T 2 m, 

(b) 

ar : ( o,t) dr 2 (o.o 

** dr "** dr • 

(c) 

lim T 2 (x’,t) = T oe , 

V—^oo 

(d) 

T l (-h,t) = T ee +(T Q -TJH(t). 

(e) 


where is the initial temperature for both materials 1, 2, and T 0 is the 
temperature at the boundary at any time t > 0, and H(t) is the Heaviside step 
function. Condition (2.2)(c) implies perfect heat transfer at the interface 0c=0). 
Let, 


0 1 (j:’,O=r 1 (A:’,o-r oo , 

(a) 

Q 2 (x\t) = T 2 (x',t)-T 00 , 

(b) 

®0 =To~ T ~- 

(c) 


Then, the two differential equations (2.1) can be written as 


dx ^ D, dt 


di 2, & 2 (x , ,t) i dQ 2 (x%t) 


, -h < x y < 0 , t > 0 , 


, 0<;t’<°°,r>0. 


dx ^ D 2 

Also, the initial and boundary conditions (2.2) become 


(a) 


(b) 


(2.4) 
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Q 1 (r,O) = 0 2 (x\O) = O, 

(a) 


©i(0,0 = © 2 (0,r) , 

(b) 


a©!^) a© 2 (o^) 

*i - k 'i -s 

ax' ax’ 

(c) 

(2.5) 

Hm © 2 (jcv) = 0 , 

X* — > oo 

(d) 


0 1 (-/r^) = 0 o //(r). 

(e) 


By applying Laplace Transform to the partial differential equations (2.4) with 
respect to t, and by using condition (2.5Xa), they may be reduced to two ordinary 

differential equations of the form 



cfie^x’p) 2 

^2 = ?1 ©i(*’,p), 

(a) 

(2.6) 

d 2 ©^) 2 

— = <?2 ©2^’^)- 

(b) 


where 



1 ^IcT 

-r 

II 

NcT 

II 


(2.7) 

The integral transform is defined by 



<=Kx\p)= f Q(x’,t)e~P‘ dt. 

J 0 


(2.8) 

The general solutions of (2.6) are 



©^Jt’p) = A x (p) coshqpc' +B l (p)sinhq l x < , 

(a) 

(2.9) 

0 2 (x'p) = A 2 (p ) + B 2 (p) . 

(b) 



where, A-^p), B^ip), A 2 (p), B 2 (p) are unknown functions to be determined from 
conditions (2.5)(b-e). By applying Laplace Transform to the conditions (2.5)(b-e) 
we obtain 
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©i(0 ,p) = 0 2 (°’P) ' 

dO : (Op) d © 2(0 ,p) 


lim Q 2 (x'f>) = 0, 

r—>oo 


(d) (2.10) 


Q l (-hj>) = — . (e) 

P 

After substitution of the conditions (2.10)(b-e), and determination of the 

constants Ayip), By(p), A 2 (p), B 2 (p), the equations (2.9) become 
_ Q 0 [ cosh qyX’-r\ sinh q^X’] 


©l (x'p) = 


p [cosh q j/t+T) sinh q^h] 


( 2 . 11 ) 


&2(x'p) = 


p [ cosh q^h+T) sinh q^h ] 


where 


*2 ,°i 
rl ‘-' l D- 2 
*1 


( 2 . 12 ) 


The temperature distribution Qyfat) , 0 2 (x’,t) can be obtained by applying 
the Inversion Theorem of Laplace Transform on equations (2.11), i.e, 

l r-y+i°o 0 q [ cosht^^x' — T| sinh l^yX’ ] 


l ry+i°° [ cosn^yX' — t| sinn j 
®l(*> f ) 2 jc * 2 [ coshC,yh+T\ sinh^yh ] 


1 fY+«“ 


[ cosh t^yh +q sinh^yh ] 


(2.13) 


where 


Cl "V z/D 1 , ^ 2 = ^l z/D 2 . 


(2.14) 


Since, the functions (e ^ Z ! D 2 *) and (sinh'l z/Dyh) are double-valued 
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functions when the angle in the complex domain is changed from 0 to 2rc, the 
integrands in (2.13) have a branch point at z = 0, and the integration contour 
should be as shown in Figure 2-2. 



Figure 2-2: Contour for evaluating the integrals in eqns. (2.13) 

So, applying the Residue Theorem in equations (2.13), and using the 

contour integral T shown in Figure 2-2, we obtain the following closed form 

solutions for ©jOc’.r) , 0 2 (x’,O 

2t| f 00 sin tfjc' + h) d% 

©0 71 ^0 [cos?’£ ) h+T[ 2 sin 2 ty] ^ 

©2(*’>0_ j _2 C^e^&lcostysin^Sx’+TisintihcosZfix’ldZ, „ (2 15) 

©0 ^0 [cos^fyi+vfsinpty] 4 

5 = ^£> 1 /£> 2 . (c) 

Defining now the Fourier Number x and the dimensionless parameters by 
x = tD x /h 2 ,Z 3 ' = h%,x* = (p+h)lh . (2.16) 

equations (2.15) become 
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;0£;t*<l, 


(a) 


®i(* » T )_ j 2rj f°° e ^ sintf'X* d^‘ 

©0 K ■'0 [cos 2 ^'+r\ 2 sin 2 ^’] 

(2.17) 

©2( ;c *> T )_^ 2 f °° g~ T ^' 2 [cog wt ^’5(x* - 1 ) + 1| sin %■ cos £,'8(x* -l)]d^ 

©o n ^0 [coP , 't ) '+T\ 2 siTp-Q] 

; l£r*£°o. (b) 

2.1. 1.2 Ramp function at the boundary 




Figure 2-3: a-Geometry of the problem (Model I) b-Temperature boundary 
condition (ramp function) 

Figure 2-3(a) represents the problem of interest, and Figure 2-3(b) 
represents the temperature change at the boundary. The differential equations 


for 1, 2 are: 

a 2 ©!(*•,*) 

a *- 2 


1 a ©!(*•,/) 

D[ dt 


;-h<X’<0,t>0, 


(a) 


a 2 0 2 (r-,f) i de 2 (x\t) 


dr- 2 


dt 


;0<r’<t»,f>0. 


The initial and the boundary conditions are: 


(2.18) 

(b) 
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0 1 (jr,O) = 0 2 (^’,O) = O, 

(a) 

©i(0,r) = © 2 ( 0 , 0 ) 

(b) 

w 30j(O,O ff 9 0 2 (O,O 
1 3*’ 2 dx- 

(c) (2.19) 

lim 0 2 (x’,O = 0 , 

JP-^oo 

(d) 

©0 ©0 

©!(-/!, 0 = JitHV—it-tJHit-tQ). 
r o r o 

(e) 

As before, by applying Laplace Transform to equations (2.18), and using 

the condition (2.19)(a), we will obtain two ordinary differential equations having 

the general solution of the form 


0^-p) =A x (p) cosh q { x< +B { (p) sinhq-pc , 

(a) 


(2.20) 

^(x-p) = A 2 (p) e*2*+B 2 rW . 

(b) 

where Aj(p), B^ip), A 2 (p), B 2 (p), are unknown functions to be determined from 

the conditions (2.19)(b-e). 


Now, applying Laplace Transform to the boundary conditions, we may 

have 


©l(0p) = © 2 (Op) , 

(b) 

d©i( Op) d0 2 (Op) 

if, ~k\ , 

1 dx> 2 dx< 

(c) 

lim © 2 (rp) = 0, 

jr — >«« 

(d) (2.21) 

0rv 1 0 a 1 

r o p b p 

(e) 

Then, by substituting from equation (2.20) 

into equation (2.21), and 


16 



solving for the constants A x (p), B^ip), A 2 (p), £ 2 (p), we find 


where 


® l (x'p)=F n (p) .B n (p), 

(a) 

@2(x'p) = F l2 (p) .B n (p). 

(b) 

0 O [ coshq^X' -q sink q±r ] 
plcoshqfi+Tisinhqfi] 

(a) 

©og^ 

F,~(p) - 

LZ plcoshq^h+^sinhq^h] 

(b) 

1 1 e~P‘o 

B n (p)=B n (p) = -(-—-). 

*0 P P 

(c) 


It is clear that, F 11( F 12 , are the same as equations (2.11). 

By applying the convolution theorem to equations (2.22), we find 

; O' =1.2) . 


where, 

/ 11 (0=L- 1 [F u (p)] 

= © ri-?H r e ~ tD ^ 2 sin ^ x ’ +h ^ 

n Jo [cos^A+q 2 sin 2 ^h] ^ 

f 12 (t)=L-HF 12 (p)] 

-0 n 2 r°°e~ lD fi 1 [cost l hsinl l dx’+T[sin^icos%bx’] dl 
n Jo [co^^A+q 2 sin 2 fyi] ^ 

b n (t) = A 12 (0 = IT 1 [G u (p)] = i[ 1 -//(r-r 0 )] . 

f 0 

From (2.16) , (2.24) and (2.25) it then follows that 
Q i(**> x ) = x^ + 2r\ r° °(g~ T fr 2 -l) + 

0 O T 0 K Jo X 0 


(a) 


(b) 

(c) 


sin E,rx* dfy 

[cop- ^+q 2 im 2 ^’] {j> 3 


;0£jc*£1 ,t<t 0 , (a) 


( 2 . 22 ) 


(2.23) 


(2.24) 


(2.25) 
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0 i ( ** ,T) _ 1 2 t| r r^(g T o^-i) t 

©0 JI Jq t 0 


sin £,vc* d%' 

[ COS 1 ^4- T| 2 ^’] ^’ 3 


; 0^a:*< 1 ,x>x 0 , (b) 


0 2 (.x*,t) _ T t 2 
©0 Xq n •'0 x 0 


(2.26) 


( cos &(x* — 1 ) + T) sin ^cos ^(x* - 1 )) dt; 

[ cos 2 T) 2 sin 2 ^’] 4’ 3 

1 < X* < °° , X<X 0 , (c) 

0 2 (a:* ) x) _ ^2 f O °e-'^’ 2 (e' c o^- 1 ) # 

©0 71 Jo x 0 


where 


( cos Ij’sm - 1 ) + r| sin %cos ^-6(x* - 1 ) ) dt; 

[ cos 2 £ 4 - r\ 2 sin 2 £•] ^’ 3 

1 < jc*<°° , x> x 0 . (d) 


x 0 = 


«6*1 *2 I & i D 1 

_, n . V_.6=V_ 

*1 


(2.27) 


2.1.2 Thermal stresses in uncracked problem 

Consider the plane strain problem in y and z directions as shown in 
Figure 2-4, i.e, 


e «y “ e iyz - e izx “ £ iyy “ 6 izz “ ® * (* 1 > 2) . (2.28) 

Since the temperatures 0- , (/ = 1, 2) are functions of (jc , t), then the 
stresses Eire functions of (x , t ) only. It is obvious that all the shear stresses are 
zero. Then, the equilibrium equation for each medium (ignoring the inertia 
effects) is : 
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Figure 2-4: Geometry of Model I (uncracked problem) 
da^ (x,t) 

- 0. 

dx 


(2.29) 

Since, we have a stress-free boundary , i.e, 0^(0 , t) 
= 0 . 

= 0 , then, 

(2.30) 

So, the Hooke’s Law becomes 



-a’0 = i[a yy - v o zz ]. 

(a) 


> 

1 

K 

£ 

-Tuj 

ii 

0 
8 

1 

(b) 

(2.31) 

where a* is the coefficient of thermal expansion, E is the Young’ s modulus, and v 

is the Poisson’s ratio of the material. From equations (2.31) we may have 

~ *-r 0C’£ 0 

a yy~ a zz~ (i_v>- 

(2.32) 

Then, the thermal stresses in materials 1, 2 are: 



a i £ l Q l(*>0 

G\yy(Xj) ( 1 — V ) * 0 

(a) 


cc 2 E 2 e 2 (x,t) 

OlyyiXl) ~ <*2 zz( X ’0 ~ j = **** 00 • 

(b) 

(2.33) 
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Let, 


a i ^1 ®o * _ x 
(1-Vi) X ~ h T " h 1 


Then, equations (2.33) become 


CTijyC^.T) _ 0 i(j:*,t) 


©n 


; 1 , 


°2yy(**’^ (X 2 ^2 0 —v l) 

“£[(1-v 2 ) 0 o 
U 1 




(2.34) 


(a) 


(2.35) 


(b) 


2.1.3 Formulation of the crack problem (perturbation problem) 


y ,k 



Figure 2-5: Crack geometry (Model I) 

Consider the two-dimensional composite medium shown in Figure 2-5, 
where each material contains a finite crack perpendicular to the interface. 
Because of symmetry, the problem will be considered for, 0 £ y < <». It can be 
shown that, the governing differential equations for the displacements in each 
medium are (Appendix A) 
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, 1 X V 72 ^ v 

(k- 1 ) V 2 w + 2 — 

dx 2 ax ay 

, ,.^72 - , d 2 !/ d 2 U 

(k-1)V 2 u+2(— — +— — 

oxdj dy* 


II 

o 

(a) 

11 

o 

(b) 


(2.36) 


where u and v are x and y components of the displacement vector, respectively. 
k=(3-4v) for plain strain, and k=(3-v)/(1 +v) for generalized plane stress, v is 
Poisson’s ratio. 

Assume that the solution of equations (2.36) may be expressed in terms of 
the following Fourier Integral, i.e, 

u(x,y) = - f f(x,ct)cosya.da+^- [ A(y,P)e“Pdp, (a) 

JtJ 0 2 tcJ „ 

2 


(2.37) 


v(x,y) = - f g(x,a)sinyada+^~ f . 

nj 0 27tJ_ 0< , 


(b) 


Because u, v, are symmetric and antisymmetric with respect to y, respectively, 
the first term of u is a Cosine transform, and the first term of v is a Sine 
transform. The second terms in equations (2.37) represent the general Fourier 
transform. Substituting equations (2.37) into equations (2.36) we obtain a 
system of ordinary differential equations for the unknown functions f , h , g , k . 


^ K x >°0 - 2 a 2 fix, a ) + a 4 f(x,a) = 0, 

(a) 

h(y$)-2 p 2 |i A(y,p) + P 4 /i(y,P) = 0 , 
dy* djr 

(b) 

j2 

— g(x, a) -2a 2 — g(x,a) + a 4 g(x, a) = 0 , 

(c) 

*(y,p)-2 p 2 -^ *(y,P) + p 4 *0>,p) = 0 . 
dy* dy 

(d) 


(2.38) 


The general solution of these four differential equations are: 
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(a) 


Ax, a) = (q + C 2 x)e xa + (C 3 + C 4 x) e *° , 

g(x, a) = (q + q*) e*° + (q + Qc) «-*“ , (b ) 

(2.39) 

fcO.P) = (D l +D^)e^+(D 3 +D^y)e~^, (c) 

A(y,P) = {D\+D' 2 y)eyW+(D' 3 +D'0)e-yM. (d) 

where 0<a<~ , -«< p< «> t o. C t - and C’, (i = l,...,4) are functions of the 
transform variable a and are linearly dependent. Similarly, D i and (/=1 ,...,4) 
are functions of the transform variable P and are not independent. 

Since, u,u are bounded as y from (2.39) it follows that D v Z> 2 » D \, 
must be zero. Thus, after eliminating C- and D'- (/=l,...,4;/=3,4) by using 
coupling relations, equations (2.39) may be written as 


fact) = (C l +C 2 x)e xa +(C i +C 4 x)e~ xa , 

(a) 


g(x,a) = (-C 1 -^C 2 -C 2 x)e xa + (C 3 -^C 4 +C 4 x)e~ xa , 

(b) 


*(y,p) = (£> 3 +D 4 y)c-ylPI, 

(c) 

(2.40) 

%,P) = ^ (^3 " ^ P 1 • 

(d) 


Defining, D 4 (P) = A(P), the equations (2.40) (c,d) can be written in the form 

^’ P)= 7^ ( 5) IPI+IPI},) ^ IP1, 

(c) 


*0'.P) = =^^(^IPI-K+|p|y)c->'lPl . 

(d) 

(2.41) 

Hence the displacements u , v are in the form 



u(x,y) = [ (q + Ctf) + (C 3 + C A x) e - * 0, ] cos yet da 

(a) 
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(2.42) 


y(*oO=-P [(-C 1 -C 2 { ±+x })^+(C 3 -C 4 { ±-x })e-^] 

TCJ q OC CX 

^(^iPl + lPly^Pl^dp. (b) 


The stress-strain relations for plane-strain problem are 
2(1 — v) v 

a ~ = ^(T^ (e « l-v^’ 

2(1 -v). v 

°yy ^ (i-2v> ^ e w ,+ i-v e;w ^’ 


(a) 

(b) (2.43) 


where p = £/2(l+v) is the shear modulus. 


Also, the strain-displacement relations sire 


du 

dx’ 


(c) 


(a) 


dv 

e ^ = a? 

du dv 

y * = Ty + Tx- 


(b) (2.44) 

(c) 


By substituting equations (2.42) into equations (2.43) and (2.44), and 


observing that at the plane of symmetry y = 0 the shear stress is zero, we 
find 


D, 


;IPI = 


K+l 


(2.45) 


AdJ)"' 1 2 

Equations (2.42), (2.43), (2.44), and (2.45) are valid for both materials. 
Therefore, the stresses and the displacement fields in terms of the Fourier 


Integral for each material may be expressed as 
Material (1) Q£x£h,0<y<°° 
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2j^<W*,>0 = |[ o °°[(C 1 a+C 2 {^L-!.+ ac})^ 

K i — 1 

+ (-C 3 a+C 4 { — ccx})e ^cosyada 

+^J i4i(-l + |p|y)e-^PIg«Pdp, 

1 2f°° k, + 3 

— o 1>7 (a:, >-) = -( q [(_c ia -C 2 {— +ac})e^ 

Kj+3 

+ (C 3 a - C 4 { — ax } ) e - *® ] cosy a da 

— oo 

- |f„“ « (-Ci«-C 2 (^ + 0«|) «» 

K. + l 

+ (-C 3 a + C 4 { — ocx }) e~ m ] sin ya da 

+-J- f iA^ye^lPle^Pdp, 


2 C°° 

u x {x,y) = [(C 1 + C 2 x)e xa + (C 3 + C li x)e~ xa ]cosyada 

»i C*o0 = U~tc-C’ 1 -C 2 { })«*“ 


K 1 

+ (C 3 -C 4 { — -x})e 

I f« Ai Ki + 1 

+ 2if-.ipi ( — ■ *W‘* K »*- 

Material (2) h <, x , 0 <, y < °° 

As x-»°o , u, v, o^, Gyy, are bounded. Then 


(a) 


(b) 


(2.46) 


(c) 


(d) 


(e) 
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2 ^ a 2 xx( j: >>’) = |f Q [-C 5 a+C 6 {^—- ax}]e *“ cosyada 

+ J_f 00 i 4 2 (-l + |P|y)^IPIe^^, (a) 

^-o 2yy (x,y) = [C 5 a-C 6 {-^— -ax}] e^cosyada 


-^o^CxoO = ^ [-C 5 a+C 6 {^ ax}] e^sinyada 


+-^J i'A 2 Pye yl P I c/p , 

2 r°° 

u 2 (x,y) = - [ C 5 + C^x ] e cosya da 


1 f°° ^ ( l|l^_ l ._p };)e -y|p| e <A : p rf p 

2rcJ^ 0 |p| P 2 K ^' 

w 2 (*o 0 = ; f [C 5 -C 6 l^-x}]e- xa sinyada 

7t J q (X 


l r~ a 2 k 2 +1 
+ 27tJ_oojp] 2 


+ |p|y)e-^PI^dp. 


(b) 


(2.47) 


(c) 


(d) 


(e) 


where the unknowns C • , 0' = 1 ,6) are functions of a, and Ay ,(/= 1 ,2) are 

functions of P, and are to be determined from the boundary and continuity 
conditions. 

The homogeneous boundary and continuity conditions of the problem 
described in Figure 2-5 are 
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CT 1 «(°>>') = 0 ; 0;Sy<<», (a) 

o iAy(°0') = 0 ; 0£y<<», (b) 

°La(*O0 = 02xx( h O') ; 0< y<~, (c) 

^i xy (h,y) = a 2xy (.h,y) ; 0£y<«, (d) (2.48) 

«i (h,y) = u 2 (h,y) ; 0£y<~, (e) 

v l (h,y) = v 2 (h,y) ; 0£y<~. (f) 

The mixed boundary conditions are 
» 1 (x,0) = 0 ; , b^xZh 

), (a) 

u 2 (x,0) = 0 ; h<>x<>a 2 , b 2 £x£oo 

(2.49) 

Olyv(^.0) = ’Oi W \ a^<x< b x 

}• (b) 


a 2xy (x,0) = -c 2 T (x) ; a 2 <x<b 2 

where, < 3 -^, a 2 T are the thermal stresses which can be determined from the 
uncracked problem. 

First, let us define the following density functions 

♦i(*)-^!^ ;0 <x<h, (a) 


w = ^ v 2 (x,0) ; h<x< oo . 


(2.50) 

(b) 


If we substitute equations (2.50) into equations (2.46)(e) and (2.47)(e) we 
obtain 


i i-oo /SAi k, + 1 

Hm — f -^-(-^r— +|p|y)«~>’ lpl ^dp=<t) 1 (j:), (a) 

y-* 0 2n}^ a |(i| 2 

(2.51) 

i r <« /BA? k? + 1 irt , . 

lim — I — (-^+|p|y)e~ ylpl ^dp = <S> 2 C*)^ (b) 

y 0 2?w | p | 2 


By applying Fourier inverse Transform, we have 


26 




(a) 



(b) 

(2.52) 

From condition (2.49)(a), it may be seen that 



$!(*) = () ;0 <x<cii , b\<x<h, 

<j> 2 (x) = 0 ; h<x< a 2 , b 2 < x<°°. 

(a) 

(b) 

(2.53) 

Then, A 1( A 2 can be put in the form 




(a) 



(b) 

(2.54) 

By applying the boundary conditions (2.48), 

and using Fourier Inverse 

Transforms and equations (2.54), with the help of the integration formulas given 

in Appendix (B), we find 



5 l“ e l^2“ fi 3 + e l^4 ” Kl + i J 

(a) 


B l + e 2 C 2 -B 3 + e 2 C A - K ^ + 1 J 

(b) 



e^B i + e 3C2 “ mB 3 + ^5^ + B5 + ^ 

= M 3'l , i (, i ) ‘ i 'i + ^TT J* 2 ' l/4 '^ < ' 2> * 2 • (c) 
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-e 4 5j + e 7 C 2 - mB 3 + e 8 C 4 +B S + e 9 C 6 

= ~ + 1 J A ^5 < t , l( r l)^ r l + K2 + 1 J M 6^2^2^2' 


(2.55) 


(d) 


e 10^1 + e ll^'2 + ^3 + e 12^4 _ ^5 _e 12^6 

- ^C Ml wm - + ^pn C M > )>A)dl2 ■ 


(e) 


-e 10 B l + e n C 2 +B 3 + e l4 C 4 -B 5 + e 15 C 6 

(f) 


where, = C^a, B 3 = C 3 a, 5 5 = C 5 a, m = an d e’s and. M’s are given in 

Appendix (C). 

Equations (2.55) form a system of six equations in eight unknowns C i , 

(/= 1 ,6) , §j,(j= 1 ,2). After lengthly calculations, we can solve the first six 

unknowns in terms of 6 , (/ = 1 ,2). The results are presented in Appendix (C). 

The two unknowns fy,(/= 1 ,2) can then be obtained by using the mixed 

boundary conditions (2.49)(b). Noting that 

a lyy (x,0) = -o 1 r (x) ; a l <x<b l . (2.56) 

By substituting from equation (2.46)(b) and taking the limits y—»0 + we 


have 

O fOO Kj +3 

lim 2p : — [( r C 1 a-C 2 {— — -+ax})e xa 

y—> 0 ttJo ^ 

+ 3 

+ (C 3 a - C 4 { — - — -ax}) e - * 0 -] cosya da 


+ 2^f A 1 (-l-|p|y)e-MPl^d(l=-a 1 7 ’W . 


(2-57) 


Or 
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2 ^J~[(^ia-C 2 {^+ca})e^ 

+(C 3 a-C 4 {-^ CLx})e~ xct ]da 

+ Hm 2^-Lj" A^-l -| (} |y) P I = -C l T (x) . (2.58) 


By using equation (2.54)(a) and taking the limit, the last integral of 
equation (2.58) will be equal to 
2 1 f*i $i( r i) 


2 1 9iVi 

JTK 1 + lJ ai 7 !-J 


dt , 


(2.59) 


Substituting into equation (2.58), and dividing by , we may 


have 


(Kj + ^f^-C^-Cjt^+OBC})^ 
Jo 


Kj + 3 

+ (C 3 a - C 4 { — occ }) e ~ xa ] rfa 

-df, =- 


r b i <t>i( f i) *(Kj + 1) 

‘J a 7^7*1 : 4 ^- <*i (*) •- 


(2.60) 


Similarly for the second condition (2.49)(b) we find 

(*2 + l)f (C 5 a-C 6 {-^- <xx})e~ xa da 

Jq i 

rb 2 ty 2 (t 2 ) 7t(K2 + l) 

+ - — -dt 2 = — a 2 '(x) ; a 2 < t 2 < b 2 , a 2 <x< b 2 . (2.61) 

J a 2 h~ x T*2 


Substituting for C,,(i = l ,6), from Appendix (C), we obtain the following 

singular integral equations for the unknowns fy,(/'= 1,2): 
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(2.62) 


f fr — +^n(x 1 ^j)]<l>i(/ 1 )Ai+ f k^ 2 {x^,t 2 )^ 2 (t 2 )dt 2 
J a l h~ x l J a 2 

JC( K l+i) T 

■ — : 4 ^—' °i C*i) ;«!<*!< 

r b 1 tb 2 1 

^2i(^2 ,? i) < t > i(^i)^i + I f ; — — + k 2 2(x 2 ,t 2 ) ^ ^2^2) ^2 

J flj J a 2 r 2“*2 

^ + 3) 

=- — <S 2 l {x 2 ) ; a 2 <x 2 <b 2 . 

where 


*n(*i> ? i) = f G n (x v t v a)da , 
J o 

(a) 

k 12 (x v t 2 ) = f G l2 (x v t 2> a)da , 
J o 

(b) 

^2i(^2’ r i) = f G 2i {x 2 ,t l> a)da , 

J o 

(c) 

k 22 (x 2 ,t 2 ) = f G 22 (x 2 ,t 2 ,a)da . 
Jo 

(d) 


^(Vj.a) = [-/fj-C^+a*!)!^ e*i a 


+ [fl’ 5 -(^-ax 1 )i/f 7 ]e^i“ , (a) 

G n (x v t 2 ,a) = (^1) {[-// 2 -(^Ii + ox 1 )i// 4 ]^« 

+ [H 6 -^-ouc 1 )^H s ]e- x i a ) , (b) 

KU + l Ko + 3 

G 21 (x 2 ,? 1 ,a) = (— ){[// 9 -(— , (c) 

Kj+3 

G^ix^t^O.) = [^10-( — 2 — — cuc 2^12l e ^ 2a * (d) 


where ,(j =1, 2,.. ..11, 12), D, can be found in Appendix (C). 


(2.63) 


(2.64) 


(2.65) 
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For internal cracks, i.e, a l > 0 , b l < h , a 2 > h , b 2 < 00 it can be shown 
that the kernels k - , (i J = 1 , 2) , are bounded for any combination of the variables 
Xj, tj, a, (j =1, 2). So, the two singular integral equations (2.62) and (2.63) are of 
the Cauchy- type. 
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2.2 Model II - Layered Medium on an Elastic Foundation 

2.2.1 Temperature distribution 

2.2. 1.1 Unit step function at the boundary 



Figure 2*6: a-Geometry of the problem (Model II) b-Temperature boundary 
condition (unit step fiinction) 

Consider two dissimilar infinite layers having thicknesses h x (coating) and 
h £ (base material) as shown in Figure 2-6(a). Figure 2-6(b) shows the 
temperature change at the boundary. Let D x , k’^ , D 2 , k' 2 be the thermal 
diffusivity and thermal conductivity for the coating and base, respectively. 


Then, the differential equations for 1, 2 are 


cP-e x (x\t) _ i dQ x (x\t) 
dx* d[ aT~ 


-h x <x'< 0 , t > 0 , 


(a) 
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(2.67) 


where 


©!(*’, 0) = 0 2 (X’,O) = 0, 

(a) 

0 1 (O,r) = 0 2 (O,f), 

(b) 

w 9©i(o,f) m a © 2 ( 0 /) 

(c) 

*1 3 k 2 ^ 

dr dr 

dx* °’ 

(d) 

e l (-h v t) = e 0 H(t). 

(e) 

® l (x\t) = T l (x\t)-T co , 

(a) 

® 2 (x’,t) = T 2 (x’,t ) - T „ , 

(b) 

I 

II 

of 

(c) 


( 2 . 68 ) 


T x is the initial temperature for both materials, and T 0 is the temperature 
at the boundary at any time, t> 0. Condition (2.67)(c) implies no resistance at 
the interface, and condition (2.67)(d) implies insulated outer boundary. 

Again, if we apply Laplace Transform to equations (2.66), and use 
condition (2.67)(a), we will end up with two ordinary differential equations that 
have the general solutions in the form 


©jCx’p) = A l (p)coshq l X’+B l (p)sinhq 1 X’, (a) 

(2.69) 

©2 ( x ’iP) — A 2 (p) cosh q-^c + B 2 (p) sinh q^X' . (b) 

where Aj(p), B-^p), A 2 (p), B 2 (p) are four constants that can be determined from 
the conditions (2.67)(b-e). So, after some manipulations equations (2.69) will 


become 


®i(*’P) = 


® 0 [cosh q 2 h 2 coshq^X'-T] sinh q 2 h 2 sinh q^X'} 
p[cosh q^hi cosh q 2 h 2 +r\sinh q l h l sinh ’ 


(a) 


0 2 (x’p) = 


® 0 [cosh q^ cosh q^X’ - sinh q 2 h 2 sinh q^c ' ] 
p[coshq^h ^ coshq 2 h 2 +t\ sinhq^h^ sinhq 2 h 2 ] ' 


(2.70) 

(b) 
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where 


«?! 


I p I p ^2 \ D \ 

-**£■ - <72=^ • n — ^7T- 


(2.71) 




Appling the Inversion Theorem the temperatures Qy(r,t), 0 2 (x , ,f), are 
found to be 


1 rY+«~©o 

e ‘ ( " )= 2*L- 


1 j* Y+i»0 o [cosh^2^2 cosh^x' -ri sinh^^ sinh^yX’]e a 

y_ tOB <. [cosh^yhy coshLfJfi^ + T| sinh^yhy sinh^h^ 


dz, (a) 


1 r V+ ia 

e ^’ ,)= ssL 


Q 0 cosh^ 2 (x’-h 2 ) e tz 


(2.72) 


y_ ioo z[cosh l^yhy cosh ( s2 h 2 +1) sink ^yhy sink 


dz. (b ) 


where 


! h^z/Dy , ; 2 =<z/D 2 . 


(2.73) 


It is clear that the integrands in equations (2.72) are single-valued functions of 

z, with simple poles at z = 0, and at the roots of the equation 

cosht^yhy cosh ^ 2 ^ 2 + 1 ) sinh ^yhy sink C, 2 h 2 = 0 . (2.74) 


Let, 

= 'Iz/Dy = 1(0,2 = -Dy(S?. 


(2.75) 


Then equation (2.74) becomes 

cos cahy cos 5co/! 2 -"n sin (ohy sin S(oh 2 = 0 . 


(2.76) 


where 



(2.77) 


It can be seen that the roots of equation (2.76) are real and symmetrically 
placed with respect to co = 0. The roots co >■ are obtained from the intersection of 
the two families of curves y = cot co hy, y = r|£an Sco/i 2 . The simple poles of 
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equation (2.76) are distributed along the real axis in the negative direction, and 
are equal to 

z m = ~ D l a i > mm 1,2,3 (2.78) 

2 

The residue at these poles (z = 0 , z = -DjO^) can be found in Appendix 
(D). Then, by using the Residue Theorem, the contour r shown in Figure 2-7, 
and letting 

*2 


x < n* 

- = (X*-1) , RmjZ , X m = A l (O m , T = — , 
"1 A 1 *1 


(2.79) 


the closed form solutions for ©j, 0 2 are found to be 



Figure 2-7: The contour for evaluation of integrals in equations (2.72) 


0i( ** ,T) = 1-2 £ 


©n 


m=l 


[cos X m (x* — 1) cor A, m 8/? +T| sin X w (x* — 1) «'n A. m 8/?] 

[(1 + 1|5/?) sin A. m cos A. m 8/? + (5/? +Tj)co.s k m sin A, m 5/?] 

(a) 
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(2.80) 




Figure 2-8: a-Geometry of the problem (Model II) b-Temperature boundary 
condition (ramp function) 

If we apply a ramp function at the boundary. Figure 2-8(b), the 
differential equations for the materials 1, 2 , and the initial and the boundary 
conditions are the same as equations (2.66), (2.67), except that, condition 
(2.67)(e) should be changed to 

©1 (rhjt) = ^ (t-t 0 ) H{t-t 0 ) . (2.82) 

r 0 f 0 

Again, by applying Laplace Transform, we obtain equations 
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Q l (xp) = F 21 (p) ,B 21 (p), 


(a) 


©2 ( X <P) ~ F 22 (p) • B^ip) • 


where 


F 2 i(p) = 


0 o [cosA <? 2 ^2 cos ^ <7i*’ - q <? 2 ^2 Sin ^ <7i*’] 

/>[c0jA? 1 A 1 coshq 2 h 2 +r\ sinhq^h^ sinhq 2 h 2 ] ' 


(b) 


(a) 


(2.83) 


F 22 (p) = 


9 0 [cosh q 2 h 2 cosh q^:.' - sinh q 2 h 2 sinh q^c<] 
p[cosh q^h j cos/t < 72 ^ 2 + T 1 sinh q fa sinh q 2 h 2 \ ' 


(b) (2.84) 


1 1 e~P‘o 

B 2l {p) = B 22 <lp) = -{— e —). 

*0 P P 


(c) 


By applying the Convolution Theorem to equations (2.83), and observing 
that F 21 , F 2 2 are the same as equations (2.70), and letting 


X’ , * „ ^2 , , tD \ t<Pi 


- = (x*-\) ,R = - , X m = h 1 co m , X = — , T 0 = 
n i n \ hf hf 


(2.85) 


the temperature distribution for strips 1, 2 become 

—s f +2 I <^"-D 

U 0 T 0 m=l 

[c<w X m (x* - 1 ) cos X m 5R + q sin X m (x* - 1 ) L8/f] 


t 0 ^ m 3 [(l + q8/?) sin X m cos X m bR + (6/? + q) cos X m sin A, m 8fi] 

; 0<**<1 , t£t 0 , (a) 

2 , „ i 2 


0 i(**,x) ~ , , , , 

= 1 -2 Y eA ( e W- 1 ) 

0 O m=l 

[coi A, m (x* - 1 ) cos X m 8/?+q s/n X m (r* - 1 ) sin A. m 8R] 
Xo^ m 3 [(l +q5/?) sin X m cos X m bR + (8 R + q) cos X m sin X m 8Z?] 


; 0£r*£l , t>t 0 , (b) 
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( 2 . 86 ) 


-V - = f +2 £ (e ' -- I> 

U 0 X 0 m=l 


cojA 8(x*-l-/?) 


T 0^m J [(l +t ) 5^) -fm A m cos A^S/? + (8R + T)) cos A^ sin A m 8/?] 

; (/?+l) , t£t 0 , (c) 


®2^ > x ) xX 2 / xX 2 

— p: — = 1-2 X e tX « (e^m -1) 

W 0 m=l 

cosA m 8(jc*-l-/?) 


x 0^m 3 [(l +t | 8#) sin A m cosX m &R + (8/? +T|) cosA m sm A m 8/f] 

; 1 ^ x* ^ (/? + 1) x > t 0 . (d) 

where 8 = ^D^/D 2 , and A m are the roots of equation (2.81) 


2.2.2 Thermal stresses in uncracked problem 

Consider a thin- walled hollow cylinder with cladding of thickness Aj, and 
base metal of thickness h 2 , shown in figure 2-9(a). 

Let us take an element from the composite shell of unit length and a very 
small angle d§, as shown in Figure 2-9(b). Assume a uniform strain over the 
shell thickness e^, e 0( , then from the Hooke’s Law we have 

% = ^ 1 (a ly- V l a lt ) + a l 0 l = ^ (a 2y-V2°2t) + a 2 e 2’ («0 

(2.87) 

e 0 1 = + a l 0 1 = £^°2/- V 2°2P + a 2 0 2 • (b ) 

where <7^, a^, <s 2yt are the stresses in axial and tangential directions for 1, 
2, respectively, v 1( J Bj, aj, v 2 , E 2 , are the Poisson’s ratio, Young’s modulus 
and thermal expansion coefficients, respectively. By integrating equations (2.87) 
over the shell thickness we obtain 
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( fl ) 


Figure 2-9: Geometry of thin- wall hollow cylinder with cladding to calculate 
thermal stresses in uncracked problem 

e 0y = ^(o ly -Viai,)+a’0i av =i-(o 2y -v 2 a 2r )+o^0^ v , 

(a) 


(2.88) 

“or = ^ (° lr- v l“ l y ) + “i ©lav=^ (“2r" V 2 “2>) + «2 0 2av 

(b) 


where 
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(2.89) 


1 f*i . 1 f*i . 

* ° lf = ^J 0 ° 1( 


- i r*i+ A 2 

o 


l r n i +ft 2 — l r"i +n 2 , 

v* • °a = /d,. ^ 


1 f fc l 1 f +/I 2 

®iav = ttj ®i(xs)dx , ©2av = r-f ® 2 (x t i)dx 
n l J 0 h 2 J h x 


1 ( k \ +h 2 
1 f*l+*2 


The resultant force in radial and axial directions must vanish giving 

• h, ph, +A, 

(a) 


r n i r n i +n 2 

j o it dx+\ a 2t dx = 0, 
J 0 J A, 


t h 1 f *1 + *2 

J„ °iy‘ i[+ j Ai 


(2.90) 


(b) 


Equations (2.90) can be put in the form 
O + ©2^2 = ® > 


CTi/i + o 2 / 2 = 0. 


(a) 


(b) 


(2.91) 


Equations (2.88) and (2.91) are six equations with six unknowns (e,^, e^, 

a ly , o lt , o 2y) Oa,). By examining these equations we can see that they would be 
unaltered if the subscripts y and t were changed. Therefore, 

(a) 


0 ly -0 lf _0 l» 

°2y = °2f = °2> 
e 0y = e 0r = e 0- 

Then, equations (2.88) and (2.91) are reduced to 


(b) (2.92) 

(c) 
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(a) 


e o = i- ( ai- v iCi ) + a i 0 i av , 

e ° = 4 ( ^ 2 “ V2?2)+a 2 02m ” (b) (2 ‘ 93) 


<J l h 1 +o 2 h 2 = 0. 

By solving equations (2.93) together we can get 


(c) 


E l ®1 ®lav E 2°2 Q 2av 

— h,+— hy 

1-Vj 1 l-v 2 ^ 

20 El E 2 

i Vi h 2 

1-Vj 1 l-v 2 2 


(2.94) 

Also, from equations (2.87) we can see that 



°ly = ^lf = 0 l > ^2y = = °2 • 


(2.95) 

Then equations (2.87) become 



1-Vi 

e o = -=^o 1 (x,f)+a , 1 Qi(x,t) , 

(a) 


l-v 2 

e 0 = — — a 2 (x,t) + e 2 (x,r ) . 

(b) 

(2.96) 


’2 


Substituting equation (2.94) into equations (2.96), the thermal stresses in 
materials 1, 2 can be written as 


<*i(*4) = yz^teo- a i ©i(*»0] , 


(a) 


E 2 

<* 2 (x,t) = — — [fio-o^ © 2 (x,r)] . 
i-v 2 


(2.97) 

(b) 


Thus, in the concentric cylindrical shells having a radial temperature 
variation, the axial and tangential stresses are given by equations (2.97) at all 
points except at the ends of the cylinder. The end surface tractions produce 
thermal end moments at the rim of the cylinder, and these generate axial 
bending stresses which disappear rapidly with the increase of distance from the 
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end. Let 


h 2 


tD, 


x = t ' R = r • x = tt 

«i «i 

Then, equation (2.94) becomes 


(2.98) 


t l ©! (x*,x) E 2 1 - vj f i+£ & 2 (x*,x) , 

— ^ dx +- 1 ~ dx 

Jq ©o E ^ 1 V2J1 ©o 


Eq - a’j ©o [ 


1 + 


^2 1 ~ v i 
E \ 1 _v 2 


(2.99) 


R 


£ Oj' Q 

By letting o 0 r = — y^T - " > from (2.97) and (2.99) the thermal stresses 

materials 1, 2 are found to be 
<Jl(**/t) _ ©i(**/t) 

a 0 T % 

(. 1 ©j^*,!) ^ # E 2 Cdj 1 - Vj j* 1 +R & 2 (,x\x) 


m 


dx* 


a; 


1 + 


e 2 1 ~ v i 
E ~i^i‘ 


;0£a:*£ 1, (a) 


( 2 . 100 ) 


o 2 (x*,x) E 2 l-v { ^ Q 2 (x*j) 


E\ 1 — Vo 


a’, 


©r 


t 1 ©i^.T) E 2 (% 1 -V t r 1 +R Q 2 (x*,x) 
— ^ dx +- — TT dx 

Jq ©0 E j t 1 V2J1 ©0 


1 + 


a 

JL 

E 2 1 ~ V 1 
£1 1 -v 2 


£ 


;l<x*=££ + l. (b) 


where the integrals 
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( 2 . 101 ) 


|»10 1 (.X*,T) 

Jo ©0 


dx* 


• l +rQ 2 (x*,t) 

— ^ — d* 
1 ©o 


are given in Appendix (E). Equation (2.100) give, essentially, the thermal 
stresses in a composite plate constrained .to remain in its own plane. 

2.2.2. 1 Stiffness of the elastic foundation 



Figure 2-10: Geometry of thin-wall hollow cylinder with cladding to calculate 
stiffness of elastic foundation % 

It is well known that the thin-wall hollow cylinder can be simulated as a 
beam on an elastic foundation of stiffness related to the radius of the neutral 
surface of the cylinder, the thickness of the cylinder , and the Young’s modulus, 
[5]. 

In case of thin- walled hollow cylinder with cladding, Figure 2-10, the 
stiffness of the elastic foundation % can be determined by taking a strip AB of 
unit width cut from the composite cylinder. Assume that x n denote the change of 
the radius R n (radius of the neutral surface, e = 0), i.e, x n = A R n , which is the 
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deflection of the strip at any section measured from neutral surface. The strain 


over the cross section e f , would then be 
x„ 

n n 

Now, the corresponding stresses are 

E \ X n 




E 2 x rt 


° 2 '~ * 


Then, the force per unit length on the strip AB is 
E,x„ EjX„ 

F ’-ir^ir^ 


The radial component of this force is 


E \ x n 


E ?n, 


F t d$ = [^h 


E l h 1 +E 2 h 2 


R. 


] x n’ 


( 2 . 102 ) 


(a) 

(2.103) 

(b) 


(2.104) 


(2.105) 


where d§ = 17i? n . 

That is, radial force is proportional to the deflection x n . So, we conclude 
that a longitudinal element of a cylindrical tube loaded symmetrically with 
respect to its axis can be regarded as a beam on an elastic foundation, having 
the stiffness 


E l h l +E 2 h 2 



(2.106) 


where R n is the radius of neutral surface of the cylinder. To obtain R n , let us 
consider the composite beam shown in Figure 2-11. 

Let p* is the radius of curvature of deflected beam and x n is the distance 
from the neutral axis, then 

e=- (2.107) 

P* 


So, the stresses in 1, 2 become 
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Figure 2-11: Geometry of composite beam to calculate the radius of neutral 
axis R„ 


c iy = 


o 2y = ' 


E &n 

I 

P ’ 

E 1?n 

P* 


The total force in y-direction must vanish, then 

J °l y ^l + | c 2y dA 2 = 0 • 


(a) 


(b) 


(2.108) 


(2.109) 


where Aj, A 2 are the cross section areas of 1, 2, respectively. Substituting 

equations (2.108) into equation (2.109) we may have 
E E 

— fx n dA l +—(x n dA 2 = 0. (2.110) 

O’ J O’ J 


Since 


j X n e ^ A 1 X 


(a) 


( 2 . 111 ) 


J X - 


^2 “ x 2n A 2‘ 


(b) 


where x ln , x 2n are the distance from the neutral axis to the center of the cross 
sectional area for strips 1 and 2 respectively. Then equation (2.110) can be 
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written as 


~ E \ X ln A l +E 2 x 2n^2 ~ ® • ( 2 . 1 12 ) 

Also 

Xln + X2n = \( h l + h 2)- ( 2 ‘ 113 ) 


Solving equations (2.112) and (2.113) for x ln , x 2n> we have 
— 1 .. . . E 2*2 


(2.114) 


Since, the beam has a unit width, i.e, Aj = h^l, A 2 = ^l* ^ ien equation (2.114) 
becomes 


— 1 ^2^2 

So, the natural axis radius R n is 
— 

R n = R i + Y +x 1"’ 


(2.115) 


(2.116) 


where R^ is the inner radius of the cylinder. 


2.2.3 Formulation of the crack problem 

Consider the problem that is shown in Figure 2-12, where % is the 
stiffness of the elastic foundation. The analysis is the same as in Model I, and 
equations (2.41), (2.42), (2.43), (2.44), and (2.45), are still valid. Then the 
displacement fields and the stresses are 
Material (1) ,0£jt£A 1 ,0^y<«> 

| /"[(qa+Cj^+o.))^ 

K.-l 

+(-C 3 a+C 4 { — - ax})e ^cosyada 


+J-f~ A^-l + lplylc-ylPl^dp, 

ZTZJ _ 00 


(a) 
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Figure 2-12: Crack geometry Model II 


- sJ 0 " t( ' c ‘ a ' CjI 'T‘ +a,|,c “ 

Ki + 3 

+ (C 3 a - C 4 { — etc}) e~ xa ] cosya da 

- 2| o "[(-C 1 a-C 2 (^ + ax|)e» 

Kj + 1 

+ (-C 3 a + C 4 { — - — - ax } ) e~ xa ] sinya da 
+2-f~ iA^ye-y^e^d^, 

271 

U]ix,y) = - f [(C 1 + C^x) e** 1 + (C 3 + C 4 x) e _;ta ] cosya da 

KJq 


1 f°° |P| K 1 -1 


;/ 


C- 


27tJ_«|pr p 2 


— Py)e MPlg^Prfp, 


(b) 


( 2 . 117 ) 


(c) 


(d) 
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i>ic*oo =-r })«** 

7C Jq (X 

+ (C,-CA—-x})e~ xct ] sinyada 
a 

+2^J~ <«) 

Material ( 2 ) , &x £ , 0 £;y < °° 

2^o ltt (xo') as |J o [(C 7 a+C 8 {^-+ca})e* a 

K 2 -I 

+(-C 5 a+C 6 {— - ca})e -Jta ] cosyada 

+^-f A 2 (-l + |Ply)«"^PIc“PdP, (a) 

^+3 

+ (C 5 a-C 6 {— ^ occ})e _;ta ] cosyada 

+l-j A 2 (-l-\$\y)e~ y W\e ix Pd$, (b) 

- i J o "[(-C 7 a-C 8 (^l + c«|)^ 

Kj+ 1 

+ (— C 5 a + C 6 { — - — — ax } ) e~ m ] sinya da 

+^-J (c) 

« 2 (x,y) = - f [(C 7 + C 8 x) + (C 5 + C^) e -Jta ] cosyada 

nJ 0 


( 2 . 118 ) 
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v 2 M = l r[(-c 7 -c 8 {-+*})**° 

kJq a 

Ko 

+ (C 5 ~C 6 {— -x})e ^]sm;ya<2a 

(e) 

where Cj , (j = 1, 2, 8) are functions of a, and A 1; A 2 are functions of p, and 

are to be determined from the boundary and continuity conditions. Eight of 
these conditions are given by (2.48), (2.49) while the two remaining conditions at 
x = h-^ + h 2 are 

°2xx( h l +h 2 >y) = -Px = - X«2 ( /! l + h 2 O0 • (g ) 

C 2xy (h 1 + h 2 ,y) = 0. (h) (2.119) 

By using the same definition of the density functions <(>^(jc), (j^Oc), as in 
equations (2.50), and by using the boundary conditions (2.48), (2.119) and the 
Fourier Inverse Transform, with the help of the integration formulas given in 
Appendix (B) we will have 

~ B l~ e l C 2~ B 3 + e l C 4~^~T' 7 J 
^l + ^2 ( -'2~ B 3 + e 2^A = 7r~rT f 

K l + 1J a 1 

e 4 B 1 + e 3^2 “ mB 3 + e 5^4 + 85 + e 6^'6~ e l0 B l + e 25^S 

■Jr !-;- £' ■ 

~ e 4 B l + e 'j^2~ mB 3 + + &5 + e g^6 + e 10^7 + e 16^8 


(a) 

(b) 

(c) 

(d) 
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(2.120) 


e 10^1 + e 11 ^ 2 + ^3 + e 12 ^ 4 “'® 5 - e 12 ^' 6 - f 10 ^ 7 - e 11^8 

= 5^7iC M ^) d h (•) 

~ e 10® 1 + e 13^ 2 + ^3 + e 14^4 + e 15^ 6 + e 10^7 + e 17^ 8 

/*' *W»1>*1 + ^jVo'W<2>‘*2 • <0 

l r*2 

- ^5 + g 2oQ + e 18^7 + e 19^8 = v . .1 ^1 1 ^ 2 ^ 2 ) ^2 ’ 0 ?) 

K 2 +1 J a 2 


(p - 1 )B 5 + (p«24 + ^22^6 + (”P ~ l) e l8^7 + (P e 23 + e 21^8 

= K2+1 J a M l2^ t 2> dt 2' 


(h) 


where Bj = a , B 3 = a C 3 , B 5 = aC 5 , B 7 = aC 7 , e’s and M’s are shown in 
Appendix (F), and p = %/2\i 2 ol. 

Equations (2.120) give eight equations in ten unknowns. After long 

manipulations, we can obtain the unknowns Cj , (j = 1, 2 8) in terms of the 

two other unknowns (j^ , The results are shown in Appendix (F). 

By using the mixed boundary conditions as in Model I, we obtain two 
singular integral equations of Cauchy-type in the same form as equations (2.62), 
(2.63), where G,y , (i ,j= 1 , 2 , ) , in this case are 

G^^a) = [-e 1 -(^ + cu 1 )G 3 ]^i“ 

K +3 

+ [Q 5 -(,- 2 (a) 
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G 12 (x v t 2 ,a) = 

+ [Qe~ 

G 2l {x 2 ,t l ,a) = 
+ [Qg~ 

G 2 2( x 2^2’ a ^ = 

+ C2 10 

in which Qj , (/• : 
a, < x t < b { , a. < f- < b i , 


Kj + 1 


Kj + 3 


*2 


TT {[-fi 2 -(-V' + ^i)24]^ a 


k 2 +3 


•(— 5 CKjXJgl^ 0 } . 


x 2 +1 . 2 b k 7 + 3 2i5 

{ [“— -(-T“ + 0«2)3-] 


*1 + 1 *>o 


k 2 + 3 

■<r=2 ox^Q n ]^) , 

Cm *2 +3 2 16 

Kj + 3 

-(-i- ca 2 )Q 12 ]e~ x 2 a . 


(b) 


(c) 


(d) 


( 2 . 121 ) 


= 1, 2, , 16), D 0 , are found in Appendix (F), 

0 = 1 , 2 ). 
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Chapter 3 

THE SINGULARITY AT THE CRACK TIP 

It is well known that the stress fields around a crack tip is proportional to 
r~ s , where r is a small distance from the crack tip at which we measure the 
stress field, and s is called the power of singularity which should be between 
zero and one, i.e, 0 < s < 1 . If s is less than zero, the stress is bounded as r— > 0 , 
and there is no singularity at the crack tip. If s is greater than one, the strain 
energy density is unbounded and goes to infinity as r — > 0, which also is 
impossible. 

The value of singularity s is dependent on the crack configuration as well 
as the material properties. For each special crack configuration, some terms of 
the kernels in the singular integral equation (2.62), (2.63) become singular. By 
using Muskhelishvili’s technique [37], the singular behavior of the stress state 
at the crack tip can then be examined for each case. 

Since, the singularity at the crack tip will be the same for both Models I 
and II, the results would be valid for both models. 

3.1 Embedded crack in both materials 

The case of a crack embedded in both materials is shown in Figure 3-1. 
The only singular terms in the integral equations (2.62), (2.63) are the dominant 
terms — - — and — - — , other kernels are bounded as a— >«. The two singular 

'i-*i *2"^ 

integral equations can thus be written in the form 
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Figure 3-1: Geometry of embedded crack in both materials 


p^i^i) 

Jojirti 


dt { + B.T. 


rcO^ + 1) _ 

— c 1 T (x l )-,a l <x l <b v 


(a) 


t b 2^ t 'i) 


dt 2 + B.T. 


7C ( k 2 + 1 ) t , . , 

- ■ 4|i ~ - o 2 \x 2 )\a 2 <x 2 <b 2 . 


(3.1) 

(b) 


where B.T. correspond to the bounded terms. 

To examine the behavior of the unknown functions around the 

irregular points (end points), following Muskhelishvili [37], we assume that the 
unknown functions <j)j , <j> 2 may be expressed as 






{t r ap*j{b r tjpj 


- = gj(tj) Wjitj) ; 7 - 1 . 2 . 


(3.2) 


where gj(tp satisfies a Holder condition in the closed interval aj<. tj<, bj , (/= 1,2) 
and g(aj) *■ 0 , gfbp * 0 , (/= 1 ,2). Also ay , P • , (/ = 1 , 2) are the singularity at the 
irregular points which should satisfy the condition 0 < Re{ ay , py) < 1 , (/ = 1,2), and 
Wy(fy) is any definite branch which varies continuously on the interval 
a j< t j < bj .(/=!. 2). 
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Define the following sectionally holomorphic function 


(3.3) 


a i J 


Substituting equation (3.2) into equation (3.3) we obtain 

1 f bj gjitpe^j 


,. (Z) = IP W g . 2 *. • i. 

7 ajttj-ajpjdj-bjpjQj-z) J ’ 


= 1 , 2 . 


(3.4) 


Following Muskhelishvili, equation (3.4) can be written as 

F(r) . 8j (ape*j gj (bj) 

7 (bj - ajftj (z - ap a j sinrnXj ( bj-ap a i(z-bph sinnfy 

+F 0 .(z) ; 7 = 1,2. (3.5) 

where F 0 y (z) is bounded everywhere except possibly at the end points 

aj ,bj,(J= 1 , 2), where it has the following behavior 

C jk 

|F 0 ;(2)I < ; J -~ ; j,k- 1,2. (3.6) 


where c ;1 =a ; , ej 2 =bj , p ;1 <Re(a.p ,pj 2 <Re($p , and Cj k ,pj k are real constants , 
i.e., .pQy (z) has singularity less than ay, Py, 0 = 1,2). 

Using the Pelemelj formula [37] 


1 fW'f) 1 + 


(3.7) 


which, by using equation (3.5), may be expressed as 
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1 f bj * i(f ?dt - g J (a P COtna j _ Sj(bj)cotK$j 

n hj tj-Xj 1 (bj - apPj (Xj - aj) a j (bj-ap a j(bj-xpPj 

+F 0j ( Xj ) ; 7 = 1,2. (3.8) 


Substituting equation (3.8) into equation (3.1) we can find 
gjtyco^ gj(bpcoti$j 

(bj - ajfij (xj - a j) a j (b r aj) a j(b r xjfj J J ’ 


(3.9) 


where \\fj (xj) , (j = 1 , 2) contain all the bounded functions. 

By multiplying equation (3.9) first by (Xj-aj) a j , and letting Xj-t dj , and 
then by (bj-xjftj , and letting Xj->bj,(j= 1,2) , we obtain the following 
characteristic equations for ctj , Py , (/= 1 ,2) 

gjiflp COtKCLj 


(b r aph 


= 0 ;y = 1 ,2, 




(a) 


(b) 


(3.10) 


Or 


cotnctj = 0 ; j= 1,2, 

(a) 

(3.11) 

cotn$j = 0 ; y-1,2. 

(b) 

The acceptable roots of this equation are = 5 . 1 

= ^ , 0 = 1 . 2). which are 


the known results in the crack problems. Hence, the fundamental functions of 
the singular integral equations are 
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1 


j= 1,2. 


(3.12) 


Wj(Xj) = - 


(Xj- a j) in (bj- x j) 


1/2 


Therefore as long as we have internal cracks , the power of singularity 
will be j . 

3.2 Edge crack 

3.2.1 Edge crack for Model I 

This case is shown in Figure 3-2, in which a 1 = 0 , < h . If aj goes to 


s f 



Figure 3-2: Geometry of edge crack (Model I) 

zero, the singular terms in the singular integral equation (2.62) are the 

dominant term — - — and some terms in the kernel *11 (*l ,fj,) • Let 

h x \ 

*n(*i-'i) = *n«(-Vi> + *ii (*iA> • ( 3 - 13 ) 

where k s ^ e is the singular term in edge crack case, and k^ is the bounded term. 
The singular part k s ^ e of the kernel is given by 
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(3.14) 


S f °° oo 

k Ue C*1 A) = G 7le (Vl*°) d0t 

Jo 


where G~ le is the asymptotic value of G n for a-> » , ^ 0 , x x -*■ 0 

(2.65) and Appendix (C) it follows that 

s f°° 1 K i Ki+3 

*li,(Vi) = J Q ^2 + Y ( " 1 +2fja)-(-L — cuj)(- 1 +2^a)] 

e~^h +x \) a da . 


or 


.s 1 6*i 

*iu — 


4XJ 2 


'l + *l aj+^i) 2 (*!+*!> 


,3 * 


Equation (3.16) can also be written as 

r * 1 U ^ d „ _ 1 

*11- - C- 2 *! 2 — - -fcrjr-nTTT 




(*i+*i) 


Then, singular integral equation (2.62) becomes 

C b i 1 j 

[ + ^u«( Jc i> r i)l <t>i(^i) dt\ + bounded terms 

J a x h~ x \ 

7t(K+l) T . 

= -~4^- a m) •*! <*!<*!• 


Again, assume that 








where (fj) , w x (t x ) have the same properties as before. 
Define sectionally holomorphic function 


Then, from 


(3.15) 


(3.16) 


(3.17) 


(3.18) 


(3.19) 
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( 3 . 20 ) 


^l(z) = - dt x 

7tJ 0 t x -z 


Substituting equation (3.19) into equation (3.20) and following [37] , we 
will have 




*i(0)e 


,jwa, 




(b^i (z) a i sin na x (b x ) a i(z-b x )$i sinn$ x 


+F 0l (z) . 


( 3 . 21 ) 


Then, as before by following [37], equation (3.21) becomes 
g l (0)cotvxt l g 1 (b l )cotn$ l 


Wi) = 


(b^i (x^i (b x ) a y (b x -x^i 


+/ ? 0iC*i) • 


( 3 . 22 ) 


and 


F i(“*i) = 


Si(0) 


(bj)Pi (jc^* 2 ! smitaj 




( 3 . 23 ) 


where F 01 (xi) , F^Otj) are similar to the Fq^x-) in equation (3.8). Substituting 

equation (3.22) and equation (3.23) into equation (3.18), we obtain 
g x {Qi)cotna i g l {b l )cotn^ l 

(bjAOcj) 0 ! ”(b 1 )“i(b 1 -*i) fi i 


r _ 2 d 2 , d 
+ [-2x x — -6*,— 
1 - l dx. 


- 1 ] 


<£c‘ 


1 


gl(0) 

(6j)Pi (*^“1 n'/inaj 




( 3 . 24 ) 


where ^J/ 1 0r' 1 ) contains all the bounded functions. Equation (3.24) can be written 
as 
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g^cotn^ g l (b l )cotn^ l 
(^Pl (*!>“! '(^“1(^1-^ 


+ 


8 1 ( 0 ) 

(ft^i (j^) 01 ! sinn^ 


[-2a 1 (a 1 + l)+6a 1 -l]=nf 1 (x 1 ). 


(3.25) 


Multiply equation (3.25) by (bj-Xj^i , and letting x 1 -tb 1 , the characteristic 
equation for Pj is found to be 

cotnfii = 0 -4 p 1 = I. (3.26) 


Similarly, multiplying equation (3.25) by (jCj) a i , and letting Xj — > 0 , the 


characteristic equation for 0 ^ is obtained as 


g 1 (0)cof7ta 1 gj(0) 


(bj)^ sirrm.^ 


[-2a 1 (a 1 + l)+6a 1 -l] =0. 


(3.27) 


which can be written as 

cus7ta 1 -2(a 1 -l) 2 +l =0. 


(3.28) 


Since 0 ^ = 1 is unacceptable, the only possible root of equation (3.28) is 
ocj =0 , giving the fundamental function as 


Wl(*l) = 


1 

Q> 


(3.29) 


3.2.2 Edge crack for Model II 

This case is shown in Figure 3-3, in which a 2 > , b 2 =h l +h 2 . 

If b 2 goes to h^h 2 , the singular terms in the singular integral equation 

(2.63) are the dominant term — - — and some terms coming from the kernel 

*2“ -*2 

* 22 (- x 2 * r 2> • Let 
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y 


E l> v 1 

E 2 , t/ 2 

a 2 6 2 


mm 

^ L 

L r - 


Li-^ — i 




4 




Figure 3-3: Geometry of edge crack (Model II) 


= hie (V2> + k 22 • 


(3.30) 


where k^ le {x^ ,t 2 ) is the singular term, (xj,^) is bounded and 

%2e = [° G~ le {x 2 j 2 , a) da . (3.31) 

J o 

G 22e is the asymptotic value of G 22 for a — > 

00 » *2 — ^ h l +h 2 ,x 2 ->h l +h 2 . From 

(2.121) and (3.31) it may be shown that 

*5a.(V2 >“ \~-[\+{l + 2(t 2 -L)a}{ 3 - + (x 2 -L)a}] * 

e -(-t 2 - x 2 + 2L)ad a . ( 3 . 32 ) 


where L=h^+h 2 . Evaluating the integrals equation (3.32) becomes 
'n.W* QL- H - h ) ' 


(3.33) 


which is similar to the equation (3.16). Thus, by doing the same analysis as in 
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section [3.2.1], we obtain the same singularities, i.e., 


^ , P 2 = 0 . The 


fundamental function would then be 


w 2 (x 2 ) = 


1 


(3.34) 


3.3 Crack tip terminating at the interface from material (1) 

This case is shown in Figure 3-4, in which a 1 > 0 , b l = h . 



Figure 3-4: Geometry of the crack terminating at the interface from material 
( 1 ) 

If b l goes to h , the only singular terms in the singular integral equation 
(2.62) are the dominant term — — and some terms coming from the kernel 

f i ~ x \ 

*n(*l’ f i) • Let 

*n (x v t t ) = *l«(Vi)+*u(Vi). ( 3 -35) 


where *iii(^i^i) is the singular term in the case of a crack terminating at the 
interface, and k b n is the bounded term. The singular term is found from 

*n;(*i>'i) = f G~ u (x v t v a)da. (3.36) 

J o 
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where G“ lt is the asymptotic value of G n for a— » °° , — » h , x 1 -> h . Thus, from 

(3.36) and (2.65) it may be shown that 




■I — K h &**i- a *‘dix , 

2 n 1 K 2 + 1 


(3.37) 


^11 ^I2(^~ x l) Gi 3 (/t-jri) 2 

llf ~ (2/i-A: 1 -r 1 ) + (2/z-x 1 -^) 2+ (2/i-x 1 -r 1 ) 3 


(3.38) 


where 


C - 3 " 1 ' 1 l mK 2~ K i 

11 2 m + Kj 2 mKj + 1 


C 12 “ “ 6 : 


C13 = 4- 


i 1 ! 

m = — . 

\h 


Equation (3.38) can also be written in the form 

*u;( x i> f i) = \.-\ c vs&- x \f L - — c \\} * 


(3.39) 


[ 1 

l r 1 -(2/*-jc 1 ) J 


(3.40) 


Therefore, the singular integral equation (2.62) can be written as 
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boimdcd terms 


JCfKj + l) 
4 ^l 


Oi r (Xi) . 


( 3 . 41 ) 


We again define 

Wi)- 




(/l-ai) a i(A-?i)Pi 


= Si('i) H»i(fi). 


( 3 . 42 ) 


Also, define the sectionally holomorphic function 
i f* fiW 


nJ a 1 r l ' 


■dt\ . 


( 3 . 43 ) 


Using equation (3.42) into equation (3.43) and following [37] we have 

g(a{)e nia i g x {h) 


F\( z ) = 


(/i-a^Pifz-a^i sinna^ (h-a{) a i(z-h)P\ sin Ttpj 

+F 01 (z) 


( 3 . 44 ) 


From (3.44) it can be shown that for a 1 <x l < h 


JW = 




F^A-Xj)— - 


{h-a x f\ (Xj-aj) 0 ! (.h-a{) a \ (h-x^ 
g iW 


+^oi(^i) • 


(A-aj) 01 ! 


+^oiC*i) • 


( 3 . 45 ) 


( 3 . 46 ) 


where F 01 (xi) , F^ (atj) are similar to F 0 (xj) in equation (3.8). Substituting from 
(3.45) and (3.46) into (3.41) we obtain 
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g l (a 1 )cotKa l g l (h)cotn^ l 

(h-a ] f\(x l -a l ) a \ (h-a^ih-x^i 


+ [~9^i3(* _x i) 2 — z~ C n( h - x \>fo — ^11 1 * 


dx\ 


-giW 


(h-a^i^h-x^x ji'/ijcpj 




where ^(xj) contains all the bounded terms. From (3.47) we have 

g^ajcotna.! g 1 (h)cotit^ l 

(h-a x f i {h-a^iQi-x^ i 

g iW 

— * 

2^13Pi(Pi + *) — ^nPi - = Vi(*i) • 


By multiplying equation (3.48) by (^-^“i and letting x r —> a 1 


characteristic equation for oij is found to be 

1 


cotnct j = 0 


->°1 = 2- 


Also, by multiplying equation (3.48) by (h-Xj ) Pi and letting x l —>h 
characteristic equation for p x becomes 

g^COtT#! gy(h) 

(h-a^i (h — a 1 ) a i smnP 1 

[ _ 2^'l3Pl(Pl + J = 


which can be written as 


(3.47) 


(3.48) 
, the 

(3.49) 
, the 

(3.50) 
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( 3 . 51 ) 


C0S7tPi — — CijPiCPi + 1 ) — ~ Cn - 0. 


Equation (3.51) has a real root in the interval (0,1). The root is dependent 
on the material constants. If material 2 is stiffer than material 1, the root will 
be less them ~ . But if material 1 is stiffer than material 2, then the root will be 


larger than ■= . The fundamental function of this case is 


^(^1) = 


1 




( 3 . 52 ) 


3.4 Crack tip terminating at the interface from material (2) 

This case is shown in Figure 3-5, in which a^ = h , b 2 > h . 


M 



Figure 3-5: Geometry of the crack terminating at the interface from material 
( 2 ) 

If 02 goes to h , the singular terms in the singular integral equation (2.63) 

are the dominant term — - — and some terms coming from the kernel k 2 2 (x 2 ,t 2 ) . 

h~ x 2 

Let 
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^22 (*2 ’ *2) ~ ^22 i ^ x 2 ’ ^2) ^22 (*2 ’ * 2 ) ' 


(3.53) 


where the singular term for this case, and /i^ 2 ( x 2 > f 2^ i s bounded. 

r oo 

* 22 * (*2 ^ 2 > = G 22 i (* 2 > r 2 > a ) • ( 3 - 54 ) 

Jo 


where G 22l - is the asymptotic value of G 22 for a-» °° , r 2 -» /i , and A . From 
(2.65) and (3.54) it may then show that 



l ”**2-*l m - 1 
2 /M + Kj /MKj+l 


{-l+2(f 2 -A)a} 


{-| + (^ 2 -/t)a}]e- (,2+J£ 2- 2,,)a da . (3.55) 


giving 


G 21 C 22 (^2~^) C 23 (x 2 -h) 

* 221 ^ ~ (t 2 +x 2 -2h) + (f 2 + * 2 - 2*) 2 + (f 2 +x 2 - 2/z) 3 


(3.56) 


where 


C 2 i - - 


G 22 “ _ 


C 23 - 4 
m = — 

H 


l mK 2- K i 3 m- 1 

(a) 

2 m + Kj 2mK2+l 

, m- 1 

1 

(b) 

m-1 
mK^ + 1 

(c) 


(d) 


(3.57) 


The kernel given by equation (3.56) is quite similar to the corresponding kernel 
given by equation (3.38) . Thus, following the procedure leading to (3.49) and 
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1 


(3.58) 


(3.51) we obtain 


* 22 iCV 2 > = [lc 23 (x 2 -h) 2 —-Ci 2 (x 2 -h)-£- + C 2] ] 
2 dx\ 6X2 


l f 2 -( 2h-x£ 

fb 2 i 

[ + JL,, (x 2 ,t 2 ) ] <f> 2 ( f 2) dt 2 + bounded terms 

J h h ~ 


]• 


*(*2+1) r 

a» l 1 


4^ 


<y 2 ( x 2>- 


1 f *2 <t> 2 ( r 2 ) g 2 (h)conux 2 £ 2 (& 2 )cofTtp 2 


lfWy. *2 

— I 

rcj/, 


h *2“*2 Q> 2 -h'^2(x 2 -h) a 2 (b 2 -h) ( h{b 2 -x 2 f , 2 


+^02(- x 2) 


-dt 2 =- 


g 2 w 


1 f b 2 ^2^2) 

71 J A 1 (b 2 - h)^i(x 2 - h) 0 ! sin 7104 

g 2 (h) cotna^ g 2 (b 2 ) cof 7 rp 2 


+F 02 (x 2 ) . 


(b 2 -h)^2(x 2 -h) < h (b 2 -h) c h.{b 2 -x 2 fi>2 
^C 23 (x 2 -h) 2 ^-C 22 (x 2 -h)^+C 2l ] 


dx 2 2 


g 2 00 


= V 2 (*2) 


(b 2 - A)P 2 (x 2 - A) 01 ! 5 in na^ 
g 2 (h)cotn 04 


Q^-hfiQ^-kyh Q^-KfhQ^-xjfr 


g 2 {h) 


(b 2 — h)$2(x 2 — hjh. sin 7104 ^ 


[o^ , 23 ct 2( ct 2 + l) + ^22 cc 2 + ^'2ll “ V2( x 2-* • 


(3.59) 

(3.60) 

(3.61) 


(3.62) 

(3.63) 
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(3.64) 


COt 7Tp 2 = 0 



and 


COS TZCL^ "f - — 1) 4* C^OC^ 4 1 C^| — 0 


(3.65) 


Equation (3.65) is identical to equation (3.51) if the material constants are 
changed. So, we have the same argument as before. If material 1 is stiffer than 
material 2, is less than ^ and if material 2 is stiffer than material 1, c ^ is 

larger than i . The fundamental function in this case is 


w 2 (x 2 ) = 


1 

(x 2 -h) CL i(b 2 -x 2 ) l/2 


(3.66) 


3.5 Crack going through the interface 

This case is shown in Figure 3-6 in which a 1 > 0 , b l = h and a*, = h , b 2 > h . 


'I 



Figure 3-6: Geometry of crack going through the interface 

In this case, the two singular integral equations (2.62) and (2.63) must be 
examined. If b l and a 2 go to the interface, we will have three irregular points 
x=a^ , h and b 2 . 
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The singular terms in the singular integral equation (2.62) are the 

dominant term — — and some other terms coming from the kernels (*i .^i) 

h x i 

and k 12 (x ± , f 2 ) • Similarly, the singular terms in the singular integral equation 

(2.63) are the dominant term — — and some other terms coming from the 

h *2 

kernels ^21 (a^,^) and k^ (x 2 , t 2 ) . The singular terms in the kernels (;tj .fj) and 
k 22 (x 2 ,t 2 ) are sill the same as equations (3.40) and (3.58). By examining the 
singular terms in the kernels k l2 (xj ,t 2 ) and k 21 (x 2 , , let us define 


*12(*l< r 2) = k l2( x l< r 2) + k 12( x l’^ > 

*21 C*2^i) = *21 C*2-'l) + *21 (*2> r l) • 


(a) 

( 3 . 67 ) 

(b) 


where k\ 2 (x 1 ,t 2 ) , (x 2 ,t{) are singular terms, and k b 12 (x v t 2 ) , ^ {x 2 ,t{) are the 

bounded terms. Since 

-00 

*12 C*i« h ) = G i2 (*i > r 2 ’ a ) da ’ ( a ) 

J o 

( 3 . 68 ) 

*21 (*2’fr) = f G 2i (*2>*i.a) da. (b ) 

J o 


where G~ 2 > G 21 are the asymptotic value of G 12 , G 21 for a-» °° , t 2 ~* h , a^-* h 
and a — > h , x 2 -> h , respectively. From (3.68) , (2.65) and Appendix (C) 

we then obtain 
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Ki + 1 


***■/, 


~ m(K 2 + l) ^3 
0 ^ mi^+l 2 


{-^ + (x 2 -h)a}-- 


e -(t 2 -x j)a da ( 

\m{ Kj+l) 


(a) 


( 3 . 69 ) 


2 m + Kj 


{l+2(r 1 -A)o}] 


g“( x 2 -, l) a dct . 


(b) 


or, evaluating the integrals we find 




d n d n Qc x -h) 


(* 2~ x l) (h~ x l> 2 


4 i(Vi) = 


^21 d22(.X2~h) 

(*2"'l) + (^ 2 - f l) 2 ’ 


where 

d _ 3 K i + 1 1 K i + 1 

11 2m+Kj 2/WC2 + 1 

K,+l K,+l 
m + Kj mK^+l 

3m(K2 + l) im(K2 + l) 
^ 21 2 m<2 + 1 + 2 m+Kj 

m( K 2 + I) m(K2 + l) 

** 22 /MK 2 +I m + Kj 

Hi 

m = — . 

H 2 


(a) 


( 3 . 70 ) 

(b) 


(a) 


(b) 


( 3 . 71 ) 

(c) 


(d) 


(e) 
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Equations (3.70) can also be written as 


t J 12 (x 1 ,r 2 ) = {rf 12 (x 1 -/ I )— +d n }— , 

1 2 1 


*21 (*2* r l) ” t ^22^1 ^21 J ' 


(a) 


(b) 


(3.72) 


The two singular integral equations (2.62) and (2.63) then become 

f {- — +*n, (^i» ? i) } < ! > l(^l) < *l + f * 12 ^ 1 ^ 2 ) ^2 

J a l ‘l - *! J h 


+ bounded terms = - 


T *21 O^l^l^l^l+f 2 { r4r + *22 ( ^2- f 2) } H f 2> dt 2 

* dj J h l 2 x 2 


"(K 1 + 1 ) r . . 

-5ST 1 ^ ' 


+ bounded terms = - 


< k 2 +1 ) t . . 
nr” 2 <Ji) - 


(a) 


(b) 


(3.73) 


Substituting equations (3.40), (3.58), and (3.72) into equations (3.73), and 
by using the same technique as before, we will end up with the same equations 
as (3.45), (3.46), (3.60) and (3.61). Also, from the two kernels k\ 2 (x v t 2 ) and 

*21 (* 2 ,^), we obtain 
1 r b 2 $2^2) 


1 rWivv 

h~ x \ 


g 2 (h) 


+F n (x{) 


{b^-h^iQi—x^fh sin na^ 

If*^*. S -^L +f * lW . 

^Jflj^l - -^2 (/i-a 1 )“i(;c 2 -/z)n sinK^ l 


(a) 


(b) 


(3.74) 


where F^O^) and F 21 (r 2 ) are similar to F 0 y(x ; ) in equation (3.8). Then the 
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singular integral equations (3.73) become 
cot raxj g x {h)cotn^ x 


(h-a^P 1 (x x -a x ) a i (A-flj)®! (A-^i) p i 


■f [ ] 




1 


-*iW 


(/* — (A-OTj)^ J/flJlPj 


* 2 <*> 




x 2 5m 71 CX 2 


d -SiW 

[ dii&l d 2 i 1 1 5 

d*2 (h-a x ) a i (x 2 -h)^\ sinn$ x 


(a) 


( 3 . 75 ) 


g 1 {h)cotna 2 g 2 (b 2 )cotn$ 2 


(b 2 — hfii (x 2 — Kfh. (b 2 -h ) a 2 (b 2 -x 2 )$2 


+ 1 2 ^23 (*2 - C 22 ( x 2 ~ + C 2 1 1 


dxt 


[ 


g 2 W 


(b 2 - hfti (x 2 -h )°2 sin na^ 


] =^ 2 (r 2 ) . 


(b) 


where 4 / 1 (r 1 ) , '¥ 2 (x 2 ) contain all the bounded functions. Then equations (3.75) 
can be written as 
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g^a^cot na x g x (h)cotn$ x 

(h-a{ft 1 (h-a^iQi-x^ 1 


*iW 


(A-aj) 0 ! (A -x^i sin TtPj 2 


[-^^(Pi+l^uPi-Cn] 


* 2 W 


(^-A^^-Xj) 0 ^ sinna^ 


[-dl 2 a 2 +dn] = 'Fi(xi) , 


(a) 


(3.76) 


-*iW 


(h-a{) a i (x 2 -h)^isinK^ l 


[“^22^1 ”421 1 


g 2 (h.)cotncL2 g 2 (b 2 )cotn$ 2 


(A 2 -A)p2(x 2 -A) 0t 2 (A 2 -A) a 2(A 2 -x 2 )p2 


* 2 <*> 


(A 2 - A)p2 (x 2 - A)“2 sin 7102 ^ 


[ x^’23 a 2(° t 2 + 1) + C 2 2 a 2 + < ^21 1 


='P 2 (x 2 ). (b) 


It is important to mention that, we just have one irregular point at the 
interface, therefore, 

Pi = «,. (3.77) 

We now multiply equations (3.76)(a) and (3.76)(b) by (*i - <3i) a i and 
(A 2 -x 2 )^2 respectively , and let x 1 — > a x and x^— >b 2 . Then the characteristic 
equations for and p 2 are found to be 


cottux x = 0 

1 

“* a l"2* 

(a) 

cot 7tP 2 = 0 

- fc-|- 

(b) 
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Also, multiply the same equations by (/i-jc^Pi and (x 2 -h)$ 1 respectively, 

and let -tj -> A and x 2 —>h. Then we find 
g,(A) i 

-[ — -CisPiCPi + ll-C^Pi-Cn + CO^TtPi] 




8lW 


(b 2 -h)^2sinn^ 1 


[-^12^1 +<^n 3 - 0 > 


(a) 


(3.79) 


*!<*) 


(h-a^i sinn^ 


l d 22^1 +d 21 3 


80(h) X 

+ s [-C23P 1 (P 1 + l) + C 2 2Pi + C2 1 + C<75nP 1 ]=0. (b) 


(b 2 -h)$2 2 


The solution of interest of (3.79) satisfies 0 < e(pj) < 1 . Thus, since g^h) 
and g 2 (h) are non-zero, the determinant of the coefficients in (3.79) must vanish, 
giving the characteristic equation for Pj as follows 

[ cosn pj + 2 ^ 23 Pl(Pl + ^) + ^22Pl + ^21 3 


[ COS7tp 1 --Ci 3 P 1 (Pi + l)~Ci 2 Pi-Cn ] 


+ [ d 2 2pi+d2i] [— d 12 Pi+dn ] — 0 . 


(3.80) 


It is clear from equation (3.79) that, g^h.) and g 2 (h) are not independent, 
and are related by 


0 h-aj) a i cosP 1 +ic2 3 Pi(Pi + l)+C 2 2Pi+C 2 i 
gl(h)— ^( h >^—^ 2 ^Pl+^21 


(3.81) 


Equation (3.81) will be necessary in order to obtain a unique solution for 
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the singular integral equations (2.62 ) , (2.63) . Equation (3.80) has always a real 
root in the interval (0,1) . The fundamental functions of the singular integral 


equations (2.62) , (2.63) are 


Wl(Xl) (;t 1 -a 1 ) 1/2 (A-x 1 )P i ’ 

(a) 

W2(X2) (x 2 -hfii(b 2 -x 2 ) 1 ' 2 ' 

(b) 
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Chapter 4 

Numerical Procedure 

The solution of the problem depends on the two unknown density 
functions ^ , <J> 2 which can be obtained by solving the two singular integral 
equations (2.62) , (2.63) numerically by using any one of the techniques 
available in [40] or [38 , 39]. In this work the expansion method described in [38 
, 39] is used. In order to solve them let us rewrite these two equations (2.62) 
and (2.63) in the following form 

f [- — ~ + ^n(*i> ? i)] + f ^ yi^ x 1^2) $2^2) ^(2 

*1 J a 2 

= p 1 (x l ) ;a 1 <x 1 <b l , (a) 

(4.1) 

[ ^21 (■*2’ ? l) < i ) l^l)^l + f [7 "*■ ^ 22 ^ x 2 '^ ^2 (*2) 

J a x J a 2 r 2“*2 


= p 2 (x 2 ) ;a 2 <x 2 <b 2 . (b) 


where 

«(K X + 1) T 

P i(*i) = 4jjj— 1 ®i‘ (* 1 ). 

(4.2) 

7t(K2+l) 

P&7> 4 ^—' ( x 2>- 


Normalizing the two singular integral equations (4.1) by using the following 
transformations 
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b r a i b i +a i 


h = — ~ r i + 2 — ’ * V -1 * n z l) 


b r a i b i +a i 


= — 2 — ^ "* 2 — '( a i^ x i^ b i~ l 


(a) 


(b) 


(4.3) 


the two singular integral equations (4.1) may be expressed as 

C+i i C+i 

[- — r + *ll(*l’ r l)W''i)dr 1 + K n {s l f 2 )y 2 (r 2 )dr 2 

J -1 r \~ 5 \ J-1 

-P\{s{) ;-l<^< 1, (a) 


f+i f+i l 

I ^21 ( s 2> r l)Vi( r i)^ / ‘i+ [- — +K 22 {s 2 ,r 2) l¥ 2 (^ 2) ^2 

J-1 J - 1 r 2~ s 2 


= p 2 {s 2 ) ; — 1 < *y 2 < 1 • (b) 


(4.4) 


where 


K ij(*r r j) = *,y (V,) I (‘ J “ 1 > 2 ) ■ 


(4.5) 


and y t (r,) , /*,(£,) are respectively the transformation of <(>•(?•) , p-(x-) (i = 1,2). Since 

(a) 






82^2) 

§2^2) = 7 " Tft " = S2^2> w 2^2) • 


(4.6) 


( f 2 “ a 2) CLl ( b 2 ~ t 2 fa 


(b) 


Then, Xj/^rj) , \y 2 (r 2 ) will be equal to 
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Vi^i) = 


(l+r 1 )“i(l-r 1 )Pr 


F 2 {r 2 ) 

V2(r2) "(l+r 2 )^(l-r 2 )fV 


(a) 


(4.7) 


(b) 


where 


F i( r i) = (r^> ai+, W r i>* 




F 2( r 2) = ( ■ 2 ^ + P2 A 2( r 2) • 
t?2 — ^2 


(a) 


(b) 


(4.8) 


and AjCrj) and A 2 (r 2 ) are the transformation of g^) and g 2 (r 2 ) respectively. 

Then, the two singular integral equations (4.4) may be written as 

f +l F^) i 

I jr[ +K n (s l f l )]dr l 

J _1 (l+r^ia-r/i r 1 -s l 




+1 F 2 (r 2 ) 


+ 1 71 ^77 - r ^ 12 ( 5 1 ,r 2 )dr 2 =F 1 (s 1 );-l< 5l < 1 , (a) 

(l+r 2 )“ 2 (l-r 2 )P 2 


j.+l F i( r i) y f+l F 2( r 2> 

J-l (l+r^Kl-r^Pi 21 ^ ,ri ^ J_! (l+r 2 )“2(l-r 2 )p2 


(4.9) 


[ + K 22 (s 2 s 2 ) ] dr 2 = P 2 (s 2 ) ; -1 < s 2 < 1. 


r 2~ s 2 


(b) 


We may define 
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K, + 1 * 


K 2 +1 _ * 

F 2^ r l)- °0 F 2 ( r 2 > ' 


(a) 


(4.10) 


(b) 


where 


T = _ 


a i’^i ©o 
1-v, ' 


(4.11) 


Then, the two singular integral equations (4.9) can be put in the form 

r+i F i( r i) 1 

J-i (l+r^ifl-r^Pi r \ 


^i^ + lf+l F 2( r 2) 


M2 


i*+i 

Kj + 1 J_! 


(l+r 2 )«2(l-r 2 )p2 


K n {s x f 2 )d r 2 


o/^i) 

= -n — ; -1 < j x < 1 , 

CT o 


(a) 


H2 K 1 +1 

ih* 2 +1 


r +1 ^fri) 

j-l (l+rj)®^!— r x )Pi 


F 2l ( S 2’ r \)^ r \ 


r+i F l^\) 

J-l (l+r 2 )“2(l-r 2 )p2 


C — + ^22 ( s 2' r 2) ] ^ r 2 

r 2~ s 2 


o 2 t (s 2 ) 

= -*-=-==- ; - 1 < 5 2 < 1 . 
°0 


(4.12) 


(b) 


Let us assume that the two functions F* (r t ) and F 2 (r 2 ) are in the form of simple 
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power series such as 


N 


F l( r l)= X a n r v 

n= 0 

(a) 

M 

W-l K'2- 

(b) 


m=0 


(4.13) 


where a , b are the coefficients of these series to be determined, the numbers 
of which are (A/ + 1) and (M+ 1) , respectively. Substituting equation (4.13) into 
equation (4.12) we obtain 


N 


M 


y a„£?i(Ji)+ £ b m E r ? 2 (s 1 )=P 1 (s 1 ) ; 

n = 0 m = 0 


(a) 


N M 

J,a n E n 21 (s 2 )+'£b m E?Z 2 (s 2 ) = P 2 (s 2 ) ; -1<* 2 <1. Cb) 


(4.14) 


n = 0 


m=0 


where 


n 

<l(^l)= f +1 _[— L_+AT 11 ( Jl ,r 1 )]dr 11 

11 1 J_i (l+r^ia-r/i r l -s l 


u, ic, + l f+l r~ 

EZ (j,) = --=-7 I =- K 12 (s^ 2 ) dr 2 - ( b ) 

12 1 M- 2 K l + lJ -l (l+r 2 )«2(l-r 2 )P2 


lt 2 K l + 1 f +1 r l „ , w 

= 7 ; ^ ^-^ 2 i( s 2 A)^i- 

^Kj + l J_! (l+r^^l-r^Pi 


(a) 


(c) 


(4.15) 


m 

> - f 1 ^ X [ -4j- + *22< V 2 > ] *2 . w ) 

J -1 (l+r 2 ) a 2 (l-r 2 ) p 2 r 2 5 2 
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PlisO =-K 


(e) 



P 2 ( s 2) = ~ n 



(f) 


It is obvious that equation (4.14)(a) has (N+l) coefficients a n , and equation 
(4.14)(b) has (M+l) coefficients b m . By choosing the number of the collocation 
points for 5 1 , s 2 between (-1 , +1) , which depend on the number of the unknown 
coefficients a n ,b m , we can obtain a system of linear equations that can be 
solved to give these coefficients. It is clear that, the number of the coefficients 
a n , b m depend on the convergence of the series in equations (4.13). 

Although there is no restriction on the choice of the collocation points, it 
was shown by Kaya [41] that a symmetric distribution with respect to the origin 
considering more points concentrated near the ends seems to help. So, the roots 
of the Chebeychev Polynomial will be used as collocation points. 

To obtain a unique solution for the singular integral equations (2.62) , 
(2.63) additional conditions are needed depending on the crack configuration. 
Each configuration will be investigated separately in the following subsection. 


4.1 Embedded Crack 

Figure 3-1 presents this case in which a 1 > 0 , b l < h , a 2 > h and b 2 < 
Also, the singularity at the end points are 

«i = Pi =02 = 1*2 = 1. (4.16) 

The collocation points for Sj and s 2 are selected as 
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= cos 


71 12 . I c= I 2 JV 

, /v 1. , , ....u ▼ , 

2N 


(a) 


^2/ 


= COS 


jc(2/-1) 

2M 


1 = 1,2 Jfcf. 


(4.17) 

(b) 


where s lt , s^ are the roots of the Chebeychev Polynomial. By substituting 

equations (4.17) into equations (4.14) we obtain the following system of ( N+M ) 

equations for {N+M +2) coefficients a„,b m 
N M 

'y. a n^[\ (^lJfc) + ^m^vi ( 5 ljt) - 1 < s^ < +1 

n=0 m=0 

; k=l,2,.N, (a) 



M 


fe,)+ s t&wr&a) 

m = 0 


^ S^< 1 


1=1,2, ..M. 


(4.18) 


(b) 


where 


*>u)-f — — +A: ii(' s u* r i)l ^ r i 

J -1 . , 

Jl-rf 


^<2 + 1 f +l r 2* „ , 

jL2 ( j u) = — 3~TT I r= *12 (^1* ^ dr 2 > 
1^2 K l + i J-1 , 

Vl-r, 


_ HoK. + l f +l r, 

^21 (^i) = ~ TTTJ I K 21 ( s 2l ’ r 0 dr l ’ 

l i 1 K 2 +W - 1 . 2 

<\-r] 


(a) 


(b) 


(c) 
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1 


(d) 


^22 ( s 2 /) = f 2 — tr — +K 22 ( s 2l ' r 2) 3 dr 2 > 

J -1 , , r 2~ s 2l 


Vi- 2 




— 71 


Oi r (*i*) 



0 2 T (S2j) 

P 2 (s 21> ~~ n t 


(e) 


(f) 


Two extra equations are needed. From the definition of the density functions 
<j>i , <j> 2 , we have the following single valuedness conditions providing the (N + \) th 


, and (M + 1)'* equations : 



S 

(a) 

(4.20) 

if- 

/- — \ 

kT 

II 

O 

(b) 


which can be written, after changing the variables r 1 
equations (4.13) , in the following form 

and t 2 , and 

using 

N 

o, 

n=0 

(a) 


M 

m=0 

(b) 

(4.21) 


where 



(a) 



n 



dr 


l’ 


G m = 




dr. 


(4.22) 


(b) 


The integrations in equations (4.19) and (4.22) can be found in appendix (G). 


4.2 Edge Crack 

Figure 3-2 shows this case in which a 1 = 0 and b±<h . Also, the singularity 
at the end points are 

0,-0 , P,-|. (4.23) 


Since, we have one edge crack, then the collocation points for ^ are 


n(2k-\) . , _ ,. T ... 

s '‘ = CM 2^riT ( y 


(4.24) 


In this case, the number of the collocation points are (N+1) which should be 
equal to the number of the unknown coefficient a n , because the single 
valuedness condition (4.20)(a) is no longer valid. Thus, by substituting equation 
(4.24) into equation (4.14)(a), considering one crack in material (1), we obtain 

the following (N+ 1) equations for (N+ 1) coefficients a n , which are 

N 

X ^1^1*) = P l (*1*) ; - 1 < *1 < 1 ;k=l,2,..(N+l). (4.25) 

n=0 

where 
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Zi' 


(lit)- \* l -p= f-~ +K U (s lk- r t>l dr l 

J-Hl-rj r l - 5 l* 


(a) 


(4.26) 


D , . °i T ( s lk) 

p i( s ik)--* f ■ 


(b) 


The integration in equation (4.26) can be found in Appendix (G) 


4.3 Crack Terminating at the Interface 

Figure 3-4 shows the case in which a { > 0 , b x = h , and the singularity at 
the end points are oc 1 = i and (3 1 , where P 1 can be obtained from equation (3.51). 

Similarly Figure 3-5 shows the case in which a 2 = h,b 2 <°° and the singularity at 
the end points are P 2 = j and , where can be obtained from equation (3.65). 


For the case shown in Figure 3-4, the collocation points and equation (4.14)(a), 
(considering one internal crack in material 1), become 

s lk = cos K ^ , — ; k = 1 , 2 , JV, (4.27) 


2N 


N 


n=0 


(4.28) 


where 


n 


(4.29) 


Equation (4.28) has N-equations with (N+ 1) coefficients a n . The {N+\) th equation 
necessary to obtain a unique solution can be obtained by using the following 
single valuedness condition; 
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» n 


(4.30) 


which can be written in the form 

N 

I^-o. 


(4.31) 


n=0 


where 


n 

f+l r , 

G n = | — i g-dr, 

J _1 (1 +r 1 ) 1 ' 2 (l-r 1 )Pi 


(4.32) 


Similarly, for the case shown in Figure 3-5, the collocation points and equation 
(4.14)(b) (considering one internal crack in material 2) will become 


rc(2/- 1) 


S V- C0S - 2 M 


; /= 1,2 M, 


(4.33) 


M 


m — 0 


(4.34) 


where 


E>2,) = f‘ 

J-l 


(l+r 2 ) <x 2 (l-r 2 ) 1 ' 2 r 2 -^ : 


[ — + f^22^ s 2l' r 2^ ] ^ r 2 • 


(4.35) 


2/ 


Equation (4.34) has M-equations for (M+l) coefficients b . The (Af+ \) ,h equation 


can be obtained by using the following single valuedness condition; 
r b 2 

I ^2 (* 2 ) ^2 = ® ■ 

* n 


(4.36) 


which can be written as 

M 

I i>„c" = o. 

m=0 


(4.37) 
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where 


f +1 r 2 

G m = —dr 2 

J-l (l+r 2 ) a 2 (l-r 2 ) 1 ' 2 


( 4 . 38 ) 


The integrations in equations (4.29), (4.32), (4.35), (4.38), can be found in 
Appendix (G). 


4.4 Crack Going Through the Interface 

Figure 3-4 presents this case in which aj > 0 , b^=h, a 2 =h , b 2 < 00 , and the 


singularity at the end points are 

= P 2 = 2 * Pi = a 2- 


( 4 . 39 ) 


where f5 1 or can be obtained from equation (3.80). In this case equations 
(4.17) and (4.18) will still be the same with 

n 

(*u)= f +1 -j ^ < a ) 

J-l (l+r 1 ) 1/z (l-r 1 )Pi 


u, K, + l i»+l r, 

Ef? 7 (i.r.)= -j „ — Ky 2 (sn c j‘ 2 )dr 2 . (b) 

12 U (i + r )P l( i_ r )i/2 12 


(I+ri) w‘(T-^ K ^' )dr " 


( 4 . 40 ) 


(c) 


m 

" f* l 7i — i 7 ± r +W2>'&*2- «> 

J-l (\+r 2 f\(\-r 2 ) iri r 2 s 2l 
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and / 5 i(5 u ),P 2 ( 5 2/) can be f° un( i ^ equation (4.19)(e , f). 

Equations (4.18) have ( N+M) equations for (N+M+2) coefficients a n , b m . 
So, the two extra conditions we need to obtain a unique solution are 

V ' i )^ f % 2 (' 2)<*'2 = 0 - ( 4 - 41 > 

a 2 



and 

(h-a^i cosn $l + 2^23Pl(Pl + !) + < ^22Pl + ^22 

** w = 


(4.42) 


After chan ging the variable and substituting from equations (4.13), equations 
(4.41) and (4.42) become 

n = 0 m= 0 

^I^( +1 )' ,+/f 3i i m (- i r=o. ( 444 ) 

m=0 m=0 

where 


D _MlK 2 +l / 6 2 -A % 
1 H2 k i + 1 h-a x 

(a) 

R 2 = (42 2 Pl +< ^2l)> 

(b) 

^,k 2 + 1 b 2 -h 

R 3= ,(. > Pl 

3 M2Kj + 1 A-aj 

(4.45) 

(2) a l _p 2 [ COS71 +ic 23 P x (P x + 1) + C 22 P : + C 21 ] , 

(c) 
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(d) 


r n_r l r l . 

1 J-Kl+r^d-r/i ri 


f +1 

' 2 ~L\ 


dr,. 


(l+r 2 )Pi(l-r 2 )^"' 2 


(e) 


The integrations in equations (4.40), and (4.45) are found in Appendix (G). 


4.5 Edge crack going through the interface 

This case is shown in Figure 4-1, in which a l - 0 ,b l =O 2 = h,b 2 >h. Also, 
the singularity at the end points are 


1 


E X ' V 1 

‘ 2 


* / I 


r~h — *h 


Figure 4-1: Geometry of edge crack going through the interface 
a l = ° ■ = \ - Pi = a 2. 


(4.46) 


where pj or may be obtained from (3.80). The collocation points for ^ and s 2 
are 
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s lk = cos 


n(2k-l) 

2(N+l) 


$ 2 / = cos- 


2 M 


; *=1,2,.. 

....(N+l), 

(a) 

; / = 1,2,.... 

,.(M). 

(b) 


( 4 . 47 ) 


Equations (4.14) may then be expressed as 

S S ^m^\ 2 ^ S lk^ =E l( s u) ^ ,...(77+ 1) , (a) 

«= 0 m =0 

( 4 . 48 ) 

N M 

5>„^l(*2/)+X ft m ^ 2 (^2/)=^2/)^=l- 2 * W- (b) 

n=0 m=0 


where 


f " 

^ll 


f +1 r l 

(*i*> = j — 

J-i (l-rj 


)Pi 


[- — +K n (s lk ^ 1 )]dr l , 
r l~ s lk 


E 21 ^2/) “ 


M2 K l + 1 f +1 r l 

(l-r/i 


E 2\( s 2l’ r \)d r \ • 


( 4 . 49 ) 


and £^ 2 (s u ) , E% 2 ( s i /) are tbe same as in equations (4.40)(b,d). 

Now equations (4.48) have (N+M+ 1) equations for (N+M+2) coefficients 
a n ,b m . Since, condition (4.41) is no longer valid, then the only extra condition 
that we need to obtain a unique solution is equation (4.41) or (4.44). Appendix 
(G) contains all the integrations needed. 
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Chapter 5 

Stress Intensity Factors 


5.1 Embedded Crack 

The stress intensity factors at the irregular points a l , , a 2 * ^2 ma Y be 

defined by 


Kai)= bm V 2 (fl 1 -x)o lyy (x, 0 ) , 

X-¥ dj 

(a) 

*0>i) = bm 'h(x-b l )o l (x, 0) , 

(b) 

k(a 2 ) = lim ^(^-x) o 2yy (x,0) , 
x—> a 2 7 

(c) 

k(b 2 ) = lim 'fo(b 2 -x)a 2 (x,0) . 

(d) 


We note that the expression for the stresses o lyj ,(x,0) and o 2yy (x,0) in equations 
(5.1) are for x outside the crack and can be obtained from the two singular 
integral equations (2.62) and (2.63), by observing that these equations give the 


stresses for y = 0 outside as well as inside the cracks. Thus, 
f [ + ^ll( JC ’ r l)] < l > l( ? l)^l + f ^ 12 (x,t 2 ) $2^2) ^2 

J <z 1 t l~ X J a 2 

«(«! + 1) 


4 ^1 


O lyy (x,0) ; 0< *< a x ; b l <x< h. 


f Ar 2 i(x,?i) < j > i(r L )<ifi+ j" [- — - + K22 (x,t 2 ) ] $2^2) ^2 

J a , J a. l 2~ x 


TtOCj + l) 
4|I 2 


<J 2 yy(.x, 0) ; h<x< 02 ; b 2 <x< °°. 


(a) 


(5.2) 


(b) 


Substituting equation (5.2)(a) into equation (5.1)(a) we find 
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{ f* 1 4»i(ri)*i+ f 2 K l2 (x,t 2 )$ 2 (t 2 )dt 2 } . 

J a 1 h~ x J a 2 


(5.3) 


Since, for embedded cracks the kernels K n (x,t x ) , K l2 (x,t 2 ) are bounded , equation 
(5.3) can be written as 

Jt(Kj + !)*-» a, 


K<*i) - ; 


lim ^2^Ja 1 -x 


{ dfj + bounded terms } . 


(5.4) 


By defining the sectionally holomorphic function 

' b i§l ( f i) 




r i~ z 




(5.5) 


and observing that 
4>i(*l)- 




gl(t{)e n ®i 


(f i - «i) a i (&i - ^i) p i (f i - ^i)" 1 (^ - ^i) p i ’ 


(5.6) 


After separating the leading terms at the end points of the cut for 0^ = Pi = 2 * 
equation (5.5) can be expressed in the following form, 

7i i 

Sl(«l)«T Sl(*>l) 


F i(r) = 


(b 1 -a 1 ) 1/2 (z-a 1 ) 1/2 (^i- a i) 1/2 ( z - ft i) 1/2 


+F 01 (z). 


(5.7) 


where F 01 (z) has the same properties as in equation (3.6). By using the Plemelj’s 
formula, equation (5.7) can be written as follows 
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( 5 . 8 ) 


1 f»i +ifti) gj( fl 1 ) 

t x -x 1 (ft 1 -a 1 ) 1 / 2 (a 1 -x) 1/2 


flgl) 


+F01W 


X< < 3 l . 


Substituting equation (5.8) into equation (5.4) we find 

4^ Sifai) 

K ° l) = 7^ Tn Um TifT, 

(K 1 + l)x-»a 1 


:) 1/2 


flgl) 

(& 1 -a 1 ) 1 ' 2 (x-b 1 ) 1/2 


+ bounded terms ] . 


(5.9) 


givmg 


„ N 4 ^1 , 2 

zrTT'i-r-ir^O’ 


Kj + l b^-a^ 


(a) 


Similarly, the stress intensity factor at the other irregular points , <z 2 > b 2 are 
found to be 




**»>- Z7 


4 ^2 . 2 


■g 2 ( a 2>< 


K2 + I ^2“°: 


(b) 

(c) 

(d) 


( 5 . 10 ) 


Equations (5.10) may also be expressed as 
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b i d\ * 

k(a{)= 

(a) 

W 1 )=-^ b -±p-o 0 T F* l (+l), 

(b) 

,b 2 ~ a 2 T * 

^ 2 >= ^~-^a 0 T F 2 (-l), 

(c) 

i m * 

fc(b 2 ) = -^^—£c 0 T F 2 (-bl). 

(d) 


5.2 Edge Crack 

In this case, shown in Figure 3-2, the stresses are bounded at = 0 , but 
have singular behavior at b l . So, the stress intensity factor at b x may be defined 
by 

k(b l )= lim 'i2'l(x-b l )o 1 (x,0) . (5.12) 

x-> b J 

where c lyy (x,0) is the stress for x outside the crack, which can be obtained from 

equation (5.2)(a), by considering one crack, i.e 

r b \ 1 71 ^ + 1 ) 

J 0 [— +^ 1 (V 1 )]<i> 1 (h)^i =-^-o lyy (x,0) . (5.13) 

Substituting equation (5.13) into equation (5.12) and observing that the kernel 
K n (x,ty) is bounded as x-» b l and t 1 -» b 1 , equation (5.12) becomes 
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4u, _ 

k(b ,) = lim 'v fa'Jx-b-, 

n(K l + l) x ^b 1 


fb^ ,(f) 

{ «|) 1 (f 1 ) dt j + bounded terms } . (5.14) 

J 0 t x -x 


The density function <J> 1 for edge crack is given by 
Sl('l) 






(5.15) 


Then, by following Muskhelishvili’s technique, the singular term in equation 


(5.14) can be expressed as 

r b \ <?i(0) 

J 0 


g iC*i) 
ix-b{) iri 


+F 0l {x) . 


(5.16) 


So, the stress intensity factor at the irregular point b x , can be written as 
4 ^l 


Kb i) = 


( K 1 + 1 ) JC — » fcj 


lim ^l2^lx— b x 


{ g - l( — — gl ^ l} -+ bounded terms } . (5.17) 

\b x ) m (x-b^V 


Therefore 

Kh) = -^r<2 g) {b x ) . (5.18) 

T I 

which can be reduced to 

^ 1 ) = -^a 0 r F*(+l) . (5.19) 
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5.3 Crack Terminating at the Interface 

The stresses have singular behavior, at the ends a j and b l = h in Figure 
3-4, and at the ends a 2 = h and b 2 in Figure 3-5. First, consider Figure 3-4 in 
which the stress intensity factor at a l and b l = h may be defined by 


k&l) = lhn 'h(a l -x)o lyy (x,0) , 


x— » a . 


k(b, = h) = lim ^2(x-h)Pi a 2vv (x,0) 
x-> h+0 yy 


(a) 


(b) 


(5.20) 


where the stresses and a 2yy (x,0) can be obtained from equations (5.2)(a) 

and (b) by considering one internal crack, i.e. 

tc(Kj + 1) 


1 [ 7~ 




-a ljy (ac,0);0<j:<a 1 , (a) 

^1 


c h 

K 21 (x^ l (t 1 )dt 1 = 

* n 


71(^2+ 1) 

-——o 2yy (x,0)-,h<x. 


(5.21) 


(b) 


Substituting from (5.21)(a) into (5.20)(a), and observing that the kernel K ll (x,t i ) 

is bounded as x^> a { , t x -* a x , equation (5.20)(a) may be written as 

4(ii 

k ( a 0 = ^--. n lim 


k(Ki + 1) x _> £J 1 


rh <j)i(f^) 

[ dt x + bounded terms ] ; a x <x< h. (5.22) 

» n tl X 


•a x ‘1 


Since, the density function ^ in this case is defined by 

Sl ( 'l) 


(5.23) 


by following Muskhelishvili’s technique [37], the singular term may be 
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(5.24) 


expressed as 

l r Mifti) 


gi( a i) 

(/ J -a 1 )Pi(a 1 -* 1 ) 1 ' 2 


g iW 

(A - a j) 1/2 (x - hfti sin rc|3 j 


+ bounded terms . 


Thus, equation (5.22) becomes 
4 Hl 


*(«!>- 


(K 1+ l) 


lim 


g i(fli) 


*— » a. 


(A-a 1 )Pi(a 1 -jc 1 ) 1/2 




(A-a^Ot-ZOP 


•+ bounded terms ] 


(5.25) 


The stress intensity factor at the end a 1 may then be expressed as 


Ka i)< 


4 ^1 <2 


(Ki + 1) (h-aj Pi 




(5.26) 


Since, by definition 


h-a, , „ „ k. + 1 , 

*i<«i> - (— ) 1Ct|S ' o 0 r f,(-D 


Equation (5.26) would be reduced to 


h-a. 


k(a ] ) = <^- o 0 T (2)^-Pi Fj(— 1) . 


(5.27) 


(5.28) 


Similarly, the stress intensity factor at b^ = h can be obtained by substituting 
equation (5.21)(b) into equation (5.20)(b), i.e. 
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( 5 . 29 ) 


k(b, =h) = — — lim 'i7.(x-hft\ 

ji(k 2 + !)*_>/, 


r h 

[ ; h 

» /j 


<r 


The singular terms of the kernel K 2X (xty) , as x — > h , t x — ► h is 
j . , ^21 d 22 (x-h) 

21 C*-'i> {x-t x ) 2 


- [d 2 2( x - h )fa- d 2l] » h<x - 


Thus, equation (5.29) may be expressed as 

4^2 


Hb x =h) 


rcO^+l)*-* h 


lim '\2{x-h)$\ * 


f n d <5>i(^i) 

[ {d 22 (x-h)—-d 2l } dt x + bounded terms ] ; h 

J a dx t x -x 


<x. 


Since the density function is defined by 

*i('i) 


<t>i('i) = 




By following Muskhelishvili’s technique, we have 
1 f * <t>i('i) . -*iW 


1 f " SW; - ^ 

n Ja 1 t i~ x 1 (h-a l ) 1/2 Qc-h)^isinn^ l 


+ bounded terms 


From equation (5.31) it then follows that 


( 5 . 30 ) 


( 5 . 31 ) 


( 5 . 32 ) 


( 5 . 33 ) 


98 


k{b{)= - ^ lixn ^2{x-hf\[{d 22 {x-h)^--d 2] ] 
v.^ 2 + U.x— » h a* 


-J_ iW 

(h-aJ^Qc-hfii sin npi 


+ bounded terms ] . 


(5.34) 


giving 


... . 4 M2 ,x ^22 Pi"*" ^21 ... 

1 “(Xj + l) (A-a 1 ) 1 ^ 2 smjcp 1 ^ 


(5.35) 


By using 

«■!<*> - c^V 2 ** <j» T ^ ^(+«- 

Equation (5.35) may also be expressed as 

1*2 K l + 1 d 22Pl +rf 21 , h -\Q „ T 
k(b,=h) = : — s (•— =— ) p i % ^(+1). 

1 Kj+l imrtpj 2 u 1 


(5.36) 


(5.37) 


Similarly, in the case of Figure 3-5, it may easily be shown that the stress 
intensity factors at the end points a 2 =h and b 2 are given by 

., ,. Mi 1C 2 +1 ~ d n a 2 +d n b 2 ~ h „ T *. .. 

k{a 2 =h) = : — : (— — )“* a 0 i F 0 (-1) , (5.38) 

L1-) Ko ”f* 1 


\\-2 i^ + i sirnza 2 


b fi 

k(b 2 ) = - yj o 0 T (2)^ -«2 F*(+l). 


(5.39) 
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5.4 Crack going through the interface 

As in the previous section [5.3] the stresses at the end point a i , b 2 have 
singular behavior, and the stress intensity factors at these points are still the 
same as equations (5.28) , (5.39) respectively, which can be put in the following 
form 

K<*l) (2) 1/2 -Pi f\ (-1). (5.40) 

k{b 2 ) = -<\ a 0 T <^~- (2)W-Pi F*(+l). (5.41) 



where ,y) and o^hy) are the components of the stresses at the interface 
which after long manipulation, are found to be 
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fb 2 

+ * i4 (h,t 2 ) <j>2 (^2) d h + bounded terms ] , (a ) 

*h 

(5.43) 

fb 2 

+ *24 0**2) $2 ( f 2) ^2 + bounded terms ] . (b ) 

• h 


where Jt 13 , * 14 , ^23 , *74 are the singular kernels given by 

, (h~h) (h-t 1 ){ 3 y 2 -(h-t 1 ) 2 } 

*13 (.h,t 1 ) = ^[A - — --B — — ], 

(h-tjf+y 1 {(h-trf+y 2 ) 2 


(a) 


(t 2 -h) (t 2 -h){ 3 y 2 -(t 2 -h) 2 } 

*4(^2) = -[" 5 o 5 — T~o ]• 

14 (t 2 -h) 2 +y 2 {(t 2 —h ) 2 + y 2 } 2 


(b) 


, y y{y 1 -'i{h-t x ) 2 } 

k 23 (h,t 1 )=A . 

(h-t^+y 2 {(h-t l ) 2 +y 2 } 2 


(c) (5.44) 


4 .(^ 2 >- g - L..* * a 




(^-hf+y 2 {(t 2 -hf+y 2 ) 2 ’ 


(d) 


m( Kj + l) /n(K 2 + l) 

A = — , 5 = 


miCj + l 


m+Kj 


(e) 


^13 » *14 , *23 » *24 can a ^ so be pat in the following forms 
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1 




-Biy{ 


[(h-tj+iy ) 2 [(h-tj-iy ] 2 


}], (a) 


BA^, 1 




+Al>{ 


1 


[(f 2 -*)+i>] 2 [(f 2 -A)-&T 


;}], (b) 


4 3 = 4 _ |> ,• , , 


-By l r + 5}, (c) 

[(h-tj+iy) 2 [(.h-tj-iy] 2 




-Ay{ 


1 1 
-+- 


l(t 2 -h)+iy] 2 [0 t 2 -h)-iy } 2 


•}• (d) 


Substituting from (5.43) into (5.42), we find 

2| x 2 fh 


k = lim y$ 


;[f ^ 3 ^A)<fi(^)^i 

> * /7 


y ->0 «(* 2 +1 ) 

f ^2 j 

+ Jfc 14 (h,t 2 ) <(>2 (t 2 ) dt 2 + bounded terms , 

Jh 


(a) 


(5.45) 


(5.46) 
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_ 2 |x, r h . 

^ = v ^ O yPl K0c, + l) C L 


y— » 0 7C(K2+ 1) 


r*2 j 

& 24 (^** 2 ) ^>2 (* 2 ) + bounded terms . (b ) 

J h 


By observing that 
<t>l('i) = 


sM 




, (j>2 (^ 2 ) ** 


#2 (* 2 ) 


(t 2 -hfi{b 2 -t 2 )W 


(5.47) 


and following Muskhelishvili’s technique, the singular terms can be written as 

1 f* 4>i( f i) _ gi{h)e~2'$\ 

■iy 1 (A-a^iyPi smjtjij 


l f* 4>i( ? i) _ A . . _ 

- — — —dt\ = s + bounded terms . 

nj^h-tj + iy “ 


(5.48) 




If* 4>i(*l) . «i(*)«2 ip i , 

aty = - + bounded terms 


u 


ui < ~ 

£2j (h-t^-iy {h-a^) a \y^\ sinnfy^ 


(5.49) 


1 f b 2 <t>2 ( f 2> #2 2 , Pl 


1 f°2 T2 \*2/ 
TtJfc (h-h) + r 


d?2 - 


h ( { 2 ~h) + iy z (fc 2 -A)P 2 yPi jmjcpj 


+ bounded terms 


(5.50) 


1 f *2 $2 (*2) 


il 


- ^2 — 


5 2 (/i)e2 ip i 


h ( f 2 ~h)— iy (ft 2 -/i)P2yPi 


+ bounded terms . 


(5.51) 


From (5.48)-(5.51) it follows that 
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nL L (h- 


1 r h i 1 

-t z Wh )d h 


Qi-tJ + iy (h-tj-iy 


gi(h) 




■+ bounded terms 


*>{ T 

nJ a t [(h-tj + iy) 2 [{h- h )-iy] 2 


d If* 1 
( 


(h-t{) + iy (h-tj-iy 


-}4>i(/i)^/i 


Pl^lW 


(A- fll )°iyPi rinlpj 


-+ bounded terms 


ir ! (_ 

TW* Oo- 


1 f*2 , 1 i , . . . . 

- + ~ — t ~ } §2 ^t) dt 2 


h {t 2 -h) + iy {t 2 -h)-iy 


g 2® 


(b 2 -/t)P2>Pi sin^Pj 


-+ bounded terms 


if\i — I 

nJh 


h [(?2 -/l ) + I >] 2 [ ( r 2 — *) — O’ ] 2 


7~2 ^ ^2 ^2^ dt 2 


d i e b 2 . 1 1 . . . , . 

y dyn\ h (t 2 -h) + iy + (t 2 -h)-iy 2 ^ * 2 




(6 2 -*)Mi j/nfPj 


-+ bounded terms . 


(5.52) 


(5.53) 


(5.54) 


(5.55) 


104 


If *[_! 

KJ ai (A-?1 )+iy 


i 

(h-t^-iy 


] ^1 (f\)dt^ 


g iW 


(/i-a^iyPi co^5p 1 


-+ bounded terms . 


( 5 . 56 ) 


If* — 


! [(/»— ri>+/y ] 2 [(A-r^-fy ] 2 


} <t»i 


dlc h . i 
y--\ [ 


dynj ai (h-tj+iy (h-tj-iy 


H<t>i('i)^i 


- Pi Si (A) 


(A-a^i/i cos^ l 


■+ bounded terms . 


( 5 . 57 ) 




82(b) 


(b 2 -hfzy$\ cos^p 1 


-+ bounded terms . 


( 5 . 58 ) 


If* 2 ( 1 

r (h-fo+i 


1 


— — } <J ) 2 (t 2 ) dt 2 


h [( t 2 -h) + iyY [( t 2 -h)-iy ] 


y dytS h l (t 2 -h) + iy {t 2 -h)-iy ] ^ l{h)dt2 


-Pi g 2 (b) 


{b 2 -hfiy$\ coi^p i 


-+ bounded terms . 


( 5 . 59 ) 


The two components of the stress intensity factor k x , k^ are then found to be 
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, .. -^2 B r A 5, 

fc. = lim jPi [(-+-) 

>,_>() k 2 +1 2 2 (A _ a jojypj S i n ll p 


Pl^lW J? A £2W 

(A-aj^yPi riiijpj 2 2 

(a) 

(5.60) 


(-Pl)SlW B A ?2( A > 

4 * ( ) 

co^^Pj 2 2 {b 2 -h)hy^\ cos^ x 


PiftW 

+A + bounded terms , 


(*2-*)Mi «»?Pi 


fc„, = lim 


* lW 


** y~* oK 2 +1 2 2 (A-a^i/i pj 


-A 


(“Pi)ftW 


+ bounded terms . 


{b 2 -hf 2/1 cos* Pi 


(b) 


These can be written as 


-2^2 

^ [{A + (l-2p 1 )B} 

2(K 2 +l)«rt|P 1 


fig) 


+ {-5-(l-2p 1 )A} 


(b 2 -/i)p2 ’ 


(a) 

(5.61) 
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xy 


2^ 

2(K 2 +l)co^p i 


[{A-(l-2p 1 )fi} 


g i(A) 


+ {5-(l-2p 1 )A} 


g 2 W ] 


(b) 


Since, by definition 

h-a 1 „ _ k, + 1 „ 

Si (*) - <— )“' * P > °o r C + l) • 

6,-/l n _ K, + l * 

fc(*) = (■ -V )tt2 2 °° Tj^ F 2 (_1) - 


(a) 


(5.62) 


(b) 


where p 1 = CX 2 and a 1 = P 2 = 1/2 , equations (5.61) may also be expressed as 

-1 1 




1 C n 


sin\^) m+1 


[{A + (1-2P 1 )5}( 


*~«l B K l + 1 
/ ^ K2+1 


— fJ(+ 1) + { -B-(l -2Pj)A } (^)PiFj(-l)] , 
M-i 1 


(a) 

(5.63) 



1 

x (2)1/2 + 2 


[{ A— (1— 2Pj)B }(— j— )^i 


4 * 1 
KjH" 1 


— F;(+l)+{5-(l-2P 1 )A}(^)PiFj(-l)]. (b) 

M-l 1 


l being the total crack length (l=b 2 -a l ) . 


107 



5.5 Edge crack going through the interface 

In this case, as in section [5.4] the stresses at the end points b 2 and 
b l =a 2 =h are unbounded (see Figure 4-1), but the stresses at <^=0 is bounded. 
Equations (5.39), (5.61) (a,b) are still valid for stress intensity factor at the crack 
tip b 2 and at the interface b l =a 2 =h respectively, provided we let a±=Q , 
a i = 0 . Pi - <* 2 . and use the definition of g 1 (h) and g 2 ( h ) in equation (5.62). It may 
then be shown that 

k(b 2 ) = -W7 o 0 T V (A_) (2)^i F* (+1) . (5.64) 


- W* °o r { A-Kl -2p,)B ) <y)Pl ^1 

sin -Pi ( 2 ) Pl 1 K 2 + 1 


£ { ' ~ BH1 ' 2PlM ' <565) 


S7 i (^) p ' • 


(5.66) 


where / is the total crack length. 
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Chapter 6 

RESULTS AND DISCUSSION 

As examples for the thermal shock crack problem two material 
combinations are considered. Material pair A corresponds to an austenitic 
stainless steel layer (material 1) and a ferritic steel base (material 2) which have 
the same elastic constants but different thermal properties [1]. Material pair B 
represents a ceramic coating (material 1) on a ferritic steel base (material 2) in 
which both the mechanical and the thermal properties are different. Since the 
problem is formulated in terms of dimensionless quantities, it is sufficient to 
consider only the ratios between the properties of the two material combinations 
shown in Table 6-1. 


6.1 Model I 


Figures 6-1 - 6-10 show the normalized transient temperature and 
thermal stresses distribution in material pair A given by equations (2.17) and 
(2.35) for t 0 = 0 , and (2.26) and (2.35) for x 0 > 0 . The temperature plots 
represent the ratio O(x*,x)/Q 0 (0=© 1 for 0 < x < h , and 0=0 2 for h < x < <») as a 
function of non-dimensional distance x*=x/h for various values of non- 


dimensional time (Fourier Number) x and x 0 defined by 
tD\ 

x = l r ’ X ° = 1F' 


( 6 . 1 ) 


where t 0 is the actual duration time of the cooling ramp on the surface (see 
Figure 2-3(b)). The non-dimensional thermal stresses Oyy/G 0 T are plotted in the 
same way, where 


T «i £ i 0 o 

'O — mT 


( 6 . 2 ) 
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From Figures 6-1 - 6-10 it is clear that, the derivative of the temperature 


distribution, that is, the heat flux and the thermal stresses (o w =o lyy for 
0<x<h, Cyy = G 2 yy for h < x < °°) are discontinuous at the interface. The 
discontinuity in the slope of the temperature, i.e., the heat flux, is due to the 
difference in the thermal conductivity of the two materials (see boundary 
condition (2.5)(c) or (2.19)(c)). Also, the discontinuity in the thermal stresses 
results from the difference in the coefficient of thermal expansion across the 
interface which is given by (see equations (2.33) or (2.35)) 

c 2yy (h,t) c C 1 E 1 Q 1 l-v l cq 


c, „,(/!, r) „ l-v, 

lw a’l 


= — = 0.75 


(6.3) 


since, there is no discontinuity in the temperature field itself Q^(h,t ) = & 2 (h,t) and 
E^=E 2 , V!=v 2 in material pair A. The influence of the ramp cooling x 0 on the 
transient temperature and thermal stresses, particularly for small values of 
time, are also shown in Figures 6-1 - 6-10. As the rate of cooling surface 
temperature (t 0 ) increases, the transient thermal stresses decrease. 

The normalized stress intensity factors (SIF) for various crack geometries, 
subjected to thermal history shown in Figures 6-1 - 6-10, in material pair A are 
given in Figures 6-11 - 6-23. Figure 6-11 shows the variation of normalized SIF, 
for various edge crack lengths in material 1 (02 = 0 , b 1 < h) and for 
T 0 =0 as a function of Fourier Number x . Note that, since the material pair A 
have the same mechanical properties, i.e. E±=E 2 , V 2 =v 2 , the normalized SIF 
would approach the uniformly loaded half plane solution of 1.1215 as x — > The 

influence of the cooling rate x 0 on the surface x = 0 for a fixed edge crack b l /h = 0.5 
is shown in Figure 6-12. It is clear that, for small time x, this influence could be 
quite considerable, and the normalized SIF decrease as x 0 increases. Figure 6-13 
shows the effect of x 0 on the normalized SIF, k(b l )/GQ r 'll l , for an edge crack 
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terminating at the interface (aj= 0 , b^ = h). Since, £ 1 =£ 2 .v 1 =v 2 in material pair 
A, the stress state at the crack tip b^=h has a square root singularity and the 
normalized SIF would again approach the uniformly loaded half plane solution 
of 1.1215 as x -» <». 

Figures 6-14 - 6-15 show the normalized SIF for an under-clad crack at 
the crack tips a 2 ,b 2 ( a 2 =h , h<b 2 <°° ) , k(a 2 )lGQ'h 2 l2 , k(b 2 )/c 0 T -il 2 / 2 , for 
various crack length and for T o =0 as a function of time (Fourier Number). Since 
E l =E 2 , v 1 =v 2 , the singularity at the crack tip a 2 =h is the square root. The 
asymptotic values of the normalized SIF as x -» » is depicted in the same 
Figures. The influence of x 0 on the normalized SIF for an under-clad crack with 
fixed length l 2 /h= 1.0 is shown in Figures 6-16 - 6-17. It appears from the 
Figures that for small values of time, this influence could be quite significant. 

The results of an edge crack crossing the interface are shown in Figures 
6-18 - 6-19. In Figure 6-18 the normalized SIF at the crack tip b 2 , ( a l =0 , 
b^=a 2 =h , h <b 2 < °°), k(b 2 )lc 0 T ^l , for various crack lengths and for x 0 =0 are 
shown as a function of time. Note that, the tensile and shear stresses at the 
interface (^ 1 =a 2 =/i) are bounded, since the mechanical properties of the 
material pair A are the same, i.e. E^=E 2 , v 1 =v 2 . Figure 6-19 shows the effect of 
x 0 on the normalized SIF for fixed crack length l/h= 4.0. The limiting values of 
the normalized SIF as x — » are also shown in Figures 6-18 and 6-19. 

Figures 6-20 - 6-23 show the results of an internal crack crossing the 
interface. The transient normalized SIFs at the crack tips a l and b 2 , (a l > 0 , 
b l = h=a 2 , b 2 < °°) , k{a^)la<f^ll2 , A:(6 2 )/o 0 r V//2 , for various crack length are 
given in Figures 6-20 - 6-21. There is no stress singularity at the interface, since 
E x =E 2 , v 1 =v 2 . Figures 6-22 - 6-23 show the influence of the ramp cooling as 
measured by x 0 on the normalized SIF at the crack tips a x and b 2 for a fixed 
crack length a 1 /A= 0.2 , b 2 /h= 2.0. 
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The results for the material pair B are given in Figures 6-24 - 6-38. The 
normalized transient temperature and thermal stress distributions for various 
ramp durations (Tq) are shown in Figures 6-24 - 6-29. Again, the discontinuity in 
the thermal stresses at the interface (x=h) is given by (see equations (2.33) or 
(2.35)) 


0 2 yy(hj) 


c^£ 2 0 2 1-Vj 
1-v, 




= 1.4019 


a ’l £ l 


(6.4) 


Note that in this case ,E 2 *E l , v 2 =Vj , and 0 2 =0 1 at the interface. 

Figure 6-30 shows the normalized SIF for an edge crack in material 1 
(^ = 0 , < h) for various crack lengths and for t 0 =0. Since the materials in pair 

B have different mechanical properties and since £ 1 > E 2 , the normalized SIF 
would approach infinity as b^-th, but for very small value of b ly it would 
approach to the half plane uniform load solution, 1.1215. Figure 6-31 shows the 
effect of x 0 for a fixed edge crack b 1 /h = 0.5. The results for an edge crack 
terminating at the interface (a^O , b^h) are shown in Figure 6-32 for different 
values of x 0 . In this case the singularity at the crack tip b 1 =h is pj = 0.552538. 

Figures 6-33 - 6-38 show the results of an edge crack going through the 
interface (a L =0 , b^=a 2 =h , h < b 2 < »). Figures 6-33 - 6-35 give the normalized 
SIFs at the interface and at the crack tip b 2 , k x /G 0 T fi i, k xy /o 0 T fi^, k(b 2 )/G 0 T ^~l, for 
various crack lengths and for t 0 =0. The values of k x and k^ control the tensile 
and shear cleavage stresses at the interface, C x=h, y > 0) (see equations (5.42) 
and (5.43)), where the singularity at the interface is Pj =02=0.01872238. The 
influence of t 0 on the normalized SIF is also shown in Figures 6-36 - 6-38 for a 
fixed crack length b 2 /h- 2.0. 

If the transient thermal stresses are sufficient to cause yielding through 
the thickness of the clad, the stress intensity factors for an under-clad crack can 
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be analyzed by using a simple plastic strip model which assumes that the 
yielding in the clad is restricted to a very thin layer in the plane of the crack. 
Thus, the solution can be obtained by changing the boundary condition (2.49)(b) 
as 

o lxy (*,0) = - o^(x)+o F ; 0 < x < h , (6.5) 

o 2yy (x,0) = - o 2 t (x) ; h<x<b 2 . 

where Op is the average plastic flow stress. The normalized stress intensity 

factor at the crack tip b 2 is given by 
k(bj) * Op „ 

— — =-w + -fW- ( 6 «> 

o 0 T 'll 2 l2 °o 

where k*(b 2 ) is the stress intensity factor for an edge crack going through the 
interface subjected to - o^Ot) (0 < x < h) and - o 2 {x) (h < x < b 2 ), and k^(b 2 ) is 
the stress intensity factor for an edge crack crossing the interface subjected to 
o F (0 < x < h). The values of k*(b 2 ) are given in Table 6-2, which shows the effect 
of clad yielding on the stress intensity factors. The SIFs decrease as the crack 
length increases. The comparison between the values of stress intensity factors 
as T — » “ for yielded clad and elastic clad for different values of crack length and 
Op/oJ are shown in Table 6-3. It appears that, as oj becomes larger, yielding 
in the clad will result in a much higher stress intensity factor at b 2 . 

6.2 Model II 

To verify the computer program, the SIF for an edge crack in a 
homogeneous strip (£^=£2 > v i~ v 2 > X = 0) under uniform load is calculated and 
compared with the results obtained by Kaya [41]. The results can be found in 
Table 6-4. 

The problem of a crack terminating or going through the interface in two 
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bonded half planes can be obtained by letting the crack length very small 
compared to h l and h 2 . The normalized SIF for these two cases are calculated for 
material pair Aluminum and Epoxy and compared with the results obtained by 
Kaya [41] and Erdogan and Biricikogh [29], which can be found in Tables 6-5 
and 6-6. 

To see the effect of the elastic foundation on the SIF for a homogeneous 
strip, the normalized SIF for x > 0 (which is related to RJL through equations 
(2.106) and (2.116)) are calculated and compared with the results for x=0, (see 
Table 6-4). Also, the effect of the elastic foundation x on the normalized SIF for 
a homogeneous strip with a fixed edge crack b x /L=0.5 are shown in Table 6-7. In 
this Table, the limiting case x=0 means that the homogeneous strip has no 
elastic foundation. The case x =0 ° correspond to uniformly loaded strip of 
thickness 2 L having symmetrically located coplanar edge cracks of depth b x on 
both surfaces. Note that in this case u(L,y)=0 , c xy (L,y)=0 , -oo<y<oo . As 
expected, the SIF is decreased for increasing x- 

To show the comparison between the homogeneous strip on an elastic 
foundation x ( which is related to RJL through the equation (2.106)) and a 
homogeneous cylinder the normalized SIFs for an edge crack under uniform 
tension are calculated and tabulated in Table 6-8 (see Nied [42] for cylinder 
results for /?,./£,= 9.0) . Also, the comparison between the transient temperature, 
thermal stress distribution and normalized SIF for the homogeneous strip on an 
elastic foundation and a hollow cylinder obtained from Nied [43] for R t /L= 9.0 , 
are shown in Tables 6-9, 6-10, and 6-11, respectively. 

Figures 6-39 - 6-62 show the normalized transient temperature and the 
thermal stress distribution in material pair A given by equations (2.80) and 
(2.100) for x 0 =0, and equations (2.86) and (2.100) for t 0 > 0. The normalization of 
the temperature and the stress is similar to in Model I, varying with the 
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nondimensional distance x*=x/h 1 and nondimensional time (Fourier Number) 
2 

x=tD l /h v The results are given for different values of h 2 /h 1 = 3.0,9.0,24.0, and 
different values of cooling rate as measured by t 0 = t 0 D l /h v where t 0 is the actual 
duration of the ramp cooling. The discontinuity in the slope of the temperature 
as well as the thermal stress have already been discussed in the previous 
section [6.1] for Model I. It is clear from the Figures that, at any instant in time 
and for fixed h 2 /h l , the transient thermal stresses decrease as t 0 increases. 

The normalized SIFs for various crack geometries and different values of 
t 0 and h 2 /h l =3.0,9.0,24.0 in material pair A are shown in Figures 6-63 - 6-107. 
These results correspond to the thermal history shown in Figures 6-39 - 6-62. 
The stiffness of the elastic foundation % are taken to be function of RJL through 
the equation (2.106), where R i is the inner radius of the cylinder and L is the 
total thickness (h 1 +h 2 ). All the results are calculated based on RJL =9.0 , which 
means that the thickness L=h l +h 2 is fixed, but the ratio h 2 /h l is varied. 

The results of the normalized SIF for an edge crack in material 1 (^=0 , 

, Kb l )/c 0 T <T l , are shown in Figures 6-63 - 6-74. The variation of 
normalized SIF with time T, for various edge crack lengths for h 2 /h l = 3.0 , 
t o =0.0,6.0 , h 2 /h l = 9.0 , T o =0.0,10.0 , and h 2 /h i =24.0 , *0=0.0, 20.0 are shown in 
Figures 6-63 - 6-64, Figures 6-67 - 6-68, and Figures 6-71 - 6-72, respectively. 
Also, Figures 6-65 - 6-66, Figures 6-69 - 6-70, and Figures 6-73 - 6-74 show the 
influence of the cooling rate x 0 for a fixed edge crack and h 2 /h l =3.0, 
^//j^O.2,0.9 , h 2 /h l =9.0, bjh x = 0.2, 0.9 , and ^^ = 24.0, 0.2, 0.9 , 

respectively. It is clear that, as the cooling rate or x 0 increases, the maximum 
normalized SIF decrease. Noting that, the total thickness L of the clad and base 
material is fixed, and as h 2 lh l increases the clad thickness h x decreases with 
respect to h 2 . So, by comparing Figures 6-63, 6-67 and 6-71 in which * o =0.0 and 
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h 2 /h l =3.0,9.0,24.0 respectively, and by normalizing the crack length by L, i.e. 
l x /L instead of / 1 //z 1 , it can be seen that the maximum SIF for fixed edge crack 
length (/j/L) would decrease as h 2 /h l is increased, (see [1]). In the case of broken 
clad the influence of x 0 for h 2 /h l =3.0,9.0,24.0 are shown in Figures 6-75 - 6-77. 

The normalized SIF for under-clad crack at the crack tips a 2 ,i> 2 (<a 2 = /i 1 , 
b 2 >h ] ) , kia 2 )/a 0 T ^l 2 /2 , k(b 2 )/a 0 T ^l 2 /2 are shown in Figures 6-78 - 6-95. The 
normalized SIF at the crack tips, for various crack length for h 2 /h 1 =9.0 and 
x 0 = 0.0, 10.0 are shown in Figures 6-80 - 6-83. It is apparent that, initially the 
two crack tips have negative SIF, and then reach maximum values at different 
times. Also, it is clear that, the maximum normalized SIF at the crack tip b 2 is 
dropping faster than that at the crack tip just touching the interface (a 2 -h ] ). 
The slope discontinuity in Figure 6-83 is due to the ramp function and takes 
place x=x 0 =10.0. Figures 6-88 - 6-91 show the normalized SIF at the crack tips 
for A 2 //ij =24.0, x 0 = 0.0, 20.0. The transient behavior of the normalized SIFs at the 
crack tips a 2 =h^ and b 2 for / 2 //ij=9.0 and t q =0.0,20.0 is shown in Figure 6-88 
and Figure 6-90. The results of the normalized SIF at the crack tip b 2 for various 
crack length for h 2 /h l = 24.0 and x 0 = 0.0, 20.0 are also shown in Figure 6-89 and 
Figure 6-91. Again the slope discontinuity which appears in Figure 6-91 is due 
to the ramp cooling and takes place at x=x o =20.0 . The influence of the cooling 
rate or x 0 on the behavior of the normalized SIF for fixed length of under-clad 
crack and for 0, / 2 / ^ = 0.004 , /t 2 //» 1 = 9.0, l 2 /h^ = \Si ,3.0 , and A 2 //i 1 = 24.0, 

=0.5, 4.0 , at the crack tips a 2 ,b 2 are shown in Figures 6-78 - 6-79, Figures 
6-84 - 6-87, and Figures 6-92 - 6-95, respectively. 

The transient normalized SIFs for an edge crack going through the 
interface (a^O , b l =a 2 =h l , b 2 > hj) are plotted in Figures 6-96 - 6-107. Figures 
6-96 - 6-97, Figures 6-100 - 6-101, and Figures 6-104 - 6-105 show the variation 
as a function of time x for various crack length and for h 2 /h l = 3.0, x 0 =0.0,6.0 , 
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h 2 /h l = 9.0, x G =0.0,10.0 , and h 2 /h l = 24.0, x o =0.0,20.0, respectively. The influence 
of x 0 on the normalized SIF for fixed crack length for h 2 /h 1 = 3.0 > b 2 / h 1 = \.2,2.0 , 
h 2 /h l =9.0, b 2 /h 1 = l.5,4-0 , and h 2 /h l = 24.0, b 2 /h l = 1.5 ,5.0, are shown in Figures 
6-98 - 6-99, Figures 6-102 - 6-103, and 6-106 - 6-107, respectively. It is seen that 
for any given crack length the normalized SIF increases until it reaches a 
maximum value and then decreases as x increases. Also, as the cooling rate or x 0 
increases the maximum value of the normalized SIF decreases for any fixed 
crack length. 

Tables 6-12 - 6-14 give the maximum of normalized SIFs for an under-clad 
crack (a 2 =h 1 , b 2 > , and their comparison with the results obtained by Nied 

[1], Similarly Tables 6-15 - 6-17 give the maximum of normalized SIFs for an 
edge crack crossing the interface (aj=0 , b 1 =a 2 =h 1 , b 2 > h lt x 0 =0.0) and their 
comparison with the results of Nied [1]. Note that, in this study the analytical 
solution of the temperature, thermal stresses and the crack problem is based on 
RJL being very large. In Nied’s work [1], RJL- 9.0 which still corresponds to a 
relatively thick cylinder. 

The normalized transient temperature and thermal stress distributions in 
material pair B for h 2 /h l =9.0 as given by equations (2.80) and (2.100) for x o =0.0 
, and equations (2.86) and (2.100) for x 0 > 0 , are shown in Figures 6-108 - 6-115. 
The stress discontinuity at the interface is due to the difference in Young’s 
modulus and coefficient of thermal expansion (£j * E 2 , a’j * a^). 

Figures 6-116 - 6-117 represent the transient normalized SIF for various 
lengths of an edge crack (^=0 , b 1 < hj) for x 0 =0.0,10.0. It is clear that, since 
E l > E 2 , the maximum normalized SIF would increase to °° as b 1 approaches the 
interface. The effect of x 0 on the normalized SIF for an edge crack having a fixed 
length b l /h 1 = 0.2 is shown in Figure 6-118. The influence of x 0 on the normalized 
SIF for broken clad (<^ = 0 , b l = h ] ) is shown in Figure 6-119, in this case the 
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singularity at the crack tip b l =h 1 is equal to =0.552538 . As expected, the 
singularity is greater than 1/2 because E l >E 2 . 

In case of under-clad crack for material pair B ( a 2 =h 1 , b 2 > h{) , the 
normalized SIF k(a 2 ) / a 0 T (l 2 / Ifh. , k(b 2 )/o 0 T 'Jl 2 /2 , at the crack tips a 2 ,b 2 for 
T 0 =0.0,10.0 are shown in Figures 6-120 - 6-121, and Figures 6-122 - 6-123, 
respectively. Since the crack touches the interface from material (2) and since 
> E 2 , the singularity is less than 1/2 , and is equal to 02 = 0.4512416. The 
effect of the cooling rate or x 0 on the normalized SIFs at the crack tips a 2 =h l and 
b 2 for a fixed crack length l 2 lh^ = 1.0 is shown in Figures 6-124 - 6-125. 

The normalized SIFs for an edge crack crossing the interface (aj=0 , 
bj=a 2 =/ti , b 2 > h]) in material pair B, for t 0 =0.0 and 10.0 are shown in Figures 
6-126 - 6-128, and Figures 6-129 - 6-131, respectively. The normalized SIFs 
k x /c 0 T , ky/o control the tensile and shear cleavage stresses at the 
interface. The singularity at the interface is equal to Pi = ot 2 = 0.01 872238 . The 
influence of t 0 on the normalized SIFs at the interface x=h^ , y > 0 , and the 
crack tip b 2 for a fixed crack length b 2 /h^ = 1.5 , is shown in Figures 6-132 - 
6-134. 

6.3 Conclusions 

Very often in transient thermal stress analysis it is assumed that the 
relevant thermal boundary condition is a step change in temperature which is 
used as a model to describe the sudden cooling or heating at the boundary. In 
the case of cladded pressure vessels, the temperature of the inner wall of the 
cylinder is suddenly brought down to the temperature of the cooling liquid (from 
288° c to 20° c). A simple calculation would show that under this idealized 
thermal boundary condition the maximum thermal stress in the clad would 
exceed the corresponding yield strength of the material (which, for austenitic 
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steel, is 390 MN/m 2 ). It is therefore, clear that in this case the unit step function 
in temperature is not a realistic representation of the actual boundary 
condition. By using a more realistic thermal boundary condition such as ramp 
function with cooling rate measured by the ramp duration x 0 , the peak values of 
the thermal stresses as well as the maximum stress intensity factors are shown 
to reduce below the corresponding values given by the step function 
temperature boundary condition. 

In Model I it is seen that the stress intensity factors are monotonically 
increasing functions of time. This implies that in this case the fracture process 
would be unstable. On the other hand in the composite plate represented by 
Model II, invariably the stress intensity factors increase, going through a 
maximum and then tend to their steady state values asymptotically. In some 
cases the under-clad crack may have initially negative stress intensity factors. 
In such cases it should be understood that the negative stress intensity factors 
may have a meaning only when they are used in a superposition with other 
loading conditions in such a way that the resulting stress intensity factor is 
positive. 

The transient thermal stresses are strongly dependent on the material 
properties. Thus, the steady state thermal stresses and consequently the 
corresponding stress intensity factors for material pairs A and B are found to be 
opposite in sign. 

The results show that the initial flaw (the edge or the internal crack) in a 
coated medium subjected to transient thermal loading generally tend to 
propagate towards the interface. Since the power of singularity for a crack 
terminating at the interface is not 1/2, the propagating crack would slow down 
and possibly be arrested if it is located in the less stiff side of the interface and 
would grow faster if it is in the stiffer material. The process of the propagating 
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crack to cross the interface would be controlled by the relative magnitudes of 
material toughness and energy release rate. 

The effect of the clad yielding on the stress intensity factors for an under- 
clad crack in Model I was investigated and it was shown in Table 6-3 that for a 
fixed crack length the stress intensity factor is higher in the yielded clad case 
than that for the elastic clad. Also, as the thermal stresses increase yielding in 
the clad would cause an increase in the stress intensity factor. 

As expected, it was found that an increase in the stiffness (x) of the elastic 
foundation would cause a decrease in the stress intensity factors as it increased. 

The results for the beam on an elastic foundation showed very good 
agreement with that obtained from the axisymmetric elasticity solution for a 
thick-walled cylinder, indicating that the model can be used quite satisfactory to 
study the highly complicated problem of a composite cylinder containing a 
circumferential crack. 

6.4 Suggestions for Future Research 

The way by which the crack problem in a composite beam on an elastic 
foundation was formulated in this study would benefit the future research as an 
extension of our work in the following areas : 

i. The cracking of layered materials perpendicular to the interface 
with a stress free boundary (x = 0) under residual and transient 
thermal stresses. 

ii. The cracking of fully constrained layered materials perpendicular 
to the interface under residual and transient thermal stresses. 

Other interesting problems can also be done in the future such as : 

iii. The cracking of layered materials along the interface with free and 
fully constrained boundaries under residual and transient thermal 
stresses. 
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iv. The cracking of layered materials perpendicular to and along the 
interface under residual and transient thermal stresses. 

v. The "free-end" problem in bonded layers under thermal cycling. 
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Table 6-1: Properties of materials pairs used in this work 




i 

a\ / a\ 

V*i 

V 2> V 1 

material pair A 

3 

3 

0.75 

1 

1 

material pair B 

3.385 

4.070 

2 . 2939 

0.6111 

1 


Table 6-2: Normalized stress intensity factor for sin embedded crack (equation 
(6.6)), obtained from the plastic strip model, x 0 =0.0, 

o 0 T =-ct\E 1 © 0 /l -v lf for Mqdel I. (Material pair A) 


l 2 /h 

* 

A: 1 (6 2 ) as t — ► oo 

w 

0.025 

9.8776 

- 8.9560 

1.0 

1 . 9008 

- 0.8354 

4.0 

1 . 4046 

- 0.2598 

6.5 

1 . 3450 

- 0.1682 


Table 6-3: Comparison of the normalized stress intensity factors obtained from 
a yielded clad and an elastic clad for various values of the crack 

length, t 0 =0.0, a 0 7 '=-a’ 1 £ 1 0 O /1 -Vj, for Model I.(Material pair A) 



k(b 2 )/* 0 T VTj 2 as r - 

► co (yield clad) 

k(b 2 ) / C 0 T V1J2 as r - 00 
(elastic clad) 

i 2 /h 

°fI a o r=1 

1/2 

1/5 

1 / 10 

0.025 

0.92157 

5.39959 

8.08640 

8 . 98201 

0.7654 

1.0 

1 .06539 

1 . 48309 

1 . 73372 

1.81726 

0 . 7740 

4.0 

1 . 14489 

1 . 27478 

1.35272 

1.37870 

0.8285 

6.5 

1 . 17684 

1 . 26092 

1.31138 

1 . 32820 

0.8532 
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Table 6-4: Normalized stress intensity factors for an edge crack in a 
homogeneous strip of thickness L under uniform load c 0 (without 
elastic foundation, x=0.0 and with elastic foundation R i /L= 9.0 and 
Xh/£ = 0.01108) 



Kay a 

present work 

bJL 

*( 6 l) / v/ *l ’ ^ =0 

*( 6 i) / > ^ =0 

k(b l )/c 0 Vf l ,R i /L = 9.0 

0.001 

1.1216 

1.1216 

1.1215 

0.1 

1 . 1892 

1.1892 

1.1574 

0.2 

1.3673 

1.3673 

1 . 2726 

0.3 

1.6699 

1.6699 

1.4267 

0.4 

2.1114 

2.1114 

1.6983 

0.6 

2.8246 

2.8246 

1.8135 

0.6 

4 . 0332 

4.0332 

2.0775 

0.7 

6 . 3649 

6.3653 

2.3801 

0.8 

11.966 

11.958 

2.7399 

0.86 

18.628 

18.636 

2.9835 

0.9 

34.632 

34.641 

s _ 

3.3660 


Table 6-5: Normalized stress intensity factors for a crack terminating at the 
interface under uniform load o 0 , (/ 2 /^i = 0.1 in the present work). 




*( a 2 ) / *o v/ y^ 

k(b 2 )/a 0 Vy~2 

Aluminum-Epoxy 

Kay a 

2.7997 

0.8826 

1/^=0 . 3 , l/^—O . 35 

^l/ /i 2 ==23 ‘° 77 

a 2 =0.3381 

present work 

2.7948 

O . 8842 

Epoxy -Aluminum 

Kay a 

O . 0955 

1.3398 

i/^O.35 , l/ 2 = 0.3 

^ 1 2 = ° * 0433 
q 2 =0.8248 

present work 

0.0983 

1 . 3403 
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Table 6-6: Normalized stress intensity factors for an embedded crack going 

through the interface under loading P 1 /P 2 = (£ 1 /£ 2 )( 1 ~ v 2 / 1 ~ v i)» 
v^O.3, v 2 = 0.35, Mi / 1*2 = 23 .077, p = Pj = 02 = 0.273692, / 1 /* 1 = 0.1, 

12/*!= 0.05 ,0.025. 


l 2 / *1 

i 

M“j) 

k 

X 

k 

*v 

Hb 2 ) 


Pj V / //2 

Pj (1/2)* 

Pj (»/2)* 

P 2 ^~2 

0.6 

present work 
Erd. & Biri. 

X . 2406 
1 . 2378 

- 0.06127 

- 0.06900 

- 0.02034 

- 0.02158 

1 . 5528 
1.6661 

0.25 

| 

present work 
Erd. & Biri. 

1 . 3409 
1.3324 

- 0.08456 

- 0.08113 

- 0.03088 

- 0.03184 

2.1320 

2.1391 


Table 6-7: Normalized stress intensity factors for an edge crack in a 
homogeneous strip of thickness L under uniform stress c 0 f° r 
fixed crack length f> 1 /£=0.5, various values of % (stiffness of the 
elastic foundation) 


X L , A’ / m 2 

k ( b l ) / ^0 V *1 

0.0 

2.8246 

io 3 

2.8238 

io 6 

2.5282 

IO® 

1 . 8940 

IO 12 

1 . 2000 

10 16 

1 . 1693 

10 20 

1.1601 

00 

1 . 1546 
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Table 6-8: Comparison between normalized stress intensity factors for an edge 
crack under unif orm loa din g a 0 in a homogeneous strip on an elastic 
foundation and a homogeneous cylinder R t /L = 9.0. 


bjL 

k ( b i) / a 0 ^ 

k ( b l)/ <T o'^ b l 


cylinder 

strip 

0.1 

1.158 

1.1574 

0.2 

1 .253 

1.2725 

0.3 

1 .392 

1.4267 

0.4 

1.568 

1 .5983 

0.5 

1 .779 

1 .8135 

0.6 

2.025 

2.0775 
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Table 6-9: Comparison between transient temperature distribution 0/ © 0 for a 
hollow cylinder and a homogeneous strip on an elastic foundation % 
due to a unit step temperature change on the inner wall, R i /L=9.0 

and xL/E=0.0im, x=tD / L 2 . 




table 6-9, continued 


x / L 

0 /©O , T=0.1 

© /© 

o 

II 

o 


Nied 

present work 

i Nied 

present 'work 

0.00 

1.000 

i 

1 . 000 I 

0.999 

1.000 

0.05 

0.908 

0.911 

0.967 

0.971 

0.10 

0.818 

0.823 

0.937 

0.942 

0.15 

0.731 

0.737 

0.906 

0.913 

0.20 

0.648 

0.665 j 

0.877 

0.885 

0.25 

0.568 

0.576 i 

0.848 

0.858 

0.30 

0.494 

O . 602 ! 

0.820 

0.832 

0.35 

0.426 

0.434 

0.794 

0.806 

0.40 

0.363 

0,371 

0.769 I 

0.782 

0.45 

0.307 

0.315 

j 0.746 

0.759 

0.50 

0.257 

0.264 

j 0.723 

0.738 

0.55 

0.213 

0.220 

0.703 

0.718 

0.60 

O. 176 

0.181 | 

0.685 

0.700 

0.65 

O. 144 

0.149 

0.668 

0.684 

0.70 

0.117 

0.121 

O . 654 

0.670 

0.75 

0.095 

0.099 

0.641 

0.657 

0.80 

0.078 

0.081 

0.631 

O . 647 

0.85 

0.065 

0.067 

0.623 

0.639 

0.90 

0.056 

0.058 

0.617 

0.634 

0.95 

0.050 

0.053 

0.614 

0.630 

1 .OO 

0.048 

0.061 

0.613 

0.629 
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Table 6-10: Comparison between Stress distribution o^/ a 0 for a hollow 
cylinder and a homogeneous strip on an elastic foundation % due to 
a unit step temperature change on the inner wall, R i /L= 9.0 and 

XL/£=0.0ll08,t = r£>/L 2 . 


x / L 

a 

/ a n T 9 T = 0-01 

yy ' 0 7 

r 

a / 
yy ! o 

, t = 0.05 


Nxed 

present work 

Nied 

present work 

0.00 

0.892 

0.8872 

0.768 

0.7477 

0.06 

0.614 

0.6108 

0.630 

0 . 6221 

0.10 

0.369 

0.3666 

0.606 

0.4995 

0.15 

O . 179 

0.1760 

0.388 

0.3829 

0.20 

0.048 

O . 0446 

0.280 

O . 2748 

0.25 

- 0.031 

- 0.0357 

0.182 

0.1769 

0.30 

- 0.074 

- 0.0789 

0.096 

O . 0905 

0.35 

- 0.094 

- 0.0996 

0.022 

0.0161 

0.40 

- O . 103 

- 0.1082 

- 0.040 

- 0.0464 

0.46 

- 0.106 

- 0.1114 

- 0.091 

- 0.0976 

0.50 

- O . 107 

- 0.1124 

- 0.131 

- 0.1386 

0.56 

- 0.107 

- 0.1127 

- 0.162 

- 0.1703 

0.60 

- 0.107 

- 0.1128 

- 0.186 

- 0.1945 

0.66 

- 0.107 

- 0.1128 

- 0.203 

- 0.2125 

0.70 

- 0.107 

- 0.1128 

- 0.216 

- 0.2264 

0.76 

- 0.107 

- 0.1128 

- 0.224 

- 0.2346 

0.80 

- 0.107 

- 0.1128 

- 0.230 

- 0.2308 

0.86 

- 0.107 

- 0.1128 

- 0.234 

- 0.2449 

0.90 

- 0.107 

- 0.1128 

- 0.237 

- 0.2474 

0.95 

- 0.107 

- 0.1128 

- 0.238 

- 0.2488 

1.00 

- 0.107 

- 0.1128 

- 0.239 

- 0.2492 
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table 6-10, continued 


1 

x / L 

I » r=01 


ff / <rJ 
yy 1 o 

, t — 0.5 


Nied 

present, work 

Nied 

present work 

0.00 


0.666 

0.6432 


0.251 


0.2361 

0.05 


0.565 

0.6642 


0.220 


0.2069 

0.10 


0.475 

0.4663 


0.189 


0.1780 

0.15 


0.388 

0.3805 


0.169 


0.1495 

0.20 


0.304 

0.2979 


O. 129 


0.1215 

0.25 


0.226 

O . 2184 


0. lOO 


0.0942 

0.30 


0.151 

0.1467 


0.073 


0.0677 

0.36 


0.083 

0.0772 


0.046 


0 . 0423 

0.40 


0.020 

0.0146 


0.021 


0.0181 

0.45 

- 

0.036 ! 

- 0.0420 

- 

0.002 

- 

O . 0047 

0.50 

- 

0.086 

- 0.0926 

- 

0.025 

- 

0.0261 

0.55 

- 

O. 130 

- O. 1369 

- 

0 . 046 

- 

O . 0459 

0.80 

- 

0.167 

- 0.1754 

- 

0.063 

- 

0.0639 

0.66 

- 

0.200 

- 0.2082 

- 

O . 080 

- 

0.0801 

0.70 

- 

0.226 

- 0.2356 

- 

0.094 

- 

O . 0943 

0.75 

- 

0.248 

- 0.2581 

- 

O. 106 

- 

O. 1065 

0.80 

- 

0.265 

- 0.2759 

- 

O. 117 

- 

0.1186 

0.85 

- 

0.279 

- 0.2894 

- 

O. 125 

- 

O. 1246 

0.90 

- 

0.288 

- 0.2988 

- 

O. 130 

- 

0.1302 

0.95 

- 

0.293 

- 0.3043 

- 

O. 134 

- 

O. 1336 

1 .OO 

i 

0.295 

- 0.3061 


0.135 


0.1347 
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Table 6-11: Comparison between normalized stress intensity factors for an edge 
crack subjected to transient thermal stresses in a hollow cylinder 
and a homogeneous strip on an elastic foundation x> Rj/L=9.0 and 

XL/£=0.01108, x-tD/L 2 , o o r =-a’£0 o /(l-v). 


bJL 

*(*i)/* 0 r 

y / T x , r = 0.01 

HbJ/'o 1 

, t = 0.05 


Nied 

present work 

Nied 

present work 

0.01 

0.962 

0.9666 

0.833 

0.8214 

0.1 

0.667 

0.6532 

0.701 

0.6913 

0.2 

0.426 

0.4316 

0.589 

0.6916 

0.3 

0.300 

0.3072 

0.601 

O . 5094 

0.4 

0.238 

0.2402 

0.432 

0.4363 

0.6 

0.206 

0.2066 

0.378 

0.3793 

0.6 

0.186 

O . 1868 

0.334 

0.3369 


bJL 

*(M/*o 7 

>i 

II 

o 

o 

II 

p 


Nied 

present work 

Nied 

present work 

0.01 

0.724 

0.7093 

0.277 

O . 2608 

O . 1 

0.633 

0.6207 

0.247 

0.2327 

0.2 

0.660 

0.6589 

O . 224 

O . 2147 

0.3 

0.502 

0.5069 

0.206 

O . 1997 

0.4 

0.453 

0.4636 

0.190 

0.1829 

0.5 

0.408 

O . 4085 

0.174 

0.1677 

0.6 

0.367 

0.3687 

0.168 

0.1629 















Table 6-12: The maximum of normalized stress intensity factors for an under- 
clad crack subjected to transient thermal stresses, = 3.0, 

R i /L= 9.0 and %L/£ 2 = 0.01 108 . (material pair A). 



Nied 

present work 

l 2 /*! 

T - 

= 0.0 

T = 

0.0 

T = 1.0 


k(a 2 ) 

*(ft 2 ) 

a 2 ) 

k(b 2 ) 

*(a 2 ) 

k(b 2 ) 

V^V 2 

T / 

V V 2 

°0 ‘2' ' 

'o^V 2 

'o Tn/ V 2 

°J^ l 2' 2 

0.004 

0.0114 

0.0111 

0 . 0068 

0 . 0065 

0.0061 

O . 0048 

0.040 

0.0101 

0.0073 

0.0066 

O . 0027 



0.200 

O . 0046 

- 0.0088 

0 . 0000 

- 0.0138 



0.400 

- 0.0022 

- 0.0277 

- 0.0071 

- 0.0333 



0.600 

- 0.0090 

- 0.0453 

- 0.0142 

- 0.0614 






Table 6-13: The maximum of normalized stress intensity factors for an under- 
clad crack subjected to transient thermal stresses, h 2 /h l =9.0, 
R i /L = 9.0 and xL/£ 2 =0.01108 . (material pair A). 



Nied 


present 

•work 



o 

d 

II 

o 

d 

II 

T= 10.0 


*(a 2 ) 

*(* 2 ) 

k(a 2 ) 

*(» 2 > 

*(“ 2 ) 

*(* 2 ) 


a T ^~2 
a 0 l 2' Z 

fj t V i it 
"0 l 2' " 

% Ty/ h^ 

^0 Tv V2 


0,01 

0.2245 

1 

0.2240 

0,2270 

0.2264 

0.2032 

O . 2027 

0.5 

0.2150 

0.1897 

0.2171 

0.1915 

0.1945 

0.1717 

1 .O 

O . 2080 

0.1589 

0.2101 

0.1603 

0.1883 

0.1437 

2.0 

O. 1951 

0.1030 

0.1950 

0.1018 

0.1749 

O . 0906 

3.0 

0.1800 

0.0523 

0.1790 

0 . 0496 

0.1604 

0.0429 

4.0 

O. 1618 

0.0066 

0.1604 

i 

0 . 0026 

0.1434 

- 0.0002 
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Table 6-14: The maximum of normalized stress intensity factors for an under- 
clad crack subjected to transient thermal stresses, /j 2 //i 1 = 24.0, 
R i /L = 9.0 and %L/£ 2 = 0.01 108 . (material pair A). 



Nied 


present, work 


*2 / *1 

T - 

= 0.0 

o 

o 

II 

t = 20.0 


k(a 2 ) 

*(» 2 ) 

k(a 2 ) 

*(» 2 ) 

*(« 2 ) 

*(» 2 ) 

'o rv y2 


% TX ~ 2 

'o TV V* 

% Ty/ V 2 

'o T ^ 

0.026 

0.3908 

0.3899 

0.3980 

0.3971 

0.3877 

0.3869 

0.6 

0.3873 

0.3701 

0.3941 

0.3766 

0.3842 

0.3677 

1 .6 

0 . 3904 

0.3369 

0.3980 

0.3436 

0.3886 

0.3366 

4.0 

0 . 4006 

0 . 2622 

0.4067 

0.2657 

0.3974 

0.2612 

9.0 

0.3919 

0.1309 

0.3942 

0.1306 

0.3874 

0.1285 

14.0 

0.3494 

0 . 0200 

0.3520 

0.0166 

0.3463 

0.0161 
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Table 6-15: The maximum of normalized stress intensity factors for an edge 
crack subjected to transient thermal stresses, /ij = 3.0, RJL= 9.0 
and = 8 . (material pair A). 



Nied 

present work 

l/\ 

o 

o 

II 

o 

o 

II 

q 

d 

II 

k. 


*(» 2 ) 

*(* 2 > 

k(b 2 ) 





1 .004 

0.4760 

0.4629 

0.3793 

1.2 

0.4097 

O . 4000 

0.3143 

1 .4 

0.3726 

0.3574 

0.2779 

1.6 

O . 3449 

0.3252 

0.2615 

2.0 

O . 3062 

0.2759 

0.2127 
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Table 6 - 16 : The maximum of normalized stress intensity factors for an edge 
crack subjected to transient thermal stresses, h 2 /h 1 =9.0, R i /L=9.0 
and %L/£ 2 = 0.01108 . (material pair A). 



Nied 

present work 

l/h 1 

q 

d 

II 

■4 

II 

q 

d 

T= 10.0 


*(* 2 ) 

k\b 2 ) 

*(* 2 ) 



T /- 

'o v ' 

1 .Ol 

O . 6090 

0.6298 

0.6898 

1 .5 

0.4740 

0 . 6050 

0.4683 

2.0 

0.4151 

0.4529 

0.4196 

3.0 

O . 3461 

0 . 4044 

0.3751 

4.0 

0.3034 

0.3867 

0.3589 

6.0 

0.2722 

0.3839 

0.3567 
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Table 6-17: The maximum of no rmaliz ed stress intensity factors for an edge 
crack subjected to transient thermal stresses, A 2 M 1 = 24.0, R t /L = 9.0 
and xL/£ 2 s 0.01108 . (material pair A). 



Nied 

present, work 

l/h, 

o 

© 

II 

T — 0.0 

r= 20.0 


*(* 2 ) 

k(b 2 ) 

*(» 2 ) 




u r * 

1 .025 

0.7383 

0.7531 

0.7372 

1 .5 

0.6136 

0 . 6307 

0.6168 

2.0 

0.5554 

0.5708 

0.5689 

5.0 

0.4071 

0.4353 

0.4282 

7.6 

0.3531 

0.4172 

0.4110 
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•40 .60 .80 1.00 



Figure 6-2: The normalized transient stress distribution in Model I for 

t q =0.0 ,x 0 =t 0 D 1 /h 2 , c 0 T =- a\ E x e 0 /(l -v^. (Material pair A) 
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I i I I 

.0 8.0 16.0 24.0 32.0 


x / h 

Figure 6-3: The normalized transient temperature distribution in Model I for 
x o =10.0 , i 0 =t 0 D l /h 2 , (Material pair A) 
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Figure 6-4: The normalized transient stress distribution a^/ c 0 r in Model I for 
x 0 = 10.0 , x 0 =t 0 DJh 2 , o 0 r =- a \ E x 0 O /(1 -Vj) (Material pair A) 
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I i i — 1 1 

.0 8.0 16.0 24.0 32.0 

x / h 

Figure 6-5: The normalized tr ans ient temperature distribution in Model I for 
t o =20.0 ,x 0 =t 0 D l /h 2 y (Material pair A) 
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Figure 6-6: The normalized transient stress distribution g^/g^ in Model I for 
t o =20.0 , z 0 =t 0 D l /h 2 , g 0 t =- ct\ E^QqK 1-v^ (Material pair A) 
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Figure 6-8: The normalized transient stress distribution o^/Oq 7 in Model I for 
t o =40.0 , x 0 =t 0 D l /h 2 , o 0 r =- ^ EjOo/a-Vj) (Material pair A) 
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Figure 6-9: The normalized transient temperature distribution in Model I for 
x o =60.0 , XQ=t Q D l /h : 2 , (Material pair A) 
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for t 0 =60.0 , -r 0 =r 0 D 1 //i 2 , c 0 T =- cl\ E 1 0 o /( 1 -v^ (Material pair A) 
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Figure 6-11: The normalized stress intensity factor k(b^) as a function of 
nondimensional time t for an edge crack of various lengths in 

Model I for t o =0.O , x 0 =t Q D y /h 2 , c 0 T =-a\ 0^/(1^) (Material 

pair A) 
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Figure 6-12: The influence of x 0 on the normalized stress intensity factor k(bj) 
as a function of nondimensional time x for an edge crack of length 

ljh= 0.5 in Model I , XQ=t 0 D 1 /fp- , a 0 T =-a\ £ 1 0 o /( l-Vj) (Material 
pair A) 
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Figure 6-13: The influence of x 0 on the normalized stress intensity factor 

as a function of nondimensional time t for a broken clad in Model 

1 , 1 0 =t 0 D l /h 2 , c 0 T ~- cx^ E^q/CI -Vj) (Material pair A) 
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Figure 6-14: The normalized stress intensity factor k (a^) as a function of 
nondimensional time t for an under-clad crack of various lengths 

in Model I for T o =0.0 , x^t^D^h 1 , o 0 r =- a\ E^Kl-v^ 
.(Material pair A) 
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Figure 6-15: The normalized stress intensity factor k{b 2 ) as a function of 
nondimensional time x for an under-clad crack of various lengths 

in Model I for t 0 =0.0 , T 0 =r 0 D 1 /A 2 , oJ=- a\ E^q/O- -Vj) 
.(Material pair A) 
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Figure 6-16: The influence of x 0 on the normalized stress intensity factor k(a 2 ) 
as a function of nondimensional time x for an under-clad crack of 

length l 2 /h= 1.0 in Model I , x 0 -tQD^h 2, , c 0 r =- 0^ ^©q/CI-Vj) 
(Material pair A) 
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Figure 6-17: The influence of x 0 on the normalized stress intensity factor k(b 2 ) 
as a function of nondimensional time x for an under-clad crack of 

length l 2 /h= 1.0 in Model I , x 0 =t 0 D l /h 2 , c 0 T =- ct’ l E l Q 0 /(l-\ 1 ) 
(Material pair A) 
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Figure 6-18: The normalized stress intensity factor k{b^) as a function of 
nondimensional time x for an edge crack of various lengths 
crossing the interface of in Model I for x o =0.0 , x 0 =t 0 D l /h 2 , 

c 0 T =- ccj E y 0 o /( 1 -Vj) (Material pair A) 
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Figure 6-19: The influence of x 0 on the normalized stress intensity factor k(b 2 ) 
as a function of nondimensional time % for an edge crack crossing 
the interface in Model I for l/h= 4.0 , Zq=IqD 1 / A 2 , 

o 0 r =- a’j £ 1 0 O /(1 — Vj) (Material pair A) 
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Figure 6-20: The normalized stress intensity factor £(u ^ as a function of 
nondimensional time z for an embedded crack of various lengths 
crossing the interface in Model I for b 2 /h =2.0 , t o = 0.0 , x 0 =r 0 D 1 //i 2 

, o 0 T =- a’j £ 1 © 0 /(l -Vj) (Material pair A) 
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Figure 6-21: The normalized stress intensity factor k{b 2 ) as a function of 
nondimensional time x for am embedded crack of various lengths 
crossing the interface in Model I for b 2 /h=2.0 , x Q =0.0 , x 0 =t Q D l /h 2 

, c 0 T =- a’ L E x © 0 /( 1 --Vj) (Material pair A) 
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Figure 6*22: The influence of x$ on the normalized stress intensity factor k (a^) 
as a function of nondimensional time x for an embedded crack 
crossing the interface, a l /h= 0.2 , bjh= 2.0 in Model I , XQ=t Q D l /h 2 

, c 0 r =- a\ E x 0 O /(1 -Vj) (Material pair A) 
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Figure 6-23: The influence of x 0 on the normalized stress intensity factor k(b 2 ) 
as a function of nondimensional time x for an embedded crack 
crossing the interface in Model I for a l /h=0.2 , b^/h= 2.0 , 

l 0 =t 0 D 1 /h 2 , c 0 T =- aj EjGq/CI-Vj) (Material pair A) 
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Figure 6-24: The normalized transient temperature distribution in Model I for 
x o =0.0 , x 0 =t 0 D l /h 2 , (Material pair B) 
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Figure 6-26: The normalized transient temperature distribution in Model I for 
t 0 =20.0 , x 0 =t Q D l /h 2 , (Material pair B) 
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for x o =20.0 , x 0 =t 0 D l /h 2 , o 0 T =- a\ E i © 0 /(l -v L ) (Material pair B) 
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for t o =60.0 , XQ=t Q DJk L , c 0 r =- ct\ E^q/ (1-Vj) (Material pair B) 


165 




.0 8.0 16.0 24.0 32.0 


r^tDjh 2 

Figure 6-30: The normalized stress intensity factor k(b^) as a function of 
nondimensional time x for an edge crack of various lengths in 

Model I for x o =0.0 , x 0 =t 0 D l /h l , a 0 r =-a , 1 E X Q 0 /(l-v{) (Material 
pair B) 
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/j/A=0.5 in Model I , t 0 =f 0 D 1 /A 2 , a 0 r =- a\ E^/Cl-v^ (Material 
pair B) 
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Figure 6 - 32 : The influence of T 0 on the normalized stress intensity factor kQ}^) 
as a function of nondimensional time t for a broken clad in Model 

I , x 0 =/ 0 D 1 //i I 2 , o 0 r = - a’ x © 0 / (1 - v i) , (Ij =0.552538 . (Material 

pair B) 
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Figure 6*33: The normalized stress intensity factor k x as a function of 
nondimensional time T for an edge crack of various lengths 
crossing the interface in Model I for Tq = 0.0 , t 0 =/ 0 £> 1 //i 2 , 

a 0 r =- aj Fj© 0 /(l — v A ) , Pi = == 0.01872238 . (Material pair B) 
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Figure 6-34: The normalized stress intensity factor as a function of 
nondimensional time t for an edge crack for various lengths 
crossing the interface in Model I for x o =0.0 , x 0 =r 0 D 1 //i 2 , 

a 0 T =- a\ E y Q q/CI-Vj) , P 1 = a 2 = 0.01872238 . (Material pair B) 
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Figure 6-35: The normalized stress intensity factor k(b 2 ) as a function of 
nondimensional time X for an edge crack of various lengths 
crossing the interface in Model I for x o =0.0 , z 0 =t 0 D 1 /h 2 , 

o 0 r =- cq £ x 0 0 Kl-v{) , Pi = 0(2 = 0.01872238 . (Material pair B) 
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Figure 6-36: The influence of t 0 on the normalized stress intensity factor k x as 
a function of nondimensional time x for an edge crack crossing the 

interface l/h= 2.0 in Model I , x Q =tQD l /h 2 > a ’i » 

Pi = a 2 = 0.01 872238 . (Material pair B) 
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Figure 6-37: The influence of t 0 on the normalized stress intensity factor as 
a function of nondimensional time T for an edge crack crossing the 

interface l/h=2.0 in Model I , % 0 =t 0 D i /h 2 , cj =- E^q/CI-v^ , 

(3 =02 = 0.01872238 . (Material pair B) 
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Figure 6-38: The influence of x 0 on the normalized stress intensity factor k(b 2 ) 
as a function of nondimensional time x for an edge crack crossing 
the interface ///t=2.0 in Model I , x 0 =r 0 D 1 //z 2 , 

o 0 r =- a’j E l O 0 /(l-v l ) , P 1 = cx 2 =0.01872238 . (Material pair B) 
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Figure 6-40: The normalized transient stress distribution o^/Oq^ in Model II 

for t o =0.0 , h 2 lh x ='l . 0 , t 0 =f 0 zyAi , o 0 r — ai^Go/a-Vj). 
(Material pair A) 
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Figure 6-41: The normalized transient temperature distribution in Model II for 
Tq= 1.0 , =3.0 , ZQ=t 0 D l /h^ , (Material pair A) 


'8 




.30 .00 .30 .60 .90 


1 



Figure 6-42: The normalized transient stress distribution in Model II 

for < Cq=1.0 , h^lh^ — 'i.O , ttj 0 q/ ( 1 ~ Vj). 

(Material pair A) 
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Figure 6-42, continued 
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Figure 6-44: The normalized transient stress distribution in Model II 

for Tq— 6.0 , h^/ hy=3.0 , Xq=?q.Dj//Ij , cx^ £j©q/( 1— Vj). 

(Material pair A) 
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Figure 6-44, continued 
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x o =0.0 , h 2 lh l = 9.0 , x Q =t Q D 1 /h\ , (Material pair A) 
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Figure 6 - 46 : The normalized transient stress distribution Gy y /GQ 1 ' in Model II 

for t o =0.0 , tyAj-9. 0 , x 0 =t Q DJh\ , o 0 r =- ct\ E^qKI-vJ. 
(Material pair A) 
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for t o =6.0 , =9.0 , 't 0 =t 0 D l /h 2 l , a 0 T =-a\E l Q 0 /(l-v l ). 

(Material pair A) 
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Figure 6-49: The normalized transient temperature distribution in Model II for 
t o =10.0 , h 1 /h l =9.0 , t 0 =t 0 D l /h 2 l , (Material pair A) 
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Figure 6-50: The normalized transient stress distribution C yy /CQ T in Model II 

for t o =10.0 , A 2 /A 1 = 9.0 , x Q =t 0 D 1 /h 2 l , C ( ?=‘-(x\E l Q 0 /(l-v l ). 
(Material pair A) 






Figure 6-52: The normalized transient stress distribution a^/c 0 r in Model II 

for t 0 =20.0 , /* 2 //* 1 = 9.0 , , O 0 T =-(X\ E^Kl-w^). 

(Material pair A) 






Figure 6-54: The normalized transient stress distribution a^/c 0 r in Model II 

for t 0 =0.0 , h 2 th x = 24.0 , x 0 = r 0 D l/ /l ? > C 0 T =- a\ E^q/0-vJ. 
(Material pair A) 
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Figure 6-56: The normalized transient stress distribution o^/Oq 7 in Model II 

for t 0 =6.0 , tyAj- 24.0 , , o 0 r =- ai E^Kl-v^ 

(Material pair A) 
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Figure 6-57: The normalized transient temperature distribution in Model II for 
t 0 =10.0 , A 2 //i 1 = 24.0 , ^Q-tQDJn^ > (Material pair A) 
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Figure 6-58: The normalized transient stress distribution c^/Gq^ in Model II 

for ^=10.0 , h 2 /h l = 24.0 , z 0 =t 0 D x lh\ , o 0 r =- a\ 

(Material pair A) 
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Figure 6-61: The normalized transient temperature distribution in Model II for 
t o =40.0 , h 2 /h l = 24.0 , x 0 =t Q D l / , (Material pair A) 
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Figure 6-62: The normalized transient stress distribution cr^/ o 0 T in Model II 

for t 0 =40.0 , V A l = 24.0 , T 0 =f 0 £> 1 /^ , a 0 T =-a\E l Q Q Kl-\ l ). 
(Material pair A) 
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Figure 6-63: The normalized stress intensity factor k{b^) as a function of 
nondimensional time i for an edge crack in Model II for t 0 =0.0 , 

h 2 /h l = 3.0 , /?,/L=9.0 and xL/E 2 = 0.01108 * ^o =t o D if h \ > 

a 0 T =- aj E x 0 o /( 1 — Vj) (material pair A) 
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Figure 6-64: The normalized stress intensity factor k(b^) as a function of 
nondimensional time i for an edge crack in Model II for Xq=6.0 , 

h 2 /h l = 3.0 , R i IL= 9.0 and x^/ £ 2 = 001108 » t o='o z V /l i > 

a 0 T =-a\ £ 1 © 0 /(l -Vj) (material pair A) 


*(»,) / % T A 


I fiifl 

I 


r o = *0^1 /** 



12.0 


16.0 


Figure 6-65: The influence of t 0 on the normalized stress intensity factor k(b^) 
as a function of nondimensional time x for an edge crack of length 
/ 1 /A 1 =0.2 in Model II , 0 , /?,/L= 9.0 and xL/£ 2 =0.01108 » 

T 0 =f 0 £) l/ /l i » a o r =- a ;£i0 o /(l- v i) (material pair A) 
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Figure 6-66: The influence of x 0 on the normalized stress intensity factor k(b x ) 
as a function of nondimensional time x for an edge crack of length 
/j/A^O.9 in Model II , h 2 /h l =3.0 , R i /L=9.0 and ‘XJL/E 2 -O.OHOS , 

x 0 =r 0 f>i/Ai , a 0 T =- a’, Ey 0 o /( 1 -Vj) (material pair A) 
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Figure 6-67: The normalized stress intensity factor k(b}) as a function of 
nondimensional time i for an edge crack in Model II for r 0 = 0.0 , 

h 2 /h l =9. 0 , R i /L-9.0 and xL/£ 2 =0.01108 , x^^D^h] , 

a 0 r =- a’j E 1 © 0 /(l -Vj) (material pair A) 
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Figure 6-68: The normalized stress intensity factor k{b^) as a function of 
nondimensional time x for an edge crack in Model II for x 0 = 10.0 , 

h 2 /h l = 9.0 , R i /L=9.0 and xZ./£ 2 =0.0H08 > t 0 = r 0 Z V /, l > 

<j 0 r =- a’, E l 0 O /(1 -v L ) (material pair A) 
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Figure 6-69: The influence of t 0 on the normalized stress intensity factor k(bj) 
as a function of nondimensional time t for an edge crack of length 
0.2 in Model II , f^/h^.O , RJL= 9.0 and xL/£ 2 =0.01108 , 

x 0 = t 0 D l /h 2 i , o 0 T =-a’ 1 E i e 0 /(\-v ] ) (material pair A) 
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Figure 6-70: The influence of x 0 on the normalized stress intensity factor k(by) 
as a function of nondimensional time t for an edge crack of length 
/j/Aj- 0.9 in Model II , f^/h^.O , Ri/L=9.0 and xL/E 2 =0.0im , 

T Q =t 0 D l lh 2 [ , cj=- a’, E l 0 O /(1 -v,) (material pair A) 
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Figure 6-71: The normalized stress intensity factor k(b{) as a function of 
nondimensional time t for an edge crack in Model II for t o =0.0 , 

A 2 /A 1 = 24.0 , R^L =9.0 and x^/^2 = 0 - 01108 » x 0 =t 0 D J h \ * 

c 0 T =- cl \ E l © 0 /(l — Vj) (material pair A) 
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Figure 6 - 72 : The normalized stress intensity factor kibj) as a function of 
nondimensional time x for an edge crack in Model II for x o =20.0 , 

/j 2 /A 1 =24.0 , R i /L = 9.0 and xL/£ 2 = 0.01 108 , x 0 =t Q D l /h\ , 

c 0 r =- a’j E j 0 q/(1 -Vj) (material pair A) 


220 


r 0 *0 



E 2' V 2 


h 2 f h x = 24.0 
lJh^O.2 


15.0 

tDJhl 


22.5 


30.0X10 


Figure 6-73: The influence of x 0 on the normalized stress intensity factor k(b^) 
as a function of nondimensional time x for an edge crack of length 
IJh^ 0.2 in Model II , /i 2 /A 1 = 24.0 , R { /L= 9.0 and xL/E 2 = 0.0110S , 

x 0 =t 0 D l /h 2 l , o 0 t =-ol\ E l Oq/CI-Vj) (material pair A) 
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Figure 6-74: The influence of x 0 on the normalized stress intensity factor k(by) 
as a function of nondimensional time x for an edge crack of length 
/j/Aj- 0.9 in Model II , A 2 /A 1 =24.0 , £,/L = 9.0 and %L/£ 2 = 0.01108 , 

- i 0 =t 0 D 1 /h \ , a 0 T =- a’, E x O 0 /(l -Vj) (material pair A) 
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Figure 6-75: The influence of x 0 on the normalized stress intensity factor k(p{) 
as a function of nondimensional time x for a broken clad in Model 

II hz/h^.O , R i /L=9.0 and xL/£ 2 =0.01108 , XQ=t Q D x ll? x , 
c 0 r =- a’j E x 0 O /(1 -v 1 ) (material pair A) 
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Figure 6-76: The influence of x 0 on the normalized stress intensity factor k(b x ) 
as a function of nondimensional time x for a broken clad in Model 

II /i 2 /A 1 =9.0 , /?,/L= 9.0 and xL/£ 2 =0.01108 * <t o='o i V /, i » 
o 0 T =- E l © 0 /(l — v : ) (material pair A) 
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Figure 6-77: The influence of Tq on the normalized stress intensity factor k(b^} 
as a function of nondimensional time i for a broken clad in Model 

II A 2 /A 1 = 24. 0 , Ri/L= 9.0 and xL/£ 2 =0.01108 « x 0 =t 0 D \ ll h * 
oJ=- o.\ £ 1 0 O /(1 -Vj) (material pair A) 
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Figure 6-78: The influence of x 0 on the normalized stress intensity factor k(a 2 ) 
as a function of nondimensional time i for an under-clad crack of 
length / 2 //* 1 = 0.004 in Model II V A l = 30 , RJL = 9.0 and 

%L/E 2 =0.01 108 , x Q =t 0 D l /h 2 l , c 0 T =- cc’j E l Q 0 /(l-v l ) (material pair 
A) 
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Figure 6-79: The influence of x 0 on the normalized stress intensity factor k(b 2 ) 
as a function of nondimensional time t for a fixed under-clad 
crack length / 2 Mi = 0.004 in Model II , RJL = 9.0 and 

XL/£ 2 =0.01108 , t 0 =f 0 Z) 1 //ii , aj=- a\ E x (material pair 
A) 
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Figure 6-80: The normalized stress intensity factor k (a 2 ) as a function of 
nondimensional time i for an under-clad crack in Model II for 

t 0 =0.0 , A 2 /A 1 =9.0 , R i /L=9.0 and xL/£ 2 =001108 » T 0 =? 0 Z) l/ /l l » 

Cq=- a’, 0 O /(1 -Vj) .(material pair A) 
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Figure 6-81: The normalized stress intensity factor k(b 2 ) as a function of 
nondimensional time x for an under-clad crack in Model II for 

x 0 = 0.0 , A 2 //i 1 = 9.0 , R i /L=9.0 and %L/E 2 = 0.01 108 , t 0 =f 0 D 1 /^ , 
o 0 r =- a’j E 1 0 o /( 1 -Vj) .(material pair A) 
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Figure 6*82: The normalized stress intensity factor k (a^) as a function of 
nondimensional time t for an under-clad crack in Model II for 

t 0 =10.0 , / t 2 / A 1 =9.0 , /?,-/£= 9.0 and x ^/£ 2 =0 - 01108 » <t o='o z V /l i » 

o 0 r =- a’j E l 0 O /(1 -Vj) .(material pair A) 
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Figure 6-83: The normalized stress intensity factor k(b^) as a function of 
nondimensional time x for an under-clad crack in Model II for 

x o =10.0 , A 2 //i 1 =9.0 , R i IL=9.0 and xL/E 2 =0.01108 , 't 0 ='o z V A i > 
o 0 r =- cc’j £j O 0 /(l -Vj) .(material pair A) 
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Figure 6-84: The influence of t 0 on the normalized stress intensity factor k(a 2 ) 
as a function of nondimensional time x for an under-clad crack of 
length / 2 //i 1 = 1.0 in Model II h 2 lh l = 9.0 , RJL= 9.0 and 

XL/£ 2 =0.01108 , x 0 =t 0 D 1 /h 2 l , o 0 r =- © 0 /(l -Vj) (material pair 

A) 
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Figure 6-85: The influence of Tq on the normalized stress intensity factor k (£> 2 ) 
as a function of nondimensional time % for an under-clad crack of 
length 1.0 in Model II h 2 /h 1 =9.0 , R i /L = 9.0 and 

XL/E 2 =0.0U08 , t 0 =/qD 1 /Aj , o 0 r =- a\ ^©q/CI-Vj) (material pair 
A) 
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Figure 6-86: The influence of x 0 on the normalized stress intensity factor k(a 2 ) 
as a function of nondimensional time x for an under -clad crack of 
length ^/h^S.O in Model II ^2/A 1 =9.0 , R i /L=9.0 and 

XL/£ 2 =0.01108 , 'i 0 =t 0 D 1 /h 2 l , a 0 T =- ct\ E x e (material pair 
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Figure 6-87: The influence of t 0 on the normalized stress intensity factor k(b 2 ) 
as a function of nondimensional time t for an under-clad crack of 
length f 2 /^l = 3- 0 i 11 Model II h 2 /h 1 = 9.0 , R i /L= 9.0 and 

%L/E 2 = 0.01108 , x 0 =t 0 D l /h 2 1 , a o r =-a’ 1 £ 1 0 o /(l-v 1 ) (material pair 
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Figure 6-88: The normalized stress intensity factor k (aj) , k(b 2 ) as a function of 
nondimensional time x for under-clad crack in Model II for 
l 2 /h x =9 . 0 , x o =0.0 , A 2 /A 1 =24.0 , R i /L=9.0 and %L/£ 2 = 0.01108 , 

x 0 =f 0 £> l/ ft i » °o r=_ a i £i© 0 /( 1-v l) -(material pair A) 
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Figure 6-89: The normalized stress intensity factor k(b^) as a function of 
nondimensional time t for an under-clad crack in Model II for 

t o =0.0 , /j 2 //i 1 =24.0 , RJL-9.0 and xE/E 2 =0.0\m > T 0 = 'o i V ,z l > 

c 0 T ~- a’. E l 0 O /(1 -Vj) .(material pair A) 
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Figure 6-90: The normalized stress intensity factor k(a 2 ) , k(b 2 ) as a function of 
nondimensional time x for an under-clad crack in Model II for 
/ 2 //j 1 =9.0 , x o =20.0 , A 2 /A 1 =24.0 , fl,/L=9.0 and xL/£ 2 =0.01108 , 

x 0 = t Q D i lh 2 l , g q t =- a’, E x 0 O /(1 .(material pair A) 
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Figure 6-91: The normalized stress intensity factor k(b 2 ) as a function of 
nondimensional time x for an under-clad crack in Model II for 

x o =20.0 , A 2 /A 1 =24.0 , R i /L=9.0 and x£/£ 2 =0 01108 » T 0 =t 0 D l/ h l > 

a 0 r =- a’j E x 0 o /( 1 -Vj) .(material pair A) 
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Figure 6-92: The influence of Iq on the normalized stress intensity factor k (tf 2 ) 
as a function of nondimensional time T for an under-clad crack of 
length 0.5 in Model II, A 2 /A 1 =24.0 , K.VL- 9.0 and 

XL/£ 2 = 0.01 108 , 'C 0 =f 0 D 1 /^ , o 0 T =- a\ (material pair 
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Figure 6-93: The influence of x 0 on the normalized stress intensity factor k(b 2 ) 
as a function of nondimensional time x for an under-clad crack of 
length / 2 //»i = 0.5 in Model II, h 2 /h 1 = 2A.O , R i /L= 9.0 and 

XL/E 2 = 0.01108 , x 0 =? 0 £> 1 //ij , o 0 T =- a’j EjQq/^-Vj) (material pair 
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Figure 6-94: The influence of Tq on the normalized stress intensity factor k ( a 2 ) 
as a function of nondimensional time T for an under-clad crack of 
length / 2 /A 1 =4.0 in Model II, h 2 /h l =24.0 , R i fL= 9.0 and 

XL/£,=0.01108 , x 0 =/qD 1 //|j , o 0 r =- a’ L EjQo/O-v^ (material pair 
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Figure 6-95: The influence of x 0 on the normalized stress intensity factor k (b 2 ) 
as a function of nondimensional time i for an under-clad crack of 
length / 2 //i 1 =4.0 in Model 13, /j 2 //i 1 =24.0 , R i /L= 9.0 and 

XL/E 2 =0.0U0& , x 0 =t 0 D 1 fh 2 l , Oo^-a^EiGo/Cl-Vi) (material pair 
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Figure 6-96: The normalized stress intensity factor k(b 2 ) as a function of 
nondimensional time x for an edge crack crossing the interface in 
Model n for x o =0.0 , 0 , /? t /L=9.0 and xL/E 2 =0.0im » 

x 0 =t 0 D l /h 2 1 , o 0 r =- a’j EjGo/fl -v^ (material pair A) 
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Figure 6-97: The normalized stress intensity factor k(b 2 ) as a function of 
nondimensional time x for an edge crack crossing the interface in 
Model II for x 0 =6.0 , /i 2 /A 1 =3. 0 , K,/L= 9.0 and xL/£ 2 =0.01108 , 

x Q =tQ D l /h^ , c 0 T =-a\ E { Q o/Cl-Vi) (material pair A) 
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Figure 6-98: The influence of T 0 on the normalized stress intensity factor k(b 2 ) 
as a function of nondimensional time t for an edge crack of length 
//Aj = 1.2 crossing the interface in Model II , k^/h^S . 0 , R t /L= 9.0 

and x^/£ 2 =0 - 01108 > <c o =r o £) i/ /l i • c o T= ~ °i ^i©o/( 1_v i) (material 
pair A) 
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Figure 6-99: The influence of x 0 on the normalized stress intensity factor k(b 2 ) 
as a function of nondimensional time x for an edge crack of length 
l/hy-2.0 crossing the interface in Model II , f^/hy^S . 0 , RJL— 9.0 

and xL/£ 2 =0.01108 , x 0 =r 0 £> 1 //ij , o 0 r =- a’jEjQo/Cl-Vj) (material 
pair A) 
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Figure 6-100: The normalized stress intensity factor k(b 2 ) as a function of 
nondimensional time x for an edge crack crossing the interface in 
Model II for x 0 =0.0 , h 2 /h x =9. 0 , tf f /L= 9.0 and x^/^ 0 - 01108 » 

x 0 =t 0 D 1 /h 2 l , o 0 r =- a’ x E x G q / (1 -Vj) (material pair A) 
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Figure 6-101: The normalized stress intensity factor k(b 2 ) as a function of 
nondimensional time x for an edge crack crossing the interface in 
Model II for x o =10.0 , h 2 lh l =9.0 , /? I /L=9.0 and xL/£ 2 =0.01108 » 

x 0 = t Q D l /h 2 l , c 0 T =-a\ E x 0 O /(1-Vj) (material pair A) 
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Figure 6-102: The influence of Xq on the normalized stress intensity factor k(b 2 ) 
as a function of nondimensional time X for an edge crack of 
length l/h^ = 1.5 crossing the interface in Model II , h 2 /h l = 9.0 , 

RJL= 9.0 and xL/E 2 =0.0U0S , X Q =t Q D l /h 2 l , c 0 T =-a\ E 1 Q 0 /(l-v l ) 
(material pair A) 
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Figure 6-103: The influence of Xq on the normalized stress intensity factor k(b 2 ) 
as a function of nondimensional time x for an edge crack of 
length l/hi=4.Q crossing the interface in Model II, A 2 /A 1 = 9.0 , 

R i !L= 9.0 and xL/£ 2 =0.01108 , x 0 =r 0 D 1 // i* , o 0 r =- ai ^©©/(l-v,) 
(material pair A) 
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Figure 6-104: The normalized stress intensity factor k(b 2 ) as a function of 
nondimensional time x for an edge crack crossing the interface in 
Model II for t 0 = 0.0 , /i 2 //z 1 =24.0 , K f ./L= 9.0 and %L/E 2 = 0.01108 , 

x Q =t Q D l lf ? l , o 0 T =- a’ x £' 1 0 O /(1 -v x ) (material pair A) 
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Figure 6-105: The normalized stress intensity factor k(b 2 ) as a function of 
nondimensional time t for an edge crack crossing the interface in 
Model II for x 0 =20.0 , * 2 /A 1 = 24.0 , RJL= 9.0 and xL/£ 2 = 0.01108 » 

x 0 =t 0 D l /h\ , o 0 r =- a\ E^qKI-v^ (material pair A) 


253 



W,)/% T St 


CU 



’=tDJh\ 

Figure 6-106: The influence of x 0 on the normalized stress intensity factor k(b 2 ) 
as a function of nondimensional time t for an edge crack of 
length l/h x = 1.5 crossing the interface in Model II , * 2 //i 1 = 24.0 , 

R J Lm 9.0 and xL/£ 2 = 0.01 108 , x 0 =/ 0 Z) 1 /a5 , c 0 T =-a\ 

(material pair A) 
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Figure 6-107: The influence of t 0 on the normalized stress intensity factor k(b^ 
as a function of nondimensional time t for an edge crack of 
length = 5.0 crossing the interface in Model II , ^ Mi =24.0 , 

R i /L = 9.0 and %L/£ 2 =0.01108 , , c 0 T =- cl\ 0 /(l-Vi) 

(material pair A) 
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Figure 6-108: The normalized transient temperature distribution in Model II 

2 

for t o =0.0 , h 2 /h l = 9.0 , t 0 =t 0 D 1 /h 1 , (material pair B) 
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Figure 6-109: The normalized transient stress distribution o^/Oq 7 in Model II 

for x 0 =0.0 , tyAj-9.0 , t 0 =/ 0 D 1 //i? , O 0 r =- Ct\ 

(material pair B) 
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Figure 6-111: The normalized transient stress distribution o^/c 0 r in Model II 

for t 0 =6.0 , h 2 /h l =9.0 , XQ=tQD^/h^ , Oq^"=— cx^ £j0q/(1— Vj). 
(material pair B) 


60 


2 < 



0 / 0 





Figure 6-113: The normalized transient stress distribution o^/o 0 ^ in Model II 

for x o =10.0 , h 2 fh l -9.0 , z 0 =t Q D l lh\ , c 0 T =- a\ E^Kl-v^ 
(material pair B) 








Figure 6-115: The normalized transient stress distribution o^/Oq^ in Model II 

for t o =20.0 , h 2 lh x =9 . 0 , , a 0 T =- a[ E^q/CI-v^. 

(material pair B) 
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Figure 6-116: The normalized stress intensity factor k(b j) as a function of 
nondimensional time x for an edge crack in Model II for x 0 =0.0 , 

h 2 /h 1 = 9.0 , R i /L= 9.0 and %LIE2=0.0\\%5 , x Q =t Q D l /h 2 1 , 

g 0 t =- © 0 /(l -Vj) (material pair B) 








H b l) / a o T ^i 



Figure 6-117: The normalized stress intensity factor k(b j) as a function of 
nondimensional time x for an edge crack in Model II for x 0 = 10.0 , 

* 2 /^ = 9. 0 , R i /L=9.0 and xLIE 2 = 0m\%5 , x Q =t Q D l /h 2 l , 

c 0 T =- 0 ^ E 1 © 0 /( 1 -Vj) (material pair B) 
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Figure 6-118: The influence of Tq on the normalized stress intensity factor k(b j) 
as a function of nondimensional tune t for an edge crack of 
length / 1 /^ 1 =0.2 in Model H , , RJL* 9.0 and 

XL/£ 2 = 0.01185 , t 0 =/ 0 D 1 //»5 , Oo^-a’j^Go/Cl-V!) (material 
pair B) 
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Figure 6-119: The influence of t 0 on the normalized stress intensity factor k(b^) 
as a function of nondimensional time x for a broken clad in Model 
II for Pj =0.552538 , A 2 //i 1 =9.0 , /?,/L= 9.0 and xL/£ 2 = 0.01185 , 

x Q =t 0 D 1 /h 2 l , c 0 r =- a’j £ : © 0 /(l -v^ (material pair B) 
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Figure 6-120: The normalized stress intensity factor k(aj) as a function of 
nondimensional time t for an under-clad crack in Model II for 
0^=0.4512416 , t 0 =0.0 , ^1^=9 . 0 , /? f /L=9.0 and xL/£ 2 = 0.01 185 , 

x 0 =t 0 D l /h] , a 0 r =- aj £ x © 0 /(l-v{) .(material pair B) 
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Figure 6-122: The normalized stress intensity factor k(a 2 ) as a function of 
nondimensional time % for an under-clad crack in Model II for 
02 = 0.4512416 , t 0 =10.0 , h 2 /h 1 =9J 0 , RJL= 9.0 and %L/£ 2 = 0.01185 , 

t 0 =t 0 DJh 2 , , c 0 T =- a\ E x © 0 /(l -Vj) .(material pair B) 
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Figure 6-123: The normalized stress intensity factor k(b 2 ) as a function of 
nondimensional time T for an under-clad crack in Model II for 
02 = 0.4512416 , ^=10.0 , A 2 /A 1 =9.0 , tf f /L= 9.0 and xL/E 2 =0.0US5 * 

x 0 =t 0 D ll^i > °o r=_ a 'i E l ©o /(! “ v i> -(material pair B) 
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Figure 6-124: The influence of t 0 on the normalized stress intensity factor k(a 2 ) 
as a function of nondimensional time T for an under-clad crack of 
length l 2 /h 1 = 1.0 in Model II, oc^O.4512416 tyh^.O , RJL= 9.0 

and xL /£ 2 = 0.01185 , , o 0 r =- a\ E { Q 0 K 1 -v^material 
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Figure 6-125: The influence of t 0 on the normalized stress intensity factor k(b 2 ) 
as a function of nondimensional time t for an under-clad crack of 
length / 2 //ii = 1.0 in Model II,oc2=0.4512416 h 2 /h l =9 . 0 , RJL= 9.0 

and = , 'i 0 =t 0 D l lh 2 1 , C 0 T =- a’ x FjGq/CI -VjXmaterial 

pair B) 
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nondimensional time x for an edge crack crossing the interface in 
Model II for p 1= 02 = 0 . 01 87223 ,x 0 =0.0 , V A l = 9.0 , « J /L= 9.0 and 

XL/E 2 = 0.01185 , x 0 =t 0 D 1 /h\ , O o r =-a’ 1 £ 1 0 o /(l-v 1 ) . (material 
pair B) 
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Figure 6-127: The normalized stress intensity factor k ^ as a function of 
nondimensional time T for an edge crack crossing the interface in 
Model II for p x = 0^=0.0187223 , x 0 =0.0 , V A l =9 - 0 , ^/L = 9.0 and 

XL/E 2 = 0.01185 , x 0 =f 0 £> 1 /Ai , c 0 r =- a’j © 0 / (1 -v x ) . (material 
pair B) 
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Figure 6-128: The normalized stress intensity factor k(b 2 ) as a function of 
nondimensional time x for an edge crack crossing the interface in 
Model II for =02 = 0.0187223 , x 0 = 0.0 , /i 2 //z 1 =9-0 , RJL= 9.0 and 

XL/£ 2 = 0.01185 , x 0 =r 0 D 1 //ii , o 0 r =- oi ^©o/Cl -v x ) . (material 
pair B) 
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Figure 6-129: The normalized stress intensity factor k x as a function of 
nondimensional time t for an edge crack crossing the interface in 
Model II for pj = 02 = 0.0187223 , >10=10.0 , A 2 //z 1 = 9.0 , RJL= 9.0 and 

XL/£ 2 =0.01185 , x 0 =r 0 D 1 /^ , c 0 r =- a\ E^qKI -Vj) . (material 
pair B) 


281 



in 

o 



.0 22.5 45.0 67.5 90.0 


r=tDJh\ 

Figure 6-130: The normalized stress intensity factor as a function of 
nondimensional time X for an edge crack crossing the interface in 
Model II for p 1 = oc 2 = 0.01 87223 , x 0 =10.0 , h 2 /h l *9.0 , RJL= 9.0 and 

XL!E 2 = 0.01185 , x : 0 =f 0 D 1 //»5 , o 0 r =- a\ E^/Q. -v{) . (material 
pair B) 
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Figure 6-131: The normalized stress intensity factor k{b 2 ) as a function of 
nondimensional time x for an edge crack crossing the interface in 
Model II for p 1 = 02 = 0.0187223 , x o =10.0 , f^/h^.O , Ri/L=9.0 and 

%LIE 2 = 0.01185 , 1 0 =/ 0 Z> 1 /*i , c 0 T =-a\ E.QoKl-^ (material 
pair B) 
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Figure 6 - 132 : The influence of % 0 on the normalized stress intensity factor k x as 
a function of nondimensional time x for an edge crack of length 
///ij = 1.5 crossing the interface, in Model II for 1^ = 02 = 0.01872238 

, /j 2 //i 1 =9.0 , R i fL=9.0 and xL/E 2 =0.01185 » T o='o I V ,, i > 

o 0 r =-a’ 1 £: 1 © 0 /(l-v 1 ) . (material pair B) 
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Figure 6-133: The influence of x 0 on the normalized stress intensity factor 

as a function of nondimensional time t for an edge crack of 
length ///«! = 1.5 crossing the interface , in Model II for 
p 1 = cx 2 =0.01872238 , f^/h^.O , RJL= 9.0 and %L/E 2 =0.01185 , 

'i Q =t 0 D l lh 2 l , c 0 T =- a\ E x © 0 /( 1 -Vj) . (material pair B) 
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Figure 6-134: The influence of Tq on the normalized stress intensity factor kip 2 ) 
as a function of nondimensional time x for an edge crack of 
length ///ij = 1.5 crossing the interface , in Model II for 
P 1 = a 2 = 0*01872238 , , JtyL* 9.0 and xL/£ 2 =0.01185 > 

XQ=t 0 D l /h^ , Gq T =- a\ E x 0 O /(1 -v x ) . (material pair B) 
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Appendix (A) 

The Governing Equations for the Displacements 

For a two dimensional, isotropic, elastic material in the absence of the 
body forces, we have the equilibrium equations 


do xx do 

+ ^-0, 
dx dy 

(a) 

do do 
dx + dy °’ 

(b) 

the strain displacement relations 


du 

e ~ = d~x’ 

(a) 

du 

£yy ~dj' 

(b) 

du du 

(c) 


( 1 ) 


and the strain-stress relations (Hooke’s law) 
for plane stress 


O 

II 

e 

II 

II 



(a) 


(b) 


( 3 ) 
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and for plane strain 


ii 

8* 

n 

II 

o 


1 “V 2 v 

e xx E (°xx J_ v a yy)’ 

(a) 

1 —v 2 v 

e yy e ^yy l-v 0xx ^' 

(b) 


( 4 ) 


By solving equations (1), (2), and (3) for the plane stress, and equations 
(1), (2), and (4) for the plane strain, we end up with the governing equations as 


follows 


(k-1)V 2 « + 2(^ + ^) = 0, 
dx 2 axoy 


(a) 


(k-1)V 2 i>+2(^ + ^) = 0, 
axdy dy 2 


(b) 


(5) 


where 

k = (3-4v) for plane strain. 
k = for generalized plane stress 

v = Poisson’s ratio 
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Appendix (B) 


Some useful integrals 


e~ px sin qxdx = — 

0 r+r 


f xe~ px sin qxdx = 

Jo (pW) 2 


f e~ px cos qxdx = P 

Jo p^+r 


r°° P 2 ~q 2 

xe px cosqxdx = — — — - 

Jo (pW) 2 


f" e-X— ) 

J-oo^+r 2 ) 2 2 r 

P e-^-^dp = ^[1-(S-X) ]e-**- x > 

J-CC (p 2 +^) 2 2 ' 

J—(p 2 +* 2 ) ' 

P 2 ~2~2 e~ tp(s ~ x)d P = 5 4 1 ' 1 • + (■*-*>] e_,(i ' x) 

J-o=(p 2 + r2) 2 2r3 

f °° 1 «“W * "*> dp = ^[ e ~* s ~ x) -1 ] 

J_oo/7 (p 2 + / 2 ) f 2 

P —P— e-^-^dp = -5 [i^] 

J-oo^+i 2 ) 2 2 * 


P f'Tr e~ ip(S ~ x) dP = ^+(5-*)]^-*)' 

J-oc^+r 2 ) 2 2 * 


; 


; s>.x 


; 5* > jc 


; j> * 


; $> x 


; s<x 


; set 
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1 


; s<x 



(p 2 +f-) 


1 

(p 2 +t ) 2 

1 

p(p 2 +f) 


-ip{s-x) ^ H e (s~x)t 

2 


e -ip(s-x) dp = 5 1 [ 1 -(J— Jc)r] 
2 fi 

e -ip(s-x) dp m 2EI| l- e (s-x)t] 

r 


; s<x 


; s< x 
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Appendix (C) 

The f unc tio ns that appear in the solution of crack problems given in 
Chapter 2. 

Kj + 1 

_ K j-i 
62 2~ 

K,-l 

e 3 = m { — - — + ha)e lha 

e 4 = me lha 

K i _1 

e 5 = 

k 2 -1 

«6 = -(-2— *«) 

K] + 1 

e n = -m ( — — + ha)e 2ha 
+ 1 

e 8 = m (— 2 Aa ) 

K, + l 

e 9 ^-{-=Y~-h a) 

. _ _2>ia 

e 10 e 

«12 = 

e 13 = -(K 1 + ha)e 2ha 
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^14 = — ( K i~ A°C) 

Cl 5 = ( K 2 -/ l a ) 

Mj =-(l-r 1 a)e w i a 

M 2 - (-t 1 a)e~h a 
M 3 = -m(t 1 -h)aeh a 
M4 = (t 2 -h)ae~( t 2- 2h ) a 
M 5 = -m[ 1 +(t l -h)a]e'i a 
M 6 = [l-(t 2 ~h)a]e~^2- 2h ) a 
k ,-1 

A / 7 = -[_ (r 1 -/i)a]^i a 

1 G >-1 

M g = [ — - — + (t 2 -h)a]e~ ( - t 2 ~ 2h ^ a 
M 9 = + (t 1 -h)a]eh a 

M l0 = {^--(t 2 -h)ale-^2-^ 


m = — 

M 2 


*1 = 

C 2 = 
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where 


- i^TT \ a + ;^TT <i> d h 

- \l 5 *»< *&>** 

®s ’ 5^1 J* , «>*i<'i>‘*i + i ^f T J‘ ! «io* 2 ('2)*2 


C 4 - 


C« = Jjij J a ’ »1, ♦,(», ><*1 *^r 1 j'Hu UV*2 




*1111 
2 ""2D 4 2D ’ 


8 


// 3 = { ^2 [ 1 + 2/»a (2r j a - 1 ) ] } e- ('i + a 
+ {-d 2 rf 3 [l+2(? 1 -/j)a]}A- 2 ^ a 
+ {d 1 d 3 [l + 2f 1 a] } e (»i -•**)« 

+ {-rf 1 d 3 }<r<'i +4/l > a 

// 4 = { - d 3 d 5 ( 2 ha) [l-2(t 2 -h)a]-d 1 d 5 }e~( t 2 +2h ) a 
+ {d 2 d 5 }e~ t 2 a 


1 11 K 1 1 

5 2 2D 3 2 D 1 


1 1 


K ll 




2D 


2 D 
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H 1 = {d 2 d A (l-2t i a)}e'h a 

+ { -^ 2^3 [ 1 -2(f 1 -/*)a] } e -(«i + 2A>« 

+ { d 2 d 3 (2/ta) [1 + 2(r j - h) a] d 4 } " 2/l > a 

+ {d 1 rf 3 }e( , i- 4/, > a 


// 8 = {-d 2 d 5 (2/ia)+d 4 d 5 [l-2(r 2 -/i)a] }e-' 2 a 
+ {-d 3 d 5 [l-2(f 2 -A)a] }f-(' 2 +2A ) a 






H n--r el + 7~n H i 


d^D 


Cl'i df. 1 

*, 2 = S [l-2(r 2 -A)c.]^-“>“+;p// 8 


and 

D = -d 2 rf 4 + [^ d A +d 2 d 3 +d 2 d 2 (2 to) 2 ] e-^-^d 3 <r 4Aa 

dj = mK 2 -K 1 

*f 2 = wik 2 + 1 
d 3 = m— 1 

<i 4 = m + Kj 
^5 = 1C 2 + 1 

4s = m ( K l + 1 ) 
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Appendix (D) 

If the function /(z) has a pole of order m at z = a , then the residue of j\z) at 
z = a , is given by 

Res(z = a)= lim - — — - {(z-a) m /(z)}] . (1) 

z ^ a {m-\)\ d/n-1 

In the case of a simple pole (m = 1) , equation (1) gives 

Res(z = a) = lim [{z-a)fiz)] . (2) 

z — > a 


If ,f(z) is expressed as the ratio 




m 

D(z)- 


(3) 


where N(a ) is finite and nonzero and D(z) vanishes at z = a in such a way that 


(z-a) 


approaches a finite limit as z-» a , then equation (2) can be written as 


Res{z = a) = lim [(z-a) 

z— > a 


m . 

D( z y 


(4) 


By using L’Hopital rule to evaluate the limits we have 

N(z) 


Res(z = a) = [ 


d . 
* D(Z) 




(5) 


Now, let equations (2.72) be written in the form 


( 6 ) 


where 
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©q e tz [cosh ^ h 2 cosh ^X’—^sinht^f^ ■ y,n ^ Cl x ’ 3 
z [ cosh C,j /ij cosh ^ t^+^sinh Ci s i n h ^2 ^2 3 


(a) 


©O e ft COS/l £ 2 (*’ - A 2 ) 

“ z [ cos/l COS/t ^ 2 ^2 + t 1 Ji'«A Cl A 1 ^ C 2 h 2 ] ’ 


( 7 ) 


(b) 


and 




( 8 ) 


2 

/j(z) , / 2 (z) have simple poles at z = 0 , and at z m = -D 1 w m , where ± w m are the 
roots of the equation 

cosw m h l cos5w m h 2 -r\ sinw m h l sinbw m h 2 = 0 (9) 


and 




**, 


D. 


( 10 ) 


2 

By applying equation (5) the residue at z = 0 and z ~ —D\ w m f° r both/^z) and/ 2 (z) 
are 


Res[z = 0 ;/ 1 (z)] = © 0 - 

(11) 

Res[z = 0 ;/ 2 (z)3 = © 0 - 

(12) 
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, Q 0 e~ ,D i w2 m 

Res[z = -D 1 w 2 m ;/j(2)]= * 

rv m 

T 

[cos W m X’ cos w m &h 2 +r\ sin w m x< sin w m 8 h 2 ] 

[(/i j +T|8/2 2 ) sin w m h l cosw m hh 2 + (hh 2 +r\h l )cosw m h l sin w m 8 h 2 \ ' 


, O 0 e ,D i w l 

Res[z = -D l w m ;/ 2 (z)]= * 

m 

~2 

cosw m b(x'-h 2 ) 

[(h j + T) 5 h 2 ) sin w m h l cosw m 6h 2 +(Sh 2 +T\h 1 ) cos w m h l sin w m hh 2 ] 


( 14 ) 
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Appendix (E) 


For z 0 = 0 


I 


1 ©!<**.*) . . ^ e~ x K 


©n 


-dx* = 1-2 X 


[ sin X m cosX m BR+r\ ( cosX m - 1 ) sin X m B R ] 

[ ( 1 +T) BR ) sin X m cos X m SR + ( bR +t| ) cos X m sin X BR ] 


f i+/? Q 2 (x*,i) ~ 

Jl ®o tx 


e TA rn 

^m 5 


sinX m BR 

m 


t (1 +T) 8/?) sin X m cos X BR + (BR +r^) cosX sin X m BR ] 


( 1 ) 


( 2 ) 


For t<t 0 


/•I ©! (AT* ,T) 

Jo ©0 


DO 





* 


[ sin X m cos X m B R + 1| (cos X m -l)sinX m BR] 

[ (1 +T| BR) sin X m cos X m BR + (BR +Ti) cosX m sin X m BR] ’ 


1 ©O ^0 ,4 

x oK 


sinX m BR 

m 

[(1 +r\BR)sin X m cos X m BR + (BR+r\) cos X m sin X m BR ] ‘ 


( 3 ) 


( 4 ) 


For x> t 0 
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,i e i (.x*,T) ^ = x _ 2 £ e^je^X-l) 

•*0 ©0 m = 1 _ 1,4 

T 0 K m 


* 


[ sin X m cosk m bR+r\ ( cos X m -1) sin X m bR) 

[ (1 +ri 6/2) sin cos + ( 87? +T|) cos X m sin X m 87? ] 


J 


1+jR ©2^ 


©r 


dx*=R - 2 £ 
m=l 


1) 

^6*1* 


( 5 ) 


sinX m bR 

[ (1 +T) 5R) sin X m cos X m 87? + (87? +T)) cos X m sin X m bR ] 


( 6 ) 


where 


„ *2 . ,0, *1,0 1 




it’, 


and X m are the roots of the following equation 
cos X m cosk m bR-r\ sin X m sin X m bR = 0 


(7) 


( 8 ) 
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Appendix (F) 

Further definitions to support the solution of the crack problems given in 

Chapter 2. e i ,(/=l ,15) , Mj(j=l ,10) are found in Appendix (C). 

+ 1 

e 16 = (— + 

e n = ( K 2+^itx)^ 2;i i a 

eis = -^ 

K~ + l 

«19 = + 

k 2 +1 

£, 20 = (~2 

K,-l 

« 21 = (A-+aL)e 2i “ 
k 2 -1 

e 22 = (— 5 aL ) 

e 23 = ctLe 2 ^ 
e 1A = a L 

e 25 = -(-^Y~ +/j i a ) g2/ ‘ 1 “ 

M n =-[l+(f 2 — L)a]e'2 a 

M 12 = “ _L ) a + P < ”3 (ifc-I) a}] eh a 
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c 2 ° spfi J„‘ Q} + ^T J^ 24 ' 1 ' 2 ( ' 2,<i ' 2 

i3 = ^ttJ 0i ' 25<, » ( ' 1> ‘* 1 + ^TT j^e*W2>*2 

Q" * 1 f ' 27<>l^l)^l + | 2^2 fr 2> df 2 

K l + iJ a l k 2 ' riJ a 2 

,3 “^J/ C9+l( ' l) * 1+ ^J/ £!l0 * 2< * ) ‘ i ' i 

I a i *• u 2 

c 6 - 7^rrr , en*i('i)‘ i 'i + ^rTj 6! ei2 | i'2<'2)*2 

K l +AJ a 1 k 2 + lJ a 2 





b zQl 6 


where 

Ci = 5(l-2ria)e- , i a -ye 3 +ie 7 
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g 3 = 5+ 1 J- { (2d 3 d 5 ) (2h l a) fi’ 13 + (~d 5 d 2 )Q i5 
+ [(d l d 5 )+ (d 3 d 5 ) (2/i ! a) (k 2 + 2/ij a) ] g’ 15 } 
q a = ^ + 1 i- { (2d 3 d 5 ) (2/i i a) fi’ 14 - + (-< d 5 t d 2 ) Q 16 
+ [(d l d 5 ) + (d 3 d 5 ) (2 h j a) (k 2 + 2/i x a) ] g’ 16 } 

e 5 = 5 £ "'“-5 e 3 + y e ’ 


1 K 1 

26 = -2 e 4 + y28 

g 7 = ^ + li- {(- 2^5 d 4 ) <2 13 + (^5 d 3) 2 13 
+ ( (^2 d 5) ( 2/l l a > + d 5X K 2 + 2,I 1 a)] ^15 

+ [ (d 3 d$) (i^ + 2/ij ct) ] 0 is ^ 



2s 'ir + 55 -((-aweM+G-wei 4 

+ [ (d 2 ^5) (2A 1 oc) + ( - ^5)^2 + 2/i 1 a) ] g 16 

+ [ ( d % ^5)0^2 + 2 /zj a) ] 2 15} 

e 9 =( ~^ ) e 3 ^ ia + 5 ( ^" 2Ai£X)ei1 


ei0 = (_ 2^ )e ^ 

+^( K 2 - 2 / l l a ) 2 i 2 


+ < 4 ) ^ Ci6£ 


2A, a 


2d. 


3 . 1 


a a ~^*^*( r - s ^ a a 


e 2 h i a 


2h, a 


+ (-3-)( K 2 + 2/j l a )p-2l5 e 1 


2i2 = (T> [1 " 2(r 2- /, i) a ] e ' <,2 ' 2,tia > +( r )e8 
“2 2 

d 2 D 0 d 2 L ) 0 
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2i3 ~ (<ftSi+q 9 s 4 -p 4 r 2 -p 5 r4)e~( 4L+t \) a 
+ (q n s 2 +q 9 s 3 -p 4 r 1 -p 5 r 3 )e~( 4L -h) a 
+ (q 7 s 3 +q s s 2 -p 4 r 3 -p 6 r{) e ~ {4L + : “ 'i) “ 
+ (?7 ^ + ^ i ~P 4 >4 -p 6 r 2 ) e ~ (4i + 2Aj + /j) a 

+ (q 9 S 1 -p 5 r 2 ) e ~ < 4L - 2h i + 1 1 ) a 
+ ( q g s 2 -p 5 r l )e~^ L ~ 2h \- t O a 
+ (<7g s 4 -p 6 r 4 )e-( 4L+4h \ +t J a 
+ (<?8 J 3 ~Pk r $) e~ i - 4L+4h \~ t \> (x 
+ (?10 s l-P7 r 2) e ~ (2L + ' i)a 
+ (<h0 s 2-P7 r l) e ~ (2L ~ ! ' )a 
+ (<hO s 3 + <7ll s 2-p7 r 3-P8 r l) e ~ (2L + 2hl ~ ti)a 
+ (<7lO s 4 + 4n s l~P7 r 4 ~P% r 2 )e~( 2L+2h \ + t 1 ) a 
+ (<hl s 3 + c ll2 s 2 ~P% r 3 -P9 r 1 )e-( 2L+4h i->O a 
+ (.<h 1 ^4+^12 J i -Ps r 4 -Pg r 2> e - (2L+4h \ + *i> a 
+ (^12 s 3 ~Pg r 3 )e-Q L+6h \- t \> a 
+ 12 s 4 ~Pg r 4 ) e ~ {2L + 6h l +t J a 
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2 14 = (07 5 5 + 09 ^6 ~Pa r 5 ~Ps r 6> e ~ (4i " ' z) ° 

+ (q-j s 9 +q g s g +q g s lQ -p A r g -p 5 r w -p 6 r g ) e~ + ' 2 > a 

+ (07 s 6 + 08 s S + 09 ^7 “P4 r 6 “^5 r 7 ~/>6 r s) £~( 4L + 2 V r 2> a 
+ (<? 7 s g + q 9 s g -p 4 r g -p 5 r 9 ) e" (41 - 2*! + / 2 ) a 

+ (07 J 10 + 08 5 9 “/>4 r \0~Pt r 9> e ~ (4L + 2A l + r 2> a 
+ (09 ~/>5 ^5) e~ (4L “ 2fc i “ * 2 ) “ 

+ (<? 7 s 7 + <? 8 j 6 -/? 4 r 7 ~ p 6 r 6 ) e ~ (4L + “ * 2) a 

+ (08 s 10-/>6 r 10> e " (4L + 4/, i + V a 
+ (09 5 8 “P5 ^g) (4i " 4 *i + r 2> a 

+ (08^7-P6 r 7) e_(4L + 6/t '" ,2)a 
+ (01O‘ S '5 _ / , 7 r 5) <?_(2i_ ' 2)a 
+ (010 ^9 + 01 1 5 8 ~Pl r 9 ~P% e ~ (U + ^ “ 

+ (010 s 8 “P7 r 8> e_ (2i _2Al +/j) “ 

+ (010 s 10 + 011 J 9 + 012 s 8 -Pi r \Q~P% r 9 -P9 r %) e ~' i2L + 2k ' + ‘ 2) 0 
+ (010 s 6 + 01 1 s 5 ~Pl r 6 ~P8 r 5> e~^ L + lh \~ h ) a 
+ (010 J 7 + 01 1 J 6 + 012 ^5 -Pi r l —P% r 6 -P9 r 5> e _(2L + 4/! i~ ‘2 ] a 

+ (01 1 *10 + 012 *9 ~PS r 10~P9 r 9> (2L + + f 2 ) a 

+ {q x J s 7 + q n s 6 -p s r 7 -p 9 r 6 ) e^ 21 + 6h \ ~9 a 

+ (012 J lO-P9 r lO>«“^ + 6 * 1 + < * )O 
+ {q\2 s i -P9 r i ) e ~ ( - 2L+ * h \ ~‘i> a 
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Q 15 = (Pi r l +P3 r 3 - <h s 2 - <73 s 3> e ~ (AL ~ h) “ 
+ (p 1 r 2 +p 3 r 4 -q 1 s 1 -q 3 s A )e~ ( - 4L+t J a 
+ (p l r 3 +p 2 r l -q 1 s 3 -q 2 s 2 )e- ( - 4L+2h i- t \'> a 
+ (Pir4+p 2 r 2 -q 1 s4-q 2 s 1 )e~( 4L ' i2h \ + , O a 
+ (P2 r 3~ ^2 s 3) e ~ +4h \~ t O a 
+ (p 2 r 4 “ ^2 ^ e_ (4i + 4Al + ^ “ 

+ (Pir l -q 3 s 2 ) e~ ( 4L ~ 2h i ~ *i) a 
+ (p 3 r 2 -q 3 s l )e~ ( - 4L ~ 2h \ Jrt -i> a 

+ (-? 4 j 1 )*-( 2L+ ' 1 )“ 

+ ( _l ?4^2) e_(2Z: ' _,l)a 
+ (~?4 ^3 ~?5 5 2 ) + lh \ _, i) a 

+ (~?4 ^4 ~<?5 «~ (2L + 2/i i + 'i> a 

+ (-<? 5 5 3 -9 6 S2) e ~^ 2£ ' + 4,Il-, ^ a 

+ (- <?5 s 4 -< ?6 s 1 )e-( 2L+4 V'i)« 

+ ( - q 6 s 3 ) e~ + 6h i ~ f 0 a 

+ (-q$ s 4 ) e~( 2L+6h \ + t O a 


310 



016 = (Pl r 5 + P3 r 6-<h s S-43 s 6) e ~ (4L - t2)a 
+ (Pi r 9 +P 2 r % +P3 r 10-<?1 s 9 _< ?2 ^8 “93 5 10> «" (4L + ' 2> ° 

+ (Pl''6 + P2 r 5 + / , 3 r 7-9l s 6-92 5 5-93^7) e_(4L+:y,1_ ' 2)a 
+ (P 1 ^ 7 +/ , 2 r 6 “ 9 1 « 7 “ 9 2 ^6> e_ (4L + ^ 1 _ <2> “ 

+ (Pi 'g +P 3 r 9 - 9i ^- 93 ^ 9 ) <r (4L-2>, i + ^ ° 

+ (Pl ^10 + P2 r 9~9i ^lO-92 ^ 9) e “ (4 ' L+2,,1 + ' 2)a 

+ (p 2 r 7 - 92 ^7) e~( 4 L+ 6 h i~ ‘i> a 
+ (P2 r 10-92*10) e-< 4L+4fc i + ' 2 ) a 
+ (P3 r 5 “ 93 $5) eT (4L - 2h i~‘2> a 

+ (p 3 r 8 “ 93 •%) ^ - 4 *i + ' 2 ) a 

+ (-^ 4 ^ 5 )(? _(2L_, 2) a 

+ (- 9 4 s 9 _ 95 s 8 ) e ' (U + r2)a 
+ (-^4- s 6~95' y 5) g_(2£ ' +2,Il-,2)a 

+(-9 4 5 8 )e- (2i - 2 V /2)a 
+ (-^4 5 7-95 s 6-96 s 5) g_(2L+4,Il ' ,2)a 

+ (~<l4 s 10-‘l5 s 9-‘l6 s &') e ~ (2L+2hl + ,2)a 
+ (-q 5 s 1 -q 6 s 6 )e-( 2 L+ 6 h i- t 2 )a 

+ (“95 J 10 “96 f 9 > (U + 4/l ' + “ 

+ (- 96 ^ 10 )^ (2L + 6Al + ,2)a 
+ ( - q 6 s 7 ) e~ ^ + 8 * 1 “ *2) a 
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D 0 = (Pi < ?7 + - c, 2^9 + / , 3^8“ < ?l/ ? 4“ < ?2^5 _ ^3/ , 6) e_4La 
+ (Pi ?8 +P2 <h-<hP6~<l2PA) e~ 2(2L + h J a 
+ (pj <? 9 +p 3 q-]—qiP$— q$p A ) e~ 2(2L ~*i> a 

+ (Pi ^10 +P3 ?1 1 - P7 “ <?3P8 “ <74^4 - <?5P5) ^ 
+ (Pi <?1 1 +P2 <?10 + P3 <?12 “ $lP8 -< ?2P7 

-^3P9-?4P6"^5P4-^6P5) g_2(Z ' +/!l)a 
+ (Pi ?i2- +P2 <7i 1 ~<hP9 ~ <hPi - qtfs ~ q(,p£ e ~ 2( - L + 
+ (P2<ls- ( l2P6) e ~* (L + hl)a 

+ (p 3 q 9 -q 3 P 5 )e-^ L - h )> a 

+(p 2 qn-q 2 P9- P^Pt) e ~ 2(1 +3Al)a 

+ (P 3 4 10 ■ - ^3 Pi - ?4P5> e ~ 2(L h ‘ )a 

+ (-?4P7> 

+(-q4P 8 ~q5P7) e ~ 2k ' a 

+ (~q4P 9 ~q 5 p 8 -q 6 p n ) e~ 4h \ a 

+ (-q5P9~<l6Ps') e ~ 6hia 

+ (-q 6 P 9 ') e ~ 8hia 

Q\* = Qx 3 e- 2h ' a 

Qu = Qu‘- 2h ' a 

Q’l5 = Ql5 e ~ 2h ' a 

Q\6 = Qi6 e ~ 2h ' a 

<?1 = (^-1){ (-2rf 6 d 5 d 4 ) + (-2d 3 rf 1 rf 4 ) + (-2rf3rf 2 rf3) 
+ (-2J 3 rf 2 ^ 3 )(2/i 1 a] 2 } 
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q 2 = ty-l ){(2d 6 d 5 d 3 )+(2d 3 d l d 3 ) } 


q 3 = {d 7 ~l){2d 3 d 2 d 4 } 

<?4 = { - 2 ^ 2 ^ 2^4 ) 

q 5 = { (2^ d 2 djtlh^j 1 + (2^ <* r <* 4 ) + (ld 2 d 2 <* 3 ) } 

= i — 2<*2 <*j <*3 } 

q 2 = (<*7 — l){ (<*g<*2<*5)(2/? 1 a) + (K2 + 2/i 1 a)[(— <* 6 <* 4 <*j) 

"t” ( <*3 <* 2 <* 4 ) + ( — <*3 <*2 <*3) ( — <*3 <*2 <*3X2/1 j Oc) 2 ] } 

7 g = (<*7— l)(K 2 + 2/j2a) { (<*^<*3<*5) + (<*3<*2<*3) } 

< 7 9 = (<*7 - 1 )(Kj + 2 h L a) { <*3 <*2 <* 4 } 
g l0 = 1 <* 7 ~ (k 2 + 2(XL)} 1 <* 2 <* 2 <* 4 } 

<*11 = 1^-( k 2 + 2 cx L )} 1 (-<* 2 <* 2 <* 3 )( 2 A 2 a ) 2 +(-<* 2 <* i <* 4 ) 

+ (— <* 2 <*2<*3) } 

<7 12 = 1^7” (K2 + 2aL)} { <*2<*2<*3 } 

/•2 = (<*7-l)l<* 6 <* 2 <* 4 } 

r 2 = (<* 7 - 1 ){ <* 6 <* 2 <* 4 } ( - 1 + 2 r : a ) 
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r 3 = fciy-lM -^6 d 2 d 3Hl+(2/i 1 a)(l + 2 f 1 a)} 
r 4 = (d 7 - 1 ){ d 6 d 2 d 3 } { 1 + 2(/ij - t^ct) 
r S = 

r 6 = (ft- l)(-d 2 ){(d 2 d 3 )( 2 h 1 a) 2 +d 1 d A +d 2 d 3 } 

r l ~ l){<*2<*j<f3} 

r g = (<*7-l){ -d 3 d 2 d A }{ -1 +2(f 2 -A 1 )a } 

r 9 = - 1 ) { [ - 1 + 2(^ - ) a ] [(<* 2 <*3 <* 3 )(2A ja) 2 +(d x d A d 3 ) 

+ (<*2 <*3 ^3) + (d A d$ <i^)] + (<f 2 <*5 dfi)(2h^Cl) } 
r 10 = (<*7 - 1){ [- 1 + 2 (f 2 - *!>«] [(-<f 3 d x d 3 ) + (-<* 3 d 6 d 5 ) ] } 


Pi = {(4 6 <* 2 )(-2<* 3 <* 5 )(2A 1 a)+(-<* 3 )[-<*7+2CL-A 1 )a] 

[(2d s d A d 6 ) + (2d 2 d 3 d 3 )(2h x a) 2 + (2 d x d 4 d 3 ) + (2rf 2 d 3 <f 3 )] } 

p 2 = { (2<* 3 d 3 )[ —dq + 2(L — Aj) a] [(.d^+id^)]} 


p 3 — { (2d 3 d 3 d 2 d^)[ d-j 4* 2(Z- /z } 
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Pa = {(-^ 6 d 2 d 3 <i 5 )(2/i 1 a)(K 2 +2A 1 a)-(d 6 d 2 d 1 rf 5 ) 

+ (<^2 d 2 ) [(d 2 d 3 )(2h | a) 2 d^ d^ A~d 2 d 2) ] 

+ [ — dq + 2 (L — hy) oc] [(dfi d ^ ^2^5X^^i^0 ("”^ 6^3 ^4 ^5X^2 ^/^oc) 
+ (-d 3 d 3 d 2 d 3 )(2A ja) 2 ^ + 2 /j L a) 

+(-d 3 d 2 )(x 2 + 2h 1 a)(d 1 d 4 +d 2 d 3 )] } 


P5 “ { (^6 ^2 ^2 ^5) _ (^2 ^2 ^2 ^4) 

+ ^3 ^2 ^4X^2 "^" jOC) [ dq + 2(L “ h j) oc] } 


p§ — { ( “ d<2 d 2 ^3) 0^2 1 ®0L — ^7 2(jL — /Z|) oc] 

[(^6 ^3 ^3 ^ 5 ) + (^3 ^3 ^1 ^ 3 )] 1 


Pq — {(-^2 ^3 ^2 ^ 4 X^ 7 + 1)} 


/?g — { (^2 ^3X^7 "^* 1)[(^2 ^3X^2 CC)^ "^” ^4 ^2 ^3 1 ) 


Pg “ {(”^2^3^! ^3X^7 + 1)} 


tJj — { (^2 ^ ^2 ^3)(^i®)( ^ 1 "t* Ot) *f (^ ^2 ^3) 

+ [-d 1 + 2(L-h 1 ) a ](<f 3 d 6 d 2 d 4 )(- 1 + 2fjO)} 


^2 — { ( 4^2 ^6 ^2 ^3)^ "*"2(*i — 

+ [-^+2(L-/j 1 )a](d 3 d 6 d 2 d 4 )} 

s 3 = {(rf 2 d 6 rf 1 £? 3 )(l +2?jCc) 

+ [-d 7 +2(L-* 1 )a](-d 3 d 6 d 2 d 3 )[l+(2A 1 a)(l+2r 1 a)]} 
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**4 " {{~"d 2 d^d^d^) + [ — dq-¥2{L— h^)<x] 
(d 3 d 6 d 2 d 3 )[l-2(t l -h 1 )a)} 


■S 5 ” { \dq + 1 + 2(t 2 '~L)(X \(— ^2 ^3 ^2 ^ 4 ) 1 

•5*5 = { \d~j + 1 + 2(/2 — Z^)oc] d 3 )[(d 2 d^ilh^a) 1 +d^d^+ d 2 d 3 ] } 

s*j = { [d~j + 1 + 2(t 2 — jL)oc] ( *“ d 2 d ^ d^ d^) } 


S 8 - { ( - d 2 d 2 d 2 ^4) + (^2 ^6 ^2 ^5) 

+ [ ~ d 2 + 2(L~ A^)oc] ( d 3 d 3 d 2 d^) [ * 1 2(t 2 } 

jp = { (d 2 d 2 )[d 2 d 3 (2h jCx)^ + d-y d 4 + d 2 d 3 ] 

+ (d 2 d 6 m 3 d 5 )(2h 1 a)(- l+2(t 2 -h l )a)-d l d 5 ] 

+ [-d 1 +2(L-h 1 )a](d 3 d 3 )[-l+2(t 2 -h 1 )a] 

[d 2 d 3 (2h jCc)^ +Jj + d 2 d 3 ] 

+ [-d 7 + 2(L-A 1 )a](-rf 3 ^ 6 ) 

[( -d 2 d 5 )(2h l a)-d 4 d 5 (- 1 + 2(/ 2 -/^a)] } 

■^10 = {(“ d 2 d 2 d^d 3 ) 

+ [-d 7 +2(L-h l )a][-l+ 2 (t 2 -h x )a] 

[( — d 3 d 3 d^d 3 ) + (— d 3 d^d 3 d 3 y ] } 


where H 3 ,H 4 ,// 7 , // 8 , Z? , d l ,d 2 ,d 3 ,d 4 ,d 5 ,d 6 can be found in Appendix (C), and 

K, + l 
*7 = P — 
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Appendix (G) 

In general, the singular integral of the form 

f +l Mrm dr 

j_l r-s 


( 1 ) 


with 


w(r) = (l+r)- a (l-r)-P. 

(2) 

n 

'C' 

(3) 


can be obtained by using a Gaussian Quadrature formula. Let us rewrite 
equation (1) in the following form. 


/W-/W 


w(s) 


J _1 r-s J_i r-s J _1 r— 


(4) 


The first integral in the right hand side of (4) is bounded as r-s , i.e. 


>+1 fir) -As) 


w(s) 


r +l w(r) 

lim w(r)[ — ]dr 

r-s 

= | Mr) [f(s) -/(5)w(5){-|^}'] dr 


(5) 


where (') , is 

dr 

Equation (4) may now be written as 

w(s ) 

f +1 w(r)^r) = r l w(r)[ ^fl ]dr+w ( s )As)log\^- (6) 

J_! r-s J_i r-s 1+5 

Since the weight function w(r) is in the form of equation (2), then by using 
Jacobi-Gauss Quadrature formula (see, Erdogan [40] or Stroud [44]), we may 
have 
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( 7 ) 


.+1 » 

w(r)g(r)dr - 
-l /=l 


where r, is the I th zero of the Jacobi polynomial 
p(-P--“) (r/ ) = 0;/=l,2, n. 

and the weighting functions A { are 

„ ( 2 /i-a-p+ 2 )r(n-a+l)r(/i-p+l), 

' (n-l)!(n-a-p+l)r(n-a-p+l) 


( 8 ) 


2 -a-p 


_£»(- P.-a) 

dr n 


(r,)P 


(-p.-o) 

n+1 


('/) 


( 9 ) 


In our case g(r) is in the form 

Ms) 


m-As) 


g(r) = ■ 


Mr) 


r-s 


or as r — » 5 , g(r) will be equal to 

200 = f(s) -AsMs) {^}'. 


( 10 ) 


( 11 ) 


In the special case, if a=P=-^, the closed form solution of the integral (1) can be 
put in the following form 

r 1 *- sv- «> 

-1 ^ll-r^fr-s) 1=0 


where 
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(13) 


- r( V } 

d- = yin — , n-i = odd, 

d i = 0 , n-i = even . 

Also, the function g(r ) in equation (7) may be in the form 

g(r)=Ar). (14) 


or 

g(r)=fl.r)K ij (s,r) ,i,j= 1,2. (15) 

where /(r) is given in equation (3), and (s , r) , (j,y = l ,2) are the kernels of the 
singular integral equations (4.4). 

The integral as in equations (2.64), i.e. 

J G i j(x,t,a)da,i,j=l,2. (16) 


may be evaluated by using Laguerre-Gauss Quadrature formula, by changing 
the integral (16) into 

j~G ij (x,t,a)da = j~[G i j(x,t,a)e a ]e- a dct,iJ= 1,2. (17) 

then, the integral becomes (see, Stroud [44]), 

f [G i :(x,t,a)e a ]e~ a da - ^A l [G ij (x,t,a l )e CL i) 

J 0 /=1 


(ij =1.2). (18) 

where is the I th zero of Laguerre polynomial 

£„(«/) = 0 , /=1 ,2 n. (19) 


and the weighting function A l given by 
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, 1=1,2, n. 
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